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ABSTRACT 

We present multi-wavelength observations and modeling of Gamma-ray Bursts (GRBs) that exhibit a simul¬ 
taneous re-brightening in their X-ray and optical light curves, and are also detected at radio wavelengths. We 
show that the re-brightening episodes can be modeled by injection of energy into the blastwave and that in all 
cases the energy injection rate falls within the theoretical bounds expected for a distribution of energy with 
ejecta Lorentz factor. Our measured values of the circumburst density, jet opening angle, and beaming cor¬ 
rected kinetic energy are consistent with the distribution of these parameters for long-duration GRBs at both 
z ~ 1 and z > 6, suggesting that the jet launching mechanism and environment of these events are similar to 
that of GRBs that do not have bumps in their light curves. However, events exhibiting re-brightening episodes 
have lower radiative efficiencies than average, suggesting that a majority of the kinetic energy of the outflow is 
carried by slow-moving ejecta, which is further supported by steep measured distributions of the ejecta energy 
as a function of Lorentz factor. We do not find evidence for reverse shocks over the energy injection period, 
implying that the onset of energy injection is a gentle process. We further show that GRBs exhibiting simul¬ 
taneous X-ray and optical re-brightenings are likely the tail of a distribution of events with varying rates of 
energy injection, forming the most extreme events in their class. Euture X-ray observations of GRB afterglows 
with Swift and its successors will thus likely discover several more such events, while radio follow-up and 
multi-wavelength modeling of similar events will unveil the role of energy injection in GRB afterglows. 
Subject headings: gamma-Tay hurst: general - gamma-ray burst: individual (GRB 100418A, GRB 100901 A, 
GRB 120326A, GRB 120404A) 


1. INTRODUCTION 

Gamma-ray bursts (GRBs) have traditionally been modeled 
as point explosions that inject ~ 10^* erg of energy into a col¬ 
limated, relativistically-expanding hreball over a period of a 
few seconds. In this model, the subsequent afterglow radia¬ 
tion is synchrotron emission produced by the interaction of 
the relativistic ejecta with the circumburst medium. Depend¬ 
ing on the density prohle of the ambient medium, usually as¬ 
sumed to be either uniform (TSM-like’) or falling with radius 
as r~^ (‘wind-like’), this model has several verihable predic¬ 
tions: smooth light curves at all frequencies from the X-rays 
to the radio, which rise and fall as the peak of the spec¬ 
tral energy distribution evolves through the observer band; 
a ‘jet break’ as the expanding ejecta decelerate and begin 
to spread sideways; and an eventual transition to the sub- 
relativistic regime where the ejecta become quasi-spherical. 
In this framework, the energetics of the explosion and the 
properties of the environment can be determined from htting 
light curves with the synchrotron model, while a measurement 
of the jet break allows for a determination of the angle of col- 
limation of the outflow and the calculation of geometric cor¬ 
rections to the inferred energy. 

Despite its simplicity, this model was quite successful in the 
study of a large number of GRB afterglows (e.g. Berger et al. 
2000; Panaitescu & Kumar 2001, 2002; Yost et al. 2003) un- 
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til the launch of Swift in 2004 (Gehrels et al. 2004). With 
its rapid-response X-ray and UV/optical afterglow measure¬ 
ments, Swift has revolutionized the study of GRB afterglows. 
Rapidly-available localizations have allowed detailed ground- 
based follow-up. X-ray observations from Swift have been 
even more revolutionary, revealing light curves with multiple 
breaks falling into a ‘canonical’ series, consisting of a steep 
decay, plateau, and normal decay, sometimes with evidence 
for jet-breaks. While the steep decay has been associated with 
the prompt emission (Tagliaferri et al. 2005; O’Brien etal. 
2006) and the normal and post jet-break decay phases are as¬ 
sociated with the afterglow, the plateaus cannot be explained 
by the standard model (Zhang et al. 2006; however, see re¬ 
cent numerical calculations by Duffell & MacEadyen 2014, 
which suggest that the plateaus may be a natural consequence 
of a coasting phase in the jet dynamics between 10*^ to 
10'® cm from the progenitor). In addition, short (Af/f ^ 1) 
flares that rise and decline rapidly and exhibit large flux 
variations (AF/F ^ 1; Chincarini et al. 2010; Margutti et al. 
2010a) are often seen superposed on the light curves. These 
flares are believed to be more closely associated with the 
GRB prompt emission than with the afterglow, and have 
been interpreted as late-time activity by the central engine 
(Ealcone et al. 2006; Romano et al. 2006; Wang et al. 2006; 
Pagani et al. 2006; Liang et al. 2006; Lazzati & Perna 2007; 
Margutti et al. 2010a; Guidorzi et al. 2015). Some GRBs also 
exhibit optical flares, which do not always correspond to flares 
in the X-ray light curves (Li et al. 2012). 

These new features of Swift X-ray light curves cannot 
be simply explained in the traditional picture. New physi¬ 
cal mechanisms such as energy injection, circumburst den¬ 
sity enhancements, structured jets, viewing angle effects, 
varying microphysical parameters, and gravitational micro- 
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lensing have been invoked to explain various features of Swift 
X-ray light curves (Zhang et al. 2006; Nousek et al. 2006; 
Panaitescu et al. 2006; Toma et al. 2006; Eichler & Granot 
2006; Granot et al. 2006; Jin et al. 2007; Shao & Dai 2007; 
Kong et al. 2010; Duffell & MacFadyen 2014; Uhm &. Zhang 
2014). However, although a wealth of information is avail¬ 
able from X-ray light curves in general, definitive statements 
on the physical origin of these features requires synergy with 
observations at other wavelengths. Ultraviolet (UV), optical, 
near infra-red (NIR), millimeter, and radio data probe distinct 
parts of the afterglow spectral energy distribution (SED), and 
the various physical mechanisms are expected to influence 
light curves in these bands differently. Thus, a detailed anal¬ 
ysis of GRB afterglows requires multi-wavelength data and 
modeling. 

Of the GRBs with plateaus in their X-ray light curves, there 
is a small class of peculiar events that additionally exhibit 
an X-ray re-brightening of a non-flaring origin (AT /T ^ 
1 ), and an even smaller class where the re-brightening ap¬ 
pears to occur simultaneously in both the optical and X-rays 
(Mangano et al. 2007; Li et al. 2012; Panaitescu et al. 2013). 
An exemplar of this latter class is GRB 120326A, which ex¬ 
hibits a peak in its well-sampled X-ray light curve at around 
0.3 d together with a simultaneous optical re-brightening. 
Urata et al. (2014) reported optical and millimeter observa¬ 
tions of the afterglow of GRB 120326A, and invoked syn¬ 
chrotron self inverse-Compton radiation from a reverse shock 
to explain the millimeter, optical, and X-ray light curves. 
By fitting energy-injection models in a wind-like circumburst 
environment to the X-ray and optical R-band light curves, 
Hou et al. (2014) proposed that a newborn millisecond pul¬ 
sar with a strong wind was responsible for the re-brightening. 
Melandri et al. (2014) present multi-band optical and NIR 
light curves of this event, and explore various physical sce¬ 
narios for the re-brightening, including the onset of the after¬ 
glow, passage of a synchrotron break frequency through the 
observing band, and geometrical effects. 

Here we report detailed radio observations of this event 
spanning 4 - 220 GHz and 0.3 to 120 d, making this the 
first GRB with an achromatic re-brightening and with such 
a rich multi-band data set. We perform the first broad¬ 
band modeling for this event using a physical GRB after¬ 
glow model (Granot & Sari 2002) using methods described 
in Laskar etal. (2013) and Laskar etal. (2014), employing 
Markov Chain Monte Carlo procedures to characterize the 
blastwave shock. Our radio observations allow us to constrain 
the synchrotron self-absorption frequency, and to unambigu¬ 
ously locate the synchrotron peak frequency, together placing 
strong constraints on the nature of the X-ray/UV/optical re¬ 
brightening. We And the clear signature of a jet break at all 
wavelengths from the radio through the X-rays, allowing us 
to constrain the true energetics of this event. 

We next consider and test various physical processes that 
may cause a re-brightening in the afterglow light curve, 
and argue that energy injection is the most plausible mech¬ 
anism. We model the re-brightening as a power-law in¬ 
crease in blastwave energy, self-consistently accounting for 
the change in the synchrotron spectrum over the injection pe¬ 
riod, and compute the fractional increase in energy during the 
re-brightening. We interpret the energy injection process in 
the context of a distribution of ejecta Lorentz factors, and pro¬ 
vide a measurement of the power law index of the Lorentz 
factor distribution. 

To place our results in context, we search the Swift X- 


ray light curve archive for all events exhibiting a similar re¬ 
brightening, and present full multi-wavelength analyses, com¬ 
plete with deduced correlations between the physical param¬ 
eters, for all events with radio detections that exhibit simul¬ 
taneous optical and X-ray re-brightenings. This selection 
yields three additional events; GRBs 100418A, 100901 A, and 
120404A. We collect, analyze, and report all X-ray and UV 
data from the X-ray telescope (XRT) and UV-optical tele¬ 
scope (UVOT) on board Swift for these three events in ad¬ 
dition to GRB 120326A. We also analyze and report previ¬ 
ously un-published archival radio observations from the Very 
Large Array (VLA), Submillimeter Array (SMA), and West- 
erbork Synthesis Radio Telescope (WSRT) for events in our 
sample, and the first complete multi-wavelength model fits 
for GRBs 100418A and 100901A. Finally, we compare our 
results to modeling efforts of a sample of GRBs ranging from 
z ^ 1 to z > 6 , as well as with a complete sample of plateaus 
in SwiftlXKY light curves, and thereby assess the ubiquity of 
the energy injection phenomenon. We infer the fractional in¬ 
crease in blastwave energy over the plateau phase using sim¬ 
ple assumptions on the afterglow properties, and determine 
the unique characteristics of these events that result in multi¬ 
band re-brightenings. We conclude with a discussion of the 
results from this X-ray-only analysis and our full broad-band 
modeling in the context of energy injection in GRBs. 

Throughout the paper, we use the following values for 
cosmological parameters: Dm = 0.27, Ha = 0.73, and Ho = 
71kms“'Mpc“^ All times are in the observer frame, uncer¬ 
tainties are at the 68 % confidence level (Icr), and magnitudes 
are in the AB system and are not corrected for galactic extinc¬ 
tion, unless stated otherwise. 

2 . GRB PROPERTIES AND OBSERVATIONS 

GRB 120326A was discovered by the Swift Burst Alert 
Telescope (BAT, Barthelmy et al. 2005) on 2012 March 26 at 
01:20:29UT (Siegel et al. 2012). The burst duration is Tgo = 
26.7 ±0.4s, with a fluence of = (1.1 ±0.1) x 10“® ergcm“^ 
(15-150keV Barthelmy et al. 2012). A bright X-ray and 
UV/optical afterglow was detected by Swift (Siegel et al. 
2012; Kuin et al. 2012) and numerous ground-based obser¬ 
vatories. Spectroscopic observations at the 10.4 m Gran 
Telescope Canarias (GTC) provided a redshift of z = 1.798 
(Tello et al. 2012). 

The burst also triggered the Fermi Gamma-ray Burst Mon¬ 
itor (GBM) at 01:20:31.51 UT (Collazzi 2012). The burst 
duration as observed by GBM is Tgo = 11.8 ± 1.8 s (50- 
300keV) with a fluence of (3.26 ±0.05) x 10“®ergcm“^ (10- 
1000keV). The time-averaged 7 -ray spectrum is well fit by 
a Band function^, with break energy, Lpeak = 43.9 ± 3.9keV, 
low energy index, a = -0.67 ± 0.19, and high-energy index, 
P = -2.33 ±0.09 Using the source redshift of z = 1.798, the 
inferred isotropic equivalent 7 -ray energy in the 1-10^ keV 
rest frame energy band is = (3.15 ±0.12) x 10^^erg. 

2.1. X-ray: Swift/XRT 

The Swift X-ray Telescope (XRT, Burrows et al. 2005b) be¬ 
gan observing the held at 69 s after the BAT trigger, leading to 
the detection of an X-ray afterglow. The source was localized 
to RA = 18h 15m37.06s, Dec = H-69d 15' 35.4" (J2000), with 
an uncertainty radius of 1.4 arcseconds (90% containment)^. 

® From the Fermi GRB catalog for trigger 120326056 at http:// 
heasarc .gsfc.nasa. gov/db-perl/W3Browse/w3query .pi. 
http : //WWW . swift . ac . uk/xrt_positions/ 00518626/ 
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XRT continued observing the afterglow for 18.7 d in photon 
counting mode, with the last detection at 5.2 d. 

We extracted XRT PC-mode spectra using the on-line tool 
on the Swift website (Evans et al. 2007, 2009) We ana¬ 
lyzed the data after the end of the steep decay at 400 s us¬ 
ing the latest version of the HEASOET package (v6.14) and 
corresponding calibration files. We used Xspec (vl2.8.1) to 
fit all available PC-mode data, assuming a photoelectrically 
absorbed power law model (tbabs x ztbabs x pow) 
and a Galactic neutral hydrogen column density of A^h.mw = 
6.3 X 10^°cm“^ (for consistency with the value used on the 
Swift website), fixing the source redshift at z = 1.798. Our 
best-fit model has a photon index of E = 1.85 ± 0.04 and ex¬ 
cess absorption corresponding to a neutral hydrogen column 
(assuming solar metallicity) of AH.int = (4.1 ±0.7) x 10^'cm“^ 
intrinsic to the host galaxy (C-stat = 497.1 for 570 degrees 
of freedom). We divided the PC-mode data into 6 roughly 
equal time bins and extracted time-resolved PC-mode spectra 
to test for spectral evolution. We do not find clear evidence for 
significant spectral evolution over this period. In the follow¬ 
ing analysis, we take the 0.3 - lOkeV count rate light curve 


from the above website together with Fx = 1.85 (correspond¬ 
ing to /3x = 1 - Fx = -0.85) for the PC mode to compute the 
1 keV flux density. We combine the uncertainty in flux cal¬ 
ibration based on our spectral analysis (2.7%) in quadrature 
with the statistical uncertainty from the on-line light curve. 
For the WT-mode, we convert the count rate light curve to a 
flux-calibrated light curve using F = 3.57 and a count-to-flux 
conversion factor of 7 x 10“"ergcm“^ct“' as reported on the 
on-line spectral analysis. 

The WT-mode X-ray light curve declines rapidly as 
f-3.3±o.i Swift switched to collecting data in PC-mode at 
150 s. The PC-mode data between 150 and 300 s continue 
to decline rapidly as f“3.8±o.4^ followed by a plateau where the 
count rate evolves as fO 22±o.05 lightcurve 

is likely dominated by the high-latitude prompt emission 
(Willingale et al. 2010), and we therefore do not consider the 
rapid decline before 300 s in our afterglow modeling. About 
0.16 d after the trigger, the X-ray count rate begins rising and 
peaks at around 0.41 d. This re-brightening is unusual for X- 
ray afterglows and we discuss this feature further in section 
3.1. The XRT count rate light curve after 1.4 d can be fit by a 
single power law with a decline rate of ax = -2.29 ±0.16. 


Table 1 

SwiftnjVtJI Observations of GRB 120326A 


(days) 

Filter 

Frequency 

(Hz) 

Flux density^ 
(AiJy) 

Uncertainty 

(/iJy) 

Detection? 

(1 = Yes) 

0.00165 

White 

8.64e+14 

8.89 

2.71 

1 

0.00354 

u 

8.56e+14 

20.7 

11.5 

0 

0.0523 

b 

6.92e+14 

39 

6.23 

1 

0.0546 

White 

8.64e+14 

15.6 

2.55 

1 

0.0558 

u\^l 

1.16e+15 

10.9 

3.29 

1 

0.0582 

u 

8.56e+14 

18.6 

4.24 

1 

0.0594 

V 

5.55e+14 

42.3 

14.2 

0 

0.115 

uvwl 

1.16e+15 

9.16 

1.54 

1 

0.126 

u 

8.56e+14 

30.1 

2.34 

1 

0.134 

b 

6.92e+14 

68.2 

8.04 

1 

0.191 

uvwl 

1.16e+15 

13.8 

2.12 

1 

0.202 

u 

8.56e+14 

50.6 

6.19 

1 

0.258 

V 

5.55e+14 

134 

11.8 

1 

0.285 

uvw2 

1.48e+15 

1.95 

0.545 

1 

0.295 

uvm2 

1.34e+15 

0.607 

0.549 

0 

0.316 

V 

5.55e+14 

150 

10.6 

1 

0.383 

uvwl 

1.16e+15 

19.1 

1.94 

1 

0.393 

u 

8.56e+14 

59.2 

3.07 

1 

0.399 

b 

6.92e+14 

124 

13.6 

1 

0.46 

uvwl 

1.16e+15 

12.3 

1.66 

1 

0.468 

u 

8.56e+14 

43.1 

4.43 

1 

0.527 

v 

5.55e+14 

97.6 

9.45 

1 

0.715 

v 

5.55e+14 

83.9 

10.4 

1 

0.715 

uvwl 

1.16e+15 

11.7 

1.35 

1 

0.737 

u 

8.56e+14 

28.9 

4.06 

1 

0.782 

uvw2 

1.48e+15 

0.968 

0.68 

0 

0.783 

h 

6.92e+14 

50.1 

4.38 

1 

0.792 

White 

8.64e+14 

26.3 

3.74 

1 

1.65 

u 

8.56e+14 

11.6 

2.31 

1 

1.65 

uvw2 

1.48e+15 

0.339 

0.668 

0 

1.66 

uvm2 

1.34e+15 

2.42 

1.71 

0 

1.66 

uvwl 

1.16e+15 

6.29 

2.29 

0 

2.49 

uvwl 

1.16e+15 

3.53 

1.67 

0 

2.49 

u 

8.56e+14 

8.3 

1.83 

1 

2.69 

b 

6.92e+14 

16 

4.01 

1 

2.7 

White 

8.64e+14 

6.59 

1.29 

1 

2.7 

V 

5.55e+14 

31.9 

8.93 

1 

3.72 

u 

8.56e+14 

5.93 

1.62 

1 

3.73 

b 

6.92e+14 

7.59 

3.43 

0 

3.73 

uvwl 

1.16e+15 

0.779 

0.626 

0 

3.73 

White 

8.64e+14 

4.26 

0.985 

1 

3.74 

V 

5.55e+14 

5.94 

7.16 

0 

4.83 

u 

8.56e+14 

2.92 

1.42 

0 


http : //www . swift . ac . uk/xrt_spectra/0 0518 62 6/ 
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Table 1 — Continued 


t-tQ 

(days) 

Filter 

Frequency 

(Hz) 

Flux densityt 

(/rJy) 

Uncertainty 

(mU) 

Detection? 

(1 = Yes) 

4.83 

b 

6.92e-Hl4 

8.72 

3.14 

0 

4.84 

White 

8.64e-i-14 

2.06 

0.768 

0 

4.84 

V 

5.55e-Hl4 

5.12 

6.37 

0 

5.73 

u 

8.56e-Hl4 

0.689 

1.39 

0 

5.73 

b 

6.92e-Hl4 

6.41 

3.13 

0 

5.74 

White 

8.64e-l-14 

0.754 

0.719 

0 

5.74 

V 

5.55e-Hl4 

0.413 

6.26 

0 

17.6 

u 

8.56e-Hl4 

0.28 

1.39 

0 

19.5 

u 

8.56e-Hl4 

2.34 

1.28 

0 


Note. — ^In cases of non-detections, we report the formal flux density 
measurement from aperture photometry. 

2.2. UV/Optical: Swift/UVOT 

The Swift UV/Optical Telescope (UVOT; Roming et al. 
2005) observed GRB 120326A using 6 filters spanning the 
central wavelength range Ac = 1928 A {W2) to Ac = 5468 A 
(v) beginning 67 s after the burst. The afterglow was de¬ 
tected in the initial exposures at RA = 18h 15m 37.13s, Dec = 
H-69d 15' 35.36" (J2000) (90% confidence, Kuin et al. 2012), 
and exhibited a clear re-brightening concomitant with the 
peak in the X-ray light curve. We analyzed the UVOT data 
using the latest version of HEASOFT (v. 6.14) and corre¬ 
sponding calibration files. We performed photometry with a 
5" aperture and used a 90" annulus with foreground sources 
masked to estimate the background. We list our derived fluxes 
in Table 1. 

2.3. Optical/NIR: Palomar Observations and GCN 
Circulars 

We observed GRB 120326A with the Wide-Field Infrared 
Camera (Wilson et al. 2003) on the Palomar 200-inch tele¬ 


scope beginning on 2012 March 30.48 UT. We acquired a 
series of nine 75 s 7-band exposures, followed by a series of 
nine 75 s fC-band exposures. We reduced the data following 
standard IR imaging techniques using a modified version of 
the WIRCSOFT pipeline and performed aperture photometry 
of the afterglow in the combined, stacked exposures relative 
to 2MASS standards in the field. 

We also carried out a series of imaging observations of 
GRB 120326A with the robotic Palomar 60-inch telescope 
(Cenko et al. 2006) on the nights of 2012 March 31, April 02, 
April 03 (r' and /'), April 04, April 08, and April 16 (r' only). 
We reduced the images using the P60 automated pipeline and 
performed aperture photometry relative to secondary standard 
stars calibrated via separate observations of Landolt (2009) 
standards. 

Finally, we collected other optical and NIR observations of 
GRB 120326A reported through the Gamma-ray Burst Coor¬ 
dinates Network (GCN) Circulars and converted all photom¬ 
etry to flux densities. We also include the optical photome¬ 
try of the afterglow published by Melandri et al. (2014) and 
Urata et al. (2014) in our analysis. We list our complete com¬ 
pilation of optical observations of this burst, together with our 
P60 and P200 observations, in Table 2. 


Table 2 

Optical Observations of GRB 120326A 


t-?0 

(days) 

Observatory 

Telescope / 
Instmment 

Filter 

Frequency 

(Hz) 

Flux density^ 
(AiJy) 

Uncertainty 

(uJy) 

Detection? 
(1 = Yes) 

Reference 

0.000291 

TNO 

ROTSE 

CR 

4.56e-Hl4 

4.04e-H03 

1.35e-H03 

0 

Rujopakam & Flewelling (2012) 

0.001 

TNO 

ROTSE 

CR 

4.56e-Hl4 

926 

308 

0 

Rujopakam & Flewelling (2012) 

0.00189 

Calem 

TAROT 

R 

4.56e-Hl4 

161 

51.2 

1 

Klotz et al. (2012) 

0.00389 

Liverpool 

RING02 

R 

4.56e-Hl4 

120 

21.6 

1 

Melandri et al. (2014) 

0.00545 

Cmi Vrh 

Cicocki 

R 

4.56e-Hl4 

95.1 

13.1 

1 

Melandri et al. (2014) 

0.00623 

Crni Vrh 

Cicocki 

R 

4.56e-Hl4 

87.6 

16.7 

1 

Melandri et al. (2014) 

0.007 

Cmi Vrh 

Cicocki 

R 

4.56e+14 

86.8 

12.9 

1 

Melandri et al. (2014) 

0.00813 

Crni Vrh 

Cicocki 

R 

4.56e-Hl4 

74.9 

14.3 

1 

Melandri et al. (2014) 

0.00824 

Liverpool 

RATCam 

R 

4.56e-Hl4 

97.8 

17.6 

1 

Melandri et al. (2014) 

0.00968 

Cmi Vrh 

Cicocki 

R 

4.56e+14 

87.6 

14.8 

1 

Melandri et al. (2014) 

0.00991 

Liverpool 

RATCam 

r’ 

4.81e-l-14 

109 

17.2 

1 

Melandri et al. (2014) 

0.0109 

Liverpool 

RATCam 

r’ 

4.81e-l-14 

91.2 

12.6 

1 

Melandri et al. (2014) 

0.0116 

Cmi Vrh 

Cicocki 

R 

4.56e+14 

82.9 

8.84 

1 

Melandri et al. (2014) 

0.0122 

Liverpool 

SkyCam 

R 

4.56e-Hl4 

90 

37.7 

1 

Melandri et al. (2014) 

0.0123 

Liverpool 

RATCam 

r' 

4.81e-i-14 

95.5 

3.58 

1 

Melandri et al. (2014) 

0.0136 

Cmi Vrh 

Cicocki 

R 

4.56e+14 

70.2 

10.4 

1 

Melandri et al. (2014) 

0.0142 

Liverpool 

RATCam 

r 

3.93e-Hl4 

86.3 

5.75 

1 

Melandri et al. (2014) 

0.0151 

Cmi Vrh 

Cicocki 

R 

4.56e-Hl4 

54.8 

13.5 

1 

Melandri et al. (2014) 

0.016 

Liverpool 

RATCam 

z' 

3.28e-i-14 

111 

14.1 

1 

Melandri et al. (2014) 

0.017 

Cmi Vrh 

Cicocki 

R 

4.56e-Hl4 

49.5 

8.92 

1 

Melandri et al. (2014) 

0.0178 

Liverpool 

RATCam 

r’ 

4.81e-i-14 

72.4 

3.41 

1 

Melandri et al. (2014) 

0.0193 

Liverpool 

RATCam 

r’ 

4.81e-i-14 

67.3 

3.82 

1 

Melandri et al. (2014) 

0.0193 

Cmi Vrh 

Cicocki 

R 

4.56e-Hl4 

54.3 

12.2 

1 

Melandri et al. (2014) 

0.0212 

Liverpool 

RATCam 

r 

3.93e-Hl4 

68.5 

3.89 

1 

Melandri et al. (2014) 

0.0228 

Calern 

TAROT 

R 

4.56e+14 

70.2 

22.3 

1 

Klotz et al. (2012) 

0.0254 

Liverpool 

RATCam 

z' 

3.28e-i-14 

103 

9.92 

1 

Melandri et al. (2014) 

0.0259 

DAO 

Skynet 

g’ 

6.29e-i-14 

17.9 

5.69 

1 

Lacluyze et al. (2012) 

0.0277 

DAO 

Skynet 

r’ 

4.56e-i-14 

41 

13 

1 

Lacluyze et al. (2012) 
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Table 2 — Continued 


t-to 

(days) 

Observatory 

Telescope / 
Instnament 

Eilter 

Frequency 

(Hz) 

Flux density'!’ 
(AiJy) 

Uncertainty 

(/iJy) 

Detection? 

(1 = Yes) 

Reference 

0.028 

Liverpool 

RATCam 

r’ 

4.8Ie+14 

63.1 

3.58 

1 

Melandri et al. (2014) 

0.0296 

Liverpool 

RATCam 

r’ 

4.81e+I4 

63.1 

2.97 

1 

Melandri et al. (2014) 

0.0311 

Liverpool 

RATCam 

r’ 

4.81e+I4 

64.3 

3.65 

1 

Melandri et al. (2014) 

0.0333 

Liverpool 

SkyCam 

R 

4.56e+14 

51.3 

43.8 

1 

Melandri et al. (2014) 

0.0344 

Liverpool 

RATCam 

r 

3.93e+14 

63.1 

4.2 

1 

Melandri et al. (2014) 

0.0358 

DAO 

Skynet 

r 

3.93e+14 

73.4 

23.4 

1 

Lacluyze et al. (2012) 

0.0396 

Liverpool 

RATCam 

z’ 

3.28e+I4 

87.1 

8.4 

1 

Melandri et al. (2014) 

0.0709 

T17 

CDK17 

CR 

4.56e+14 

77 

15.6 

1 

Hentunen et al. (2012) 

0.32 

LOAO 

Im 

R 

4.56e+14 

272 

15.4 

1 

Jang et al. (2012) 

0.324 

LOAO 

Im 

R 

4.56e+14 

298 

11.2 

1 

Jang et al. (2012) 

0.327 

LOAO 

Im 

R 

4.56e+14 

277 

10.4 

1 

Jang et al. (2012) 

0.331 

LOAO 

Im 

R 

4.56e+14 

224 

6.28 

1 

Jang et al. (2012) 

0.335 

LOAO 

Im 

R 

4.56e+14 

212 

5.94 

1 

Jang et al. (2012) 

0.338 

LOAO 

Im 

R 

4.56e+14 

180 

5.03 

1 

Jang et al. (2012) 

0.367 

McDonald 

CQUEAN 

r’ 

4.81e+I4 

173 

2.2 

1 

Urata et al. (2014)* 

0.371 

McDonald 

CQUEAN 

i’ 

3.93e+14 

181 

1 

1 

Urata et al. (2014)^ 

0.375 

McDonald 

CQUEAN 

z’ 

3.28e+14 

233 

1.6 

1 

Urata et al. (2014)* 

0.379 

McDonald 

CQUEAN 

y 

2.94e+14 

264 

4 

1 

Urata et al. (2014)t 

0.383 

McDonald 

CQUEAN 

r’ 

4.8Ie+14 

164 

0.8 

1 

Urata et al. (2014)t 

0.387 

McDonald 

CQUEAN 

V 

3.93e+14 

185 

0.5 

1 

Urata et al. (2014)* 

0.391 

McDonald 

CQUEAN 

z’ 

3.28e+14 

231 

2.1 

1 

Urata et al. (2014)* 

0.398 

McDonald 

CQUEAN 

r’ 

4.8Ie+14 

151 

2.7 

1 

Urata et al. (2014)t 

0.402 

McDonald 

CQUEAN 

i’ 

3.93e+14 

190 

2.3 

1 

Urata et al. (2014)* 

0.405 

McDonald 

CQUEAN 

z’ 

3.28e+14 

225 

7.3 

1 

Urata et al. (2014)* 

0.418 

McDonald 

CQUEAN 

g' 

6.29e+I4 

115 

3.39 

1 

Urata et al. (2014)t 

0.425 

McDonald 

CQUEAN 

r’ 

4.8Ie+14 

148 

1 

1 

Urata et al. (2014)t 

0.429 

McDonald 

CQUEAN 

V 

3.93e+14 

166 

3.9 

1 

Urata et al. (2014)t 

0.737 

GMG 

2.4m 

R 

4.56e+14 

101 

9.79 

1 

Zhao et al. (2012) 

0.742 

GMG 

2.4m 

V 

5.55e+14 

148 

14.2 

1 

Zhao et al. (2012) 

0.747 

Lulin 

LOT 

V 

3.93e+14 

99.1 

2.1 

1 

Urata et al. (2014)t 

0.747 

Lulin 

LOT 

g’ 

6.29e+I4 

55.3 

0.212 

1 

Urata et al. (2014)^ 

0.75 

GMG 

2.4m 

R 

4.56e+14 

92.6 

8.93 

1 

Zhao et al. (2012) 

0.754 

GMG 

2.4m 

V 

5.55e+14 

135 

13 

1 

Zhao et al. (2012) 

0.758 

Lulin 

LOT 

r’ 

4.81e+I4 

85.2 

1.11 

1 

Urata et al. (2014)t 

0.758 

Lulin 

LOT 

z’ 

3.28e+14 

132 

8.8 

1 

Urata et al. (2014)t 

0.794 

Lulin 

LOT 

V 

3.93e+14 

98.2 

2.8 

1 

Urata et al. (2014)* 

0.806 

Lulin 

LOT 

z’ 

3.28e+14 

109 

8.5 

1 

Urata et al. (2014)^ 

0.81 

Xinglong 

TNT 

R 

4.56e+14 

115 

36.7 

1 

Xin et al. (2012a) 

0.842 

Maisoncelles 

30cm 

R 

4.56e+14 

111 

10.7 

1 

Sudilovsky et al. (2012) 

0.92 

BBO 

320mm 

R 

4.56e+14 

140 

9.33 

1 

Melandri et al. (2014) 

1.36 

McDonald 

CQUEAN 

r’ 

4.81e+I4 

49.5 

0.36 

1 

Urata et al. (2014)* 

1.36 

McDonald 

CQUEAN 

i’ 

3.93e+14 

62.9 

0.3 

1 

Urata et al. (2014)^ 

1.37 

McDonald 

CQUEAN 

z’ 

3.28e+14 

75.6 

0.6 

1 

Urata et al. (2014)* 

1.37 

McDonald 

CQUEAN 

y 

2.94e+14 

83.4 

2.5 

1 

Urata et al. (2014)t 

1.38 

McDonald 

CQUEAN 

g’ 

6.29e+I4 

31 

0.72 

1 

Urata et al. (2014)t 

1.39 

LOAO 

Im 

R 

4.56e+14 

89.3 

15.8 

1 

Urata et al. (2014)* 

1.43 

FLWO 

PAIRITEL 

K 

1.37e+14 

248 

64.1 

1 

Melandri et al. (2014) 

1.43 

FLWO 

PAIRITEL 

H 

I.84e+14 

308 

137 

1 

Melandri et al. (2014) 

1.43 

FLWO 

PAIRITEL 

J 

2.38e+14 

121 

53.7 

1 

Melandri et al. (2014) 

1.66 

Ishigakijima 

MITSuME 

R 

4.56e+14 

52.8 

4.56 

1 

Melandri et al. (2014) 

1.66 

Ishigakijima 

MITSuME 

g’ 

6.29e+I4 

21.3 

2.93 

1 

Melandri et al. (2014) 

1.71 

Ishigakijima 

MITSuME 

R 

4.56e+14 

47.7 

4.12 

1 

Melandri et al. (2014) 

1.71 

Ishigakijima 

MITSuME 

g’ 

6.29e+I4 

20.5 

3.04 

1 

Melandri et al. (2014) 

1.82 

Xinalona 

TNT 

R 

4.56e+14 

57.3 

18.2 

1 

Xin et al. (2012a) 

1.9 

Hanle 

HCT 

R 

4.56e+14 

41.2 

2.74 

1 

Sahu et al. (2012) 

2.01 

BBO 

320mm 

R 

4.56e+14 

56.8 

8.42 

1 

Melandri et al. (2014) 

2.15 

Ishigakijima 

MITSuME 

I 

3.93e+I4 

36.6 

7.4 

1 

Melandri et al. (2014) 

2.34 

McDonald 

CQUEAN 

r’ 

4.81e+I4 

23.8 

3.68 

1 

Urata et al. (2014)^ 

2.35 

McDonald 

CQUEAN 

i’ 

3.93e+14 

33.3 

5.34 

1 

Urata et al. (2014)t 

2.35 

McDonald 

CQUEAN 

z’ 

3.28e+14 

42.7 

15.4 

1 

Urata et al. (2014)* 

2.38 

McDonald 

CQUEAN 

y 

2.94e+14 

59.6 

5.4 

1 

Urata et al. (2014)* 

2.4 

LOAO 

Im 

R 

4.56e+14 

41.5 

7.08 

1 

Urata et al. (2014)1 

2.4 

McDonald 

CQUEAN 

g' 

6.29e+I4 

18.6 

1.81 

1 

Urata et al. (2014)1 

2.42 

FLWO 

PAIRITEL 

K 

1.37e+14 

168 

43.6 

1 

Melandri et al. (2014) 

2.42 

FLWO 

PAIRITEL 

H 

I.84e+14 

123 

71.8 

1 

Melandri et al. (2014) 

2.42 

FLWO 

PAIRITEL 

J 

2.38e+14 

43.8 

25.6 

1 

Melandri et al. (2014) 

2.43 

McDonald 

CQUEAN 

r' 

4.81e+I4 

26.1 

1.3 

1 

Urata et al. (2014)1 

2.43 

McDonald 

CQUEAN 

i’ 

3.93e+14 

30.5 

2.71 

1 

Urata et al. (2014)1 

2.44 

McDonald 

CQUEAN 

z’ 

3.28e+14 

42.8 

5.47 

1 

Urata et al. (2014)1 

2.6 

Ishigakijima 

MITSuME 

R 

4.56e+14 

30.6 

5.19 

1 

Melandri et al. (2014) 

2.6 

Ishigakijima 

MITSuME 

g' 

6.29e+I4 

14.3 

4.91 

1 

Melandri et al. (2014) 
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Table 2 — Continued 


t-to 

(days) 

Observatory 

Telescope / 
Instrument 

Eilter 

Frequency 

(Hz) 

Flux density'!’ 
(AiJy) 

Uncertainty 

(/7Jy) 

Detection? 

(1 = Yes) 

Reference 

3 

BBO 

320mm 

R 

4.56e+14 

28 

9.24 

1 

Melandri et al. (2014) 

3.32 

McDonald 

CQUEAN 

r 

3.93e-Hl4 

21.2 

0.864 

1 

Urata et al. (2014)* 

3.33 

McDonald 

CQUEAN 

z 

3.28e-Hl4 

28.1 

2.03 

1 

Urata et al. (2014)* 

3.33 

McDonald 

CQUEAN 

Y’ 

2.94e-Hl4 

44.5 

5.3 

1 

Urata et al. (2014)t 

3.34 

McDonald 

CQUEAN 

r’ 

4.81e-Hl4 

16.8 

0.868 

1 

Urata et al. (2014)* 

3.38 

McDonald 

CQUEAN 

g’ 

6.29e-tl4 

11.6 

1.33 

1 

Urata et al. (2014)* 

3.43 

McDonald 

CQUEAN 

i’ 

3.93e-Hl4 

21.2 

5.02 

1 

Urata et al. (2014)t 

3.49 

LOAD 

Im 

R 

4.56e-Hl4 

23 

3.74 

1 

Urata et al. (2014)* 

3.71 

Ishigakijima 

MITSuME 

R 

4.56e-Hl4 

16.7 

2.47 

1 

Melandri et al. (2014) 

3.73 

Ishigakijima 

MITSuME 

g’ 

6.29e-tl4 

6.31 

1.21 

1 

Melandri et al. (2014) 

3.73 

Ishigakijima 

MITSuME 

I 

3.93e-tl4 

19.4 

4.35 

1 

Melandri et al. (2014) 

3.77 

Ishigakijima 

MITSuME 

R 

4.56e-Hl4 

16.8 

2.32 

1 

Melandri et al. (2014) 

4.4 

McDonald 

CQUEAN 

z’ 

3.28e-tl4 

13.7 

2.3 

1 

Urata et al. (2014)t 

4.41 

McDonald 

CQUEAN 

V 

3.93e-Hl4 

10.8 

1.01 

1 

Urata et al. (2014)t 

4.42 

McDonald 

CQUEAN 

r’ 

4.81e-Hl4 

7.88 

0.742 

1 

Urata et al. (2014)t 

4.43 

McDonald 

CQUEAN 

g' 

6.29e-tl4 

5.2 

0.32 

1 

Urata et al. (2014)* 

4.43 

Palomar 

P200 

j 

2.38e-Hl4 

31.2 

3.65 

1 

This work 

4.45 

Palomar 

P200 

K 

1.37e-Hl4 

33.1 

5.97 

1 

This work 

5.25 

P60 


/’ 

3.93e-Hl4 

6.79 

0.519 

1 

This work 

5.26 

Palomar 

P60 

r’ 

4.81e-Hl4 

6.14 

0.409 

1 

This work 

5.38 

McDonald 

CQUEAN 

z’ 

3.28e-tl4 

6.35 

3.62 

1 

Urata et al. (2014)t 

5.39 

McDonald 

CQUEAN 

i’ 

3.93e-Hl4 

6.39 

1.06 

1 

Urata et al. (2014)t 

5.41 

McDonald 

CQUEAN 

r’ 

4.81e-Hl4 

4.65 

0.735 

1 

Urata et al. (2014)* 

5.42 

McDonald 

CQUEAN 

g’ 

6.29e-tl4 

3.1 

0.42 

1 

Urata et al. (2014)t 

5.92 

CRO 

51cm 

R 

4.56e-Hl4 

7.56 

3.17 

1 

Melandri et al. (2014) 

5.92 

CRO 

51cm 

CR 

4.56e-Hl4 

5.84 

0.563 

1 

Voinova et al. (2012) 

6.4 

McDonald 

CQUEAN 

r' 

4.81e-tl4 

2.96 

0.0931 

1 

Urata et al. (2014)* 

6.41 

McDonald 

CQUEAN 

V 

3.93e-Hl4 

4.13 

0.526 

1 

Urata et al. (2014)* 

6.42 

McDonald 

CQUEAN 

z’ 

3.28e-Hl4 

5.11 

2.44 

1 

Urata et al. (2014)t 

6.83 

LOAD 

Im 

R 

4.56e-Hl4 

5.75 

0.88 

1 

Urata et al. (2014)t 

7.24 

P60 


i’ 

3.93e-Hl4 

6.43 

0.885 

1 

This work 

7.25 

Palomar 

P60 

r’ 

4.81e-Hl4 

2.81 

0.759 

1 

This work 

7.41 

McDonald 

CQUEAN 

r’ 

4.81e-tl4 

1.97 

0.3 

1 

Urata et al. (2014)^ 

7.51 

McDonald 

CQUEAN 

/’ 

3.93e-Hl4 

2.44 

0.29 

1 

Urata et al. (2014)t 

8.29 

P60 


V 

3.93e-Hl4 

3.37 

0.5 

1 

This work 

8.3 

Palomar 

P60 

r’ 

4.81e-Hl4 

2.47 

0.392 

1 

This work 

9.29 

Palomar 

P60 

r’ 

4.81e-tl4 

2.44 

0.236 

1 

This work 

13.3 

Palomar 

P60 

r’ 

4.81e-tl4 

1.98 

0.4 

1 

This work 

21.4 

Palomar 

P60 

r’ 

4.81e-tl4 

1.64 

0.11 

1 

This work 


Note. — tin cases of non-detections, we report the 3(7 upper limit on the 
flux density, t We note that the photometry reported in Table 1 of Urata et al. 
(2014) yields an unphysical spectral index of /3 Ri -4.5. We assume that 
these numbers have been scaled by the same factors as reported in Figure 4 
of their paper, and divide the gr'Rizy light curves by a factor of 0,5, 1, 1, 2,4, 
and 6, respectively. We further assume that the data have not been corrected 
for galactic extinction. 


2.4. Sub-millimeter: CARMA and SMA 

We observed GRB 120326A with the Combined Array for 
Research in Millimeter Astronomy (CARMA) beginning on 
2012 March 30.56 UT (4.52 d after the burst) in continuum 
wideband mode with w 8 GHz bandwidth (16 windows, 487.5 
MHz each) at a mean frequency of 93 GHz. Following an 
initial detection (Perleyetal. 2012), we obtained two addi¬ 
tional epochs. Weather conditions were excellent for each 
epoch, with the first observation being taken in C configu¬ 
ration (maximum baseline ~370 m), and the following obser¬ 
vations in D configuration (maximum baseline ~145 m). For 
all observations, we utilized 3C371 for bandpass, amplitude 
and phase gain calibration, and MWC349 for flux calibration. 
We reduced the data using standard procedures in MIRJAD 
(Sault et al. 1995). In the last epoch on April 10, we also ob¬ 
served 3C273 and 3C279 and utilized these observations as 


independent checks on the bandpass and flux calibration. We 
summarize our observations in Table 3. 

The Submillimeter Array (SMA; Ho et al. 2004) observed 
GRB 120326A at a mean frequency of 222 GHz (1.3 mm; 
lower sideband centered at 216 GHz, upper sideband at 228 
GHz; PI: Urata). Six epochs of observations were obtained 
between 2012 March 26.43 UT (0.37 d after the burst) and 
2012 April 11.64 UT (16.6d after the burst). These obser¬ 
vations have been reported in Urata et al. (2014). We carried 
out an independent reduction of the data using standard MIR 
IDL procedures for the SMA, followed by flagging, imaging 
and analysis in MIRIAD and the Astronomical Image Pro¬ 
cessing System (AIPS; Greisen 2003). We utilized 3C279 
for bandpass calibration in all but the last epoch, where we 
utilized J1924-292. We utilized J1800H-784 and J1829H-487 
for gain calibration®. We utilized MWC349a for flux cali¬ 
bration, determining a flux density of 1.42 Jy for J1800-1-784, 
consistent with flux values in the SMA catalog at similar 
times. We utilized this flux throughout all observations (not 
every epoch contained useful flux calibrator scans) and scaled 
the gains appropriately. We note an uncertainty in the abso- 

^ The gain calibrators were separated by more than 20 degrees from the 
source; thus some decoherence, resulting in a systematic reduction of the 
observed flux is possible. 
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lute flux density scale of r; 15%. We detect a source at RA 
= 18hl5m37.15s (±0.02)s, Dec = +69d 15'35.41" (±0.14; 
J2000). Two epochs obtained on March 31.40 and April 6.51 
showed poor noise characteristics, and we do not include them 
in our analysis. We measure poorer noise statistics in the ob¬ 
servations than reported by Urata et al. (2014), and also do 
not detect the source in the epoch at 3.53 d, contrary to the 
previous analysis of this dataset. We report the results of our 
analysis in Table 3. 

2.5. Radio: VIA 

We observed the afterglow at C (4-7 GHz), K (18-25 GHz), 
and Ka (30-38 GHz) bands using the Karl G. Jansky Very 


Large Array (VLA) starting 5.45 d after the burst. We de¬ 
tected and tracked the flux density of the afterglow over eight 
epochs until 122 d after the burst, until the source had either 
faded below or was barely detectable at the 3(T-level at all fre¬ 
quencies. Depending on the start time of the observations, we 
used either 3C286 or 3C48 as the flux and bandpass calibra¬ 
tor; we used J1806H-6949 or J1842H-6809 as gain calibrator de¬ 
pending on the array configuration and observing frequency. 
We carried out data reduction using the Common Astronomy 
Software Applications (CASA). We list our VLA photometry 
in Table 3. 


Table 3 

Millimeter and Radio Observations of GRB 120326A 


Date 

(UT) 

t-to 

(days) 

Observatoiy* 

Frequency 

(GHz) 

Integration time 
(min) 

Flux density 
(mJy) 

Uncertainty^ 

(mJy) 

Detection? 

(1 = Yes) 

2012 Mar 30.56 

4.55 

CAItMA 

92.5 

60.3 

3.38 

0.87 

1 

2012 Apr 5.53 

10.51 

CARMA 

92.5 

73.0 

1.46 

0.42 

1 

2012 Apr 10.40 

15.51 

CARMA 

92.5 

291.9 

0.476 

0.17 

1 

2012 Mar 26.43 

0.504 

SMA* 

222 

386 

3.3 

0.90 

1 

2012 Mar 27.53 

1.55 

SMA* 

222 

227 

2.4 

0.80 

1 

2012 Mar 29.50 

3.53 

SMAt 

222 

252 

< 1.8 

0.60 

0 

2012 Apr 11.64 

16.6 

SMA* 

222 

201 

< 2.4 

0.80 

0 

2012 Mar 31.55 

5.66 

VLA/C 

5.0 

8.54 

0.674 

0.017 

1 

2012 Mar 31.55 

5.66 

VLA/C 

7.1 

4.87 

0.369 

0.018 

1 

2012 Mar 31.53 

5.64 

VLA/C 

19.2 

9.35 

1.25 

0.031 

1 

2012 Mar 31.53 

5.64 

VLA/C 

24.5 

7.70 

1.49 

0.038 

1 

2012 Apr 04.48 

9.42 

VLA/C 

5.0 

5.37 

0.411 

0.024 

1 

2012 Apr 04.48 

9.42 

VLA/C 

7.1 

6.00 

0.561 

0.017 

1 

2012 Apr 04.46 

9.41 

VLA/C 

19.2 

7.16 

0.825 

0.028 

1 

2012 Apr 04.46 

9.40 

VLA/C 

24.5 

8.46 

0.936 

0.039 

1 

2012 Apr 04.45 

9.40 

VLA/C 

33.5 

7.16 

0.94 

0.041 

1 

2012 Apr 10.55 

15.49 

VLA/C 

5.0 

7.59 

0.313 

0.019 

1 

2012 Apr 10.55 

15.49 

VLA/C 

7.1 

6.72 

0.352 

0.016 

1 

2012 Apr 10.53 

15.48 

VLA/C 

19.2 

7.87 

0.667 

0.035 

1 

2012 Apr 10.53 

15.48 

VLA/C 

24.5 

7.88 

0.686 

0.045 

1 

2012 Apr 10.52 

15.46 

VLA/C 

33.5 

7.87 

0.706 

0.053 

1 

2012 Apr 26.42 

31.37 

VLA/C 

5.0 

5.39 

0.409 

0.022 

1 

2012 Apr 26.42 

31.37 

VLA/C 

7.1 

4.55 

0.562 

0.020 

1 

2012 Apr 26.41 

31.35 

VLA/C 

19.2 

13.34 

0.458 

0.025 

1 

2012 Apr 26.41 

31.35 

VLA/C 

24.5 

9.43 

0.401 

0.032 

1 

2012 Apr 26.39 

31.34 

VLA/C 

33.5 

7.86 

0.277 

0.034 

1 

2012 May 27.49 

62.44 

VLA/B 

5.0 

6.11 

0.217 

0.020 

1 

2012 May 27.49 

62.44 

VLA/B 

7.1 

5.10 

0.283 

0.017 

1 

2012 May 27.46 

62.40 

VLA/B 

19.2 

20.25 

0.183 

0.017 

1 

2012 May 27.46 

62.40 

VLA/B 

24.5 

21.98 

0.142 

0.019 

1 

2012 Jul 26.17 

122.1 

VLA/B 

5.0 

9.81 

0.057 

0.017 

1 

2012 Jul 26.17 

122.1 

VLA/B 

7.1 

9.77 

0.072 

0.014 

1 

2012 Jul 26.12 

122.1 

VLA/B 

19.2 

31.85 

0.052 

0.016 

1 

2012 Jul 26.12 

122.1 

VLA/B 

24.5 

26.57 

<0.07 

0.023 

0 


Note. — * The letter following ‘VLA’ indicates the array configuration, 
t 1(7 statistical uncertainties from AIPS task JMFIT (CARMA and SMA 
observations) or CASA task IMFIT (VLA observations), tpi: Y. Urata. 

3. BASIC CONSIDERATIONS 

We consider the X-ray to radio afterglow of GRB 120326A 
in the context of the standard synchrotron afterglow model, 
where the afterglow emission is produced by synchrotron ra¬ 
diation from a non-thermal distribution of electrons. The elec¬ 
tron energy distribution is assumed to be a power law with 
electron number density, n{'-f) oc 7“^, for 7 > 7min; here 7n,i„ 
is the minimum Lorentz factor and the electron energy in¬ 
dex, p is expected to lie between 2 and 3. This electron 
distribution results in a synchrotron spectrum that can be de¬ 
scribed by a series of power law segments connected at ‘break 


frequencies’: the self-absorption frequency, below which 
synchrotron self-absorption suppresses the flux, the character¬ 
istic synchrotron frequency, i/m, which corresponds to emis¬ 
sion from electrons with 7 = 7tnin, and the cooling frequency, 
t'c, above which synchrotron cooling is important; the flux 
density at r'u, sets the overall flux normalization. The syn¬ 
chrotron model is described in detail in Sari et al. (1998), 
Chevalier & Li (2000) and Granot & Sari (2002). 


3.1. X-ray/UV/Optical Re-brightening at 0.4 d 

A prominent feature of the afterglow light curve of 
GRB 120326A is a re-brightening at about 0.4 d. We quan¬ 
tify the shape of the light curve at the re-brightening by fitting 
the X-ray data with a smoothly-joined broken power law of 
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Figure 1. Broken power law fits to the X-ray (grey, solid), t/-band (blue, 
dashed), and z-band (red, dotted) light curves for GRB 120326A near the 
re-brightening around 0.35 d. X-ray points before 0.15 d are not included in 
the fit. En'orbai's at z-band are typically smaller than the size of the plotted 
symbols. The lines correspond to the independent fits at the three frequencies. 
The best fit parameters are listed in Table 4. 


the form 


F,,=Fu 




-l/v 


( 1 ) 


where fb is the break time, F\, is the flux at the break time, 
ai and a 2 are the temporal indices before and after the break, 
respectively, and y is the sharpness of the break The X- 
ray data before 0.15 d exhibit a plateau and we therefore re¬ 
strict the fit to span 0.15 to 1.25 d. We use the Python func¬ 
tion curve_fit to estimate these model parameters and 
the associated covariance matrix. Our best-fit parameters are 
fb = (0.41 ±0.02)d, ai = 0.85 ±0.19, a 2 = -1.22 ±0.18, and 
y = 5 ± 4 (Figure 1 and Table 4). 

The optical data are more sparsely sampled than the X- 
ray observations. We fit the U -band data between 0.05 d and 
1.25 d after fixing y = 5 as suggested by the fit to the better- 
sampled X-ray light curve. Our derived value of the peak 
time, fb = 0.31 ±0.04 is marginally earlier than, but close to 
the peak time of the X-ray light curve. The rise and decay 
rate in t/-band are also statistically consistent with those de¬ 
rived from the X-ray light curve. Similarly, we fit the z-band 
light curve between 0.04d and 1.4d, fixing y = 5. Since this 
time range includes a single point prior to the peak, the peak 
time and rise rate are degenerate in the fit. We therefore fix 
fb = 0.31 as derived from the t/-band fit. The best-fit parame¬ 
ters are listed in Table 4, and are consistent with those derived 
for the UV and X-ray light curves. 

Finally, we fit the X-ray, UV, and interpolated optical/NIR 
data jointly, where the three light curves are constrained to 
the same rise and decay rate, time of peak, and sharpness of 
the break, with independent normalizations in the three bands. 


Using a Markov Chain Monte Carlo simulation using EMCEE 
(Foreman-Mackey et al. 2013), we find fb = 0.40 ±0.01, ai = 
0.52 ±0.06, and az = -1-10 ±0.10. 

In summary, the X-ray/UV/optical light curves exhibit a 
prominent peak, nearly simultaneously in multiple bands (X- 
rays through the optical), which, therefore, cannot be related 

We impose a floor of 12% on the uncertainty of each data point, as 
explained in Section 4. 

to the passage of a synchrotron break frequency. We explore 
various explanations for this behavior in Section 5.2.4. 

3.2. X-ray/radio Steep Decline: and fjet 

The X-ray data at > 1.4d exhibit a steep decline, with 
ctx = “2.29 ±0.16. In the standard afterglow model a steep 
decline of a « -p < -2 at frequencies above I'm is expected 
after the ‘jet-break’ (fjet), when the bulk Lorentz factor, F, de¬ 
creases below the inverse opening angle of the jet, 1 /0jet, and 
the edges of the collimated outflow become visible. The r'- 
band light curve also exhibits a shallow-to-steep transition at 
Ri 1.5 d, with a post-break decay rate of or = -2.27 ±0.11, 
consistent with the X-rays. These observations suggest that 
the jet break occurs at about 1.5 d. We now consider whether 
this interpretation is consistent with the radio observations. 

The multi-wavelength radio SED at 9d is well-fit by a 
single power law with spectral index, /3 = 1/3 from 7 GHz 
to 93 GHz (Figure 2 and Section 3.4), indicating that > 
93 GHz at this time. In the absence of a jet break, we would 
expect the flux density below t'ni to remain constant (for a 
wind-like circumburst environment) or rise with time (for a 
constant density circumburst environment). However, the ra¬ 
dio light curves decline after 4.6 d at all frequencies from 
15 GHz to 93 GHz. The combination of Vm > 93 GHz at 9d 
and the declining light curve in the radio bands is only possi¬ 
ble in the standard afterglow model if a jet break has occurred 
before 4.6 d. This is consistent with the steepening observed 
in the X-ray and optical light curves at Ri 1.5 d, suggesting 
that fjet ~ 1.5 d. 

We note that the flux at a given frequency decays steeply 
following a jet break only once has crossed the observ¬ 
ing frequency; thus the steepening in the radio light curves is 
expected to be delayed past that of the steepening in the X- 
ray and optical light curves until passes through the radio 
band. We find that the 7 GHz light curve is consistent with be¬ 
ing flat to 50 d, after which it declines rapidly with aradio ~ “2. 
This suggests that crosses the 7 GHz band at around 50 d. 
Since is expected to decline as f“^ following the jet break, 
we have ~ 8 x 10*^ Hz at 1.5 d if we take fjet ~ 1.5 days. 
Thus Vjn is below both the X-rays and the optical frequen¬ 
cies at f = fjet, consistent with the steepening being observed 
around the same time in the X-ray and optical bands. 

To summarize, a simultaneous steepening in the X-ray and 
optical light curves between 1 and 2d indicates that fjet ~ 
1.5 d. Taken together with a similar steepening observed in 
the 7 GHz radio light curve around 50 d, we find f^m ~ 7 GHz 
at50d. 


Table 4 

Parameters for broken power law fit to X-rayAJV/optical re-brightening for GRB 120326A 



XRT 

t/-band 

z-band 

Jointt (MCMC) 

Break time, t\j (d) 

0.41 ±0.02 

0.31 ±0.04 

0.3lt 

0.40 ±0.02 
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Table 4 — Continued 


XRT 

Flux density at ( 5 , F\, (^Jy) 1.23 it 0.05 

Riserate, O] 0.85±0.19 0. 

Decay rate, 02 -1.22±0.18 -1 

Smoothness, y 5.0 it 4.0 


Note. — t Parameters fixed during fit (see text for details). ^ We 

used a flat prior on all parameters, with the ranges % S [0.05,1.0], Fb x S 

[1 X 10^“*. 1 X 10^2] mJy, Fb,uv £ [5 x lO^El x 10-']mjy, Fb,z S [1 x 

10-2,2]mJy,oi e [0.1,5.0],a2 S [-5.0,0.0], andy e [0.5,30]. 

3.3. The Circumburst Density Profile and the Location of Vc 

The density profile in the immediate (sub-pc scale) environ¬ 
ment of the GRB progenitor impacts the hydrodynamic evolu¬ 
tion of the shock powering the GRB afterglow. The evolution 
of the shock Lorentz factor is directly reflected in the after¬ 
glow light curves at all frequencies and multi-band modeling 
of GRB afterglows therefore allows us to disentangle different 
density profiles. Since the progenitors of long-duration GRBs 
are believed to be massive stars, the circumburst density struc¬ 
ture is expected to be shaped by the stellar wind of the pro¬ 
genitor into a profile that falls of as p(r) = Ar~^ with radius 
r. Here A = M„/4ttV„ = 5 x 10"A* gcm“’ is a constant pro¬ 
portional to the progenitor mass-loss rate M„ (assumed con¬ 
stant), for a given wind speed, y„, (Chevalier & Li 2000), and 
A* is a dimensionless parametrization, corresponding to M„ = 

1 X lO“^M 0 yr“' and V^, = lOOOkms"'. Alternatively, the 
shock may directly encounter the uniform interstellar medium 
(ISM). Both wind- and ISM-like environments have been in¬ 
ferred for different events from previous observations of GRB 
afterglows (e.g., Panaitescu & Kumar 2000; Harrison et al. 
2001; Panaitescu & Kumar 2001, 2002; Yostetal. 2002; 
Frail et al. 2003; Chandra etal. 2008; Cenko et al. 2010, 
2011; Laskaret al. 2013, 2014). Here, we explore the con¬ 
straints on the progenitor environment of GRB 120326A. 

The spectral index between the PAIRITEL K-band obser¬ 
vation at 1.4 d (Morgan 2012) and the X-rays is /3nir,x = 
-0.96 ± 0.05 (Figure 3), which is consistent with the X-ray 
spectral index of fix = -0.85 ± 0.04 (Section 2.1) at 2a, sug¬ 
gesting that the NIR, optical, and X-ray bands are located on 
the same power law segment of the afterglow SED at 1.4 d. At 
the same time, the spectral index within the NIR/optical bands 
is fi = -1.80 ±0.16, indicating that extinction is present. 

Since /?nir,x « /3x ~ -0.90, the cooling frequency, r'c must 
lie either below the NIR or above the X-rays at 1.4 d. For 
Fc < Fir, we would infer an electron energy index, p Ri 1.8 
and a light curve decay rate of a ~ -0.85 regardless of the cir¬ 
cumburst density profile. On the other hand, > fx requires 
p « 2.8 with a Ri -1.35 for a constant density environment 
and a Ri -1.85 for a wind-like environment. Actual measure¬ 
ments of the light curve decay rate at r:; 1.4 d are complicated 
by the presence of the jet break at around this time. The op¬ 
tical and X-ray light curves decline as ctr = -2.05 ±0.13 and 
ax = 2.29 ± 0.16 after the jet break, with the expected decay 
rate being aK,-p (Rhoads 1999; Sari et al. 1999). This in¬ 
dicates p K. 2 (Section 3.2) and Fc < F„ for an ISM model. 
Upon detailed investigation (Section 4), we find that a p ^2 
model with an ISM-like environment fits the data after the re¬ 
brightening well. This model additionally requires < f^ 
(fast cooling) until Ri 2 d. For completeness, we present our 


[/-band z-band Jointt (MCMC) 

66±7 255±10 — 

89±0.21 0.59±0.03 0.52±0.06 

.10 it 0.29 -0.94 ±0.05 -1.10 ±0.10 

nt nt 1 ”7 _L 1A 



Figure 2. Afterglow SED for GRB 120326A at 9 d from 5 GHz to 92.5 GHz 
together with a best fit power law fit to the data above 7 GHz. The spectrum 
is optically thin with a spectral index of iy0.35±0.03 7 to 92.5 GHz. 

The 5 GHz observation shows evidence of synchrotron self-absorption. 



Figure 3. Afterglow SED for GRB 120326A at 0.33 d from the NIR to the 
X-rays together with the best fit forward shock model (red, solid). The optical 
and UV data exhibit a clear decrement due to extinction in the host galaxy. 
The solid line is the SED from the highest-likelihood model and the dashed 
curve indicates the SED in the absence of absorption along the line of sight 
to the GRB by the Milky Way and the GRB host (Section 5.1). The spectral 
break apparent in the model SED at 5 ^ Hz is z/m- The data above 
9 X Hz (Ly-CK in GRB rest frame) are likely affected by absorption by 
the intergalactic medium (IGM) along the line of sight. 

investigation of the wind model in Appendix B. 

3.4. Location of 

The radio SED from 7 GHz to 93 GHz at 9 days is opti¬ 
cally thin with a spectral slope of /3 = 0.33 ±0.04 (Figure 2), 
suggesting that lies above 93 GHz and lies below 7 GHz 
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at this time. This is consistent with the passage of through 
7 GHz at 50 d inferred in Section 3.2. The spectral index be¬ 
tween 5 GHz and 7 GHz at 9d, ^ = 0.92 ±0.23, is steeper than 
v^/^. This may suggest that the synchrotron self-absorption 
frequency, ~ 5 GHz. However, this spectral index does not 
show a monotonic trend with time. Since this part of the ra¬ 
dio spectrum is strongly affected by interstellar scintillation 
(ISS) in the ISM of the Milky Way, a unique interpretation 
of the observed spectral index is difficult at these frequencies. 
This difficulty in constraining results in degeneracies in the 
physical parameters. We return to this point in Section 5.1. 

4. MULTI-WAVELENGTH MODELING 

Although the panchromatic peak at 0.4 d is a unique feature 
of GRB 120326A, the X-ray, optical, and radio light curves of 
this event exhibits standard afterglow features after this time, 
with evidence for an un-broken power law spectrum extend¬ 
ing from the optical to the X-rays (Section 3.3), a spec¬ 
trum in the radio (Section 3.4), and evidence for a jet break at 
fjet ~ 1.5 d. (Section 3.2). We therefore determine the physi¬ 
cal properties of this event by using observations after the X- 
ray/optical peak. We model the data at > 0.4 d as arising from 
the afterglow blastwave, using the smoothly-connected power 
law synchrotron spectra described by Granot & Sari (2002). 
We compute the break frequencies and normalizations using 
the standard parameters: the fractions of the blastwave energy 
imparted to relativistic electrons (ce) and magnetic fields (e^), 
the kinetic energy (Ek^so), and the circumburst density (no). 
We also use the SMC extinction curve" (Pei 1992) to model 
the extinction in the host galaxy (Ay)- Since the E-band light 
curve flattens at « 10 days due to contribution from the host, 
we fit for the E-band flux density of the host galaxy as an 
additional free parameter. 

The various possible orderings of the spectral break fre¬ 
quencies (e.g., Vm < Vc- ‘slow cooling’ and z/c < z/m: ‘fast 
cooling’) give rise to five possible shapes of the afterglow 
SED (Granot & Sari 2002). Due to the hydrodynamics of 
the blastwave, the break frequencies evolve with time and the 
SED transitions between spectral shapes. To preserve smooth 
light curves when break frequencies cross and the spectral 
shape changes, we employ the weighting schemes described 
in Laskar etal. (2014) to compute the afterglow SED as a 
function of time. To efficiently and rapidly sample the avail¬ 
able parameter space, we carry out a Markov Chain Monte 
Carlo (MCMC) analysis using a python implementation of the 
ensemble MCMC sampler EMCEE (Eoreman-Mackey et al. 
2013). Eor a detailed discussion of our modeling scheme, 
see Laskar etal. (2014). To account for heterogeneity of 
UV/Optical/NIR data collected from different observatories, 
we usually impose an uncertainty floor of 5% prior to fitting 
with our modeling software. Eor this GRB, we find that the 
fit is driven by the optical data at the expense of fits at the 
radio and X-ray bands, which we mitigate by increasing the 
uncertainty floor further to 12%. We additionally correct for 
the effect of inverse Compton cooling (Appendix A). 

5. RESULTS EOR GRB 120326A 

5.1. Multi-wavelength model at > 0.4 d 

In confirmation of the basic analysis presented in Section 
3.3, our highest likelihood model (Eigure 4) has p « 2.09, 

* * Our previous work shows negligible changes in the blastwave parame¬ 
ters with LMC and Milky Way-like extinction curves (Laskar et al. 2014). 


Ce ~ 0.33, €b ~ 0.33, iiQ ^ 0.27 cm“^, Ek.Iso ~ 1.4 x 10^^ erg, 
fjet ~ 1.5 d. Ay « 0.48 mag, and an r'-band flux density of 
2.3 /iJy for the host galaxy. The spectrum is in fast cooling un¬ 
til 1.4 d. During the fast cooling phase, z/a splits into two dis¬ 
tinct frequencies: z/jc and Usa (Granot & Sari 2002). The spec¬ 
trum has the Rayleigh-Jeans shape below z^ac, and a slope 
of z/"/^ between r'ac and Eor the highest likelihood model, 
the break frequencies^^ are located at z/ac ~ 2.8 x 10® Hz, 
z/sa « 4.9 X 10® Hz, z^c « 9.3 x 10*2, and z^n, « 1.5 x 10'^ Hz 
at 1 d and the peak flux density is Ri 18 mJy at z/m evolves 
as f“2 following the jet break at r; 1.5 d to r 7.7 GHz at 50 d, 
consistent with the basic considerations outlined in Section 

3.2. The Compton y-parameter is r 0.6, indicating that cool¬ 
ing due to inverse-Compton scattering is moderately signifi¬ 
cant. 

We present histograms of the marginalized posterior density 
for each parameter in Eigure 5. Most of our radio observa¬ 
tions are after the transition to slow cooling, at which time z/^ 
lies below the lowest radio frequency observed, resulting in 
degeneracies between the model parameters (Eigure 6). We 
summarize the results of our MCMC analysis in Table 5. 

Using the relation 


EK,iso,52\ * /fjet/(l+z)\ 

no J 1 1 ^ J 


( 2 ) 


for the jet opening angle (Sari et al. 1999), and the distribu¬ 
tions of Ekaso, no, and fjet from our MCMC simulations, we 
And Hjsi — 4.6 ± 0.2 degrees. Applying the beaming correc¬ 
tion, E = £iso(l -cos0jet), we And £-, = (1.0 ± 0.1) x 10^® erg 
(l-lO^^keV; rest frame) and Ek = (4.6^g j) x 10^® erg. 

5.2. X-ray/UV/optical Re-brightening 

We now consider physical explanations of the unusual re¬ 
brightening between 0.1 and 1 day observed in the X-ray, UV, 
and optical bands. 

Scattering by dust grains in the host galaxy of GRBs has 
been suggested as a potential explanation for the shallow de¬ 
cay phase of X-ray light curves (Shao & Dai 2007). However, 
this mechanism is expected to cause a significant softening of 
the X-ray spectrum with time (Shen et al. 2008), and cannot 
produce light curves that rise with time, as is observed in the 
case of GRB 120326A at about 0.4 d. The bump can also 
not be caused by the passage of a spectral break frequency, 
since it occurs almost simultaneously in the X-ray, UV, and 
optical bands. Eour remaining potential models for the bump 
are: (1) onset of the afterglow, (2) geometric effects due to an 
observer located outside the jet (off-axis scenario), (3) a den¬ 
sity enhancement in the circumburst environment, and (4) a 
refreshed shock (energy injection scenario). We now explore 
these possibilities in turn. 


5.2.1. Afterglow onset 

In this scenario, the peak of the light curve emission cor¬ 
responds to emission from a reverse shock due to deceler¬ 
ation of the ejecta by the surrounding material. The ejecta 
are assumed to be composed of a conical shell segment at 
a single Lorentz factor. The reverse shock light curves be¬ 
fore and after the deceleration time depend on the properties 


Uc here is reported as where i /3 and u\\ are (differently- 

normalized) expressions for the cooling frequency in the slow cooling and 
fast cooling regimes, respectively (Granot & Sari 2002), such that = 

Wefindi /3 = 1.2 X 10^^ and t^n = 7.4 x 10 ^^ 1 d. 
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Figure 4. X-ray, UV (top left), optical (top right), and radio (bottom left) light curves of GRB 120326A in the ISM scenario, with the full afterglow model (solid 
lines), including energy injection before 0.4 d. The X-ray data before 0.004 d is likely dominated by high-latitude prompt emission and we do not include these 
data in our analysis. The dashed envelopes around the radio light curves indicate the expected effect of scintillation at the la level. The data prior to the end of 
the re-brightening at 0.4 d (open symbols) are not used to determine the parameters of the forward shock in the MCMC analysis. The Sw/ft/UVOT data in the 
uvw2 and uvm2 bands are strongly affected by IGM absorption and are not included in the analysis. Bottom right: blastwave Lorentz factor (green, dashed; upper 
sub-panel) and isotropic equivalent kinetic energy (red, solid; upper sub-panel) as a function of time, together with the energy distribution across ejecta Lorentz 
factors (black, solid; lower sub-panel) as determined from fitting the X-rayAJV/optical re-brightening at 0.4 d. 



0.34 0.36 0.38 0.40 0.42 0.44 0.46 


8.0 xlO^ 
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0.270 0.285 0.300 0.315 0.330 







Figure 6. Icr (red), 2cr (green), and Scr (black) contours for congelations be¬ 
tween the physical parameters, £'K,iso> and eg for GRB 120326A, in 

the ISM model from Markov chain Monte Carlo simulations, together with 
the maximum-likelihood model (blue dot). We have restricted Ce < 

< V^- on-line version of this Figure for additional plots of corre¬ 

lations between these parameters and p, tjet, ^jet. Fk, Ay, and Fi^^hosi,R- 


Figure 5. Posterior probability density functions for the physical parameters 
of GRB 120326A in the ISM model from MCMC simulations. We have 
restricted Ce < Vs and es < V^- 


of the environment (ISM or wind) and the ejecta. Reverse 


shock light curves for a constant density environment depend 
upon whether the ejecta shell is thick (Aq > or thin 

(Aq > where Aq is the initial shell width, Fq is the 

initial Lorentz factor of the ejecta and / = (3£’/47rnompC^)^/^ 
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is the Sedov length (Kobayashi 2000). Similarly, when the 
circumburst medium has a wind-like profile, the light curves 
again depend on whether Ao > £/ 47 rAc^ro (thick shell) or 
Ao < E / 47 rAc^rQ (thin shell; Kobayashi & Zhang 2003). 

In the thick shell case, the reverse shock crosses the ejecta 
in a time comparable to the burst duration, t^, and the light 
curves at all frequencies are expected to peak at that time 
scale. In the case of GRB 120326A, Tgo in the S'wi/f/BAT 
band is ~ 27 s, which is much shorter than the observed time 
at which the light curves peak, fpeak ~ 0.4 d. Thus, the re¬ 
brightening cannot be explained by reverse shock emission in 
the thick shell case. 

In the thin shell case, temporal separation is expected be¬ 
tween the GRB itself and the peak of the reverse shock, 
which occurs when the RS crosses the shell at = 

90£’j^'^;jq 52 '^o'^^royoo ® (ISM environment; Kobayashi 2000) 
or t = 2.9 X 10^(1 +z)£’K,iso, 52 ro%oA“'s (wind environment; 
Zou et al. 2005). A reverse shock peak at 0.4 d would imply 
a rather low burst Lorentz factor of Fq ^ 17. Additionally, 
the light curves at all frequencies would be expected to rise 
until t.y and decline thereafter (except above where the flux 
is constant until and disappears after f^-)- The observed R- 
band light curve, on the other hand, declines between 100 s 
and 500 s, followed by a flat segment until 3500 s, followed 
by the rise into the bump at 0.4 d. This combination of a flat 
portion followed by a rise cannot be explained as being due to 
a reverse shock. 

Thus the onset of the afterglow does not provide a viable ex¬ 
planation for the re-brightening regardless of density profile, 
and whether the ejecta are in the thick- or thin-shell regime. 

5.2.2. Ojf-axis model 

We now investigate the re-brightening in the context of 
viewing geometry effects. We consider a jet with an open¬ 
ing half-angle, 0jet, viewed from an angle 0obs- While small 
offsets in the observer’s viewing angle relative to the jet axis 
do not cause significant changes in the light curves provided 
that 0 obs ^ djet, it is possible to obtain a rising light curve when 
^obs ^ 20jet (Granot et al. 2001, 2002). In this case, the time 
of the peak of the light curve, fp is related to the jet break time 
for an on-axis viewer as fp « (5-l-21n0)0^fjet > 5fjet, where 
0 = 0obs/0jet“ 1^1 (Granot et al. 2002). This implies that 
for a light curve that peaks at r; 0.4 d, the on-axis jet break 
time must have occurred earlier, at Ri 2 hours. 

In this case, the radio light curves are expected to rise un¬ 
til the observer’s line of sight enters the beaming cone of 
the jet and then approximately converge with the predicted 
on-axis light curves after 0.4 d. After the jet break, the flux 
density declines as f“‘/^ for ly < and as f“^ for v > Vm 
(for an on-axis observer; Sari et al. 1999). Thus, the radio 
light curves should be declining at all frequencies after 0.4 d. 
However, the flux density at the 15.75 GHz AMI-LA band 
rises from (0.34±0.14)mJy at 0.31d to (0.77 ± 0.08) mJy 
at 7.15d, which is inconsistent with this expectation. This 
argues against the off-axis model as an explanation for the 
X-ray/UV/optical re-brightening. 

5.2.3. Density enhancement 

If the blastwave encounters an enhancement in the local 
density as it propagates into the circumburst environment, an 
increase in the flux density is expected. Bumps in afterglow 
light curves have been ascribed to this phenomenon in the past 


(e.g., Berger et al. 2000; Lazzati et al. 2002; Schaefer et al. 
2003). However, the flux density at frequencies above z/c (for 
slow-cooling spectra) and above ly^ (for fast cooling spec¬ 
tra) are insensitive to variations in the circumburst density 
(Nakar & Granot 2007). In the case of GRB I20326A, a re¬ 
brightening is observed both in the UV/optical and in the X- 
rays. The X-rays are located above both ly^ and r'u, in the 
p Ki 2.1 (fast cooling) model. Thus in our favored afterglow 
model, the X-ray re-brightening cannot be due to a density 
enhancement. 

5.2.4. Energy injection model 

An alternative model for a re-brightening of the afterglow 
is the injection of energy into the blastwave shock due to 
prolonged central engine activity, deceleration of a Poynt- 
ing flux dominated outflow, or the presence of substantial 
ejecta mass (and hence kinetic energy) at low Lorentz fac¬ 
tors (Dai (fe Lu 1998; Rees & Meszaros 1998; Kumar & Piran 
2000; Sari & Meszaros 2000; Zhang & Meszaros 2001, 2002; 
Granot & Kumar 2006; Zhang et al. 2006; Dall’Osso et al. 
2011; Uhm et al. 2012). This mechanism has been invoked to 
explain plateaus in the observed light curves of a large fraction 
of GRB X-ray afterglows (Nousek et al. 2006; Liang et al. 
2007; Yu & Dai 2007; Margutti et al. 2010b; Hascoetetal. 
2012; Xin et al. 2012). In this section, we explore the 
energy-injection model as a possible explanation for the X- 
ray/UV/optical re-brightening at 0.4 d. 

Models involving a transfer of energy from braking radia¬ 
tion of millisecond magnetars into the forward shock predict 
plateaus in X-ray light curves but do not generally lead to a re¬ 
brightening. In particular, they require an injection luminos¬ 
ity, Locf^ (corresponding to an increase in the blastwave en¬ 
ergy with time as £ oc f *''‘‘^), with o' < 0. In our p « 2.1 model, 
the X-rays are located above the cooling frequency. In this 

regime, the flux density above ly^ is £,y>i/, oc 
(Granot & Sari 2002). For p « 2.1, this reduces to £i/>i/c OC 
■ During energy injection, £ oc f'”, such that £y>jy^ oc 
^i.03m-i.i 'j’jjg steep rise {a\ =0.85 ±0.19; Section 3.1) 
requires m = 1.89 ± 0.18 or q = m-\ = 0.89 ±0.18. En¬ 
ergy injection due to spin-down or gravitational wave radi¬ 
ation from a magnetar is expected to provide at best a con¬ 
stant luminosity {q < 0). Thus the observed re-brightening 
cannot be explained by energy injection from a magnetar. 
In the case of energy injection due to fail-back accretion 
onto a black hole, the expected accretion rate is M cx 
(Phinney 1989), which is also insufficient to power the ob¬ 
served plateaus. Similarly, central engine activity is usu¬ 
ally associated with flaring behavior in X-ray (Burrows et al. 
2005a; Fan & 'Wei 2005; Falcone et al. 2006; Proga &. Zhang 
2006; Perna et al. 2006; Zhang et al. 2006; Chincarini et al. 
2007, 2010; Margutti et al. 2010a; Bernardini et al. 2011; 
Margutti et al. 2011a) and optical (Li et al. 2012) light curves, 
but it involves shorter characteristic time scales (Af/f ^ 1) 
than observed for GRB 120326A at 0.4 d. 

In the standard afterglow model, the ejecta are assumed 
to have a single Lorentz factor. We now relax this assump¬ 
tion and consider models including a distribution of Lorentz 
factors in the ejecta as a possible explanation for the late re¬ 
brightening in GRB 120326A. Provided they are released over 
a time range small compared to the afterglow timescale, the 
ejecta arrange themselves in homologous expansion with the 
Lorentz factors monotonically increasing with distance from 
the source (Kumar &. Piran 2000). We follow the formal- 
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ism of Sari & Meszaros (2000), where the ejecta are assumed 
to possess a continuous distribution of Lorentz factors such 
that a mass, M(> E) cx E"* is moving with Lorentz factors 
greater than E with corresponding energy, E{> E) oc 
down to some minimum Lorentz factor, Emm (also known as 
the ‘massive-ejecta model’). Additionally, we posit a maxi¬ 
mum Lorentz factor for this distribution, Emax, corresponding 
to the Lorentz factor of the blastwave Ebw at the onset of 
energy injection. In this model, there is a gap between the 
blastwave, which is powered by an initial shell of fast-moving 
ejecta, and the leading edge of the remaining ejecta travel¬ 
ing at Emax- When these slower ejecta shells catch up with 
the blastwave at a time when Ebw(0 ~ ^max, they begin de¬ 
positing energy into the blastwave, and energy injection com¬ 
mences. This proceeds until the lowest energy ejecta located 
at Emin have transferred their energy to the blastwave, and the 
subsequent evolution proceeds like a standard afterglow pow¬ 
ered by a blastwave with increased energy. The Lorentz fac¬ 
tor of the blastwave at the end of energy injection is therefore 
Tmin, which Can be determined from modeling the subsequent 
afterglow evolution and invoking the standard hydrodynami- 
cal framework. 

Lor an ambient medium with a density profile, n oc r~^, the 
energy of the blastwave increases as £ oc during the pe¬ 
riod of energy injection, where*^ m = Lor a constant 

density medium (k = 0) this translates to m = which is 

bounded (0 < m < 3 for s S [1, oo)) while for a wind-like en¬ 
vironment (k = 2) we have m= which is also bounded 
(0 < m < 1 for s e [l,oo)). 

Lor GRB 120326A, the optical and X-ray light curves are 
located above both z/m and and the spectrum is fast cool¬ 
ing (z/c < i^m', Section 5.1). In this regime, the flux den¬ 
sity depends on the blastwave energy and time as, /v cx 
]7(2+p)/4f(2-3p)/4_ Writing E = Eoit/to)'” for the energy injec¬ 
tion episode, this implies q(. f [(2+/j)»!+2-3p]/4 ^ = 2. 

Since the optical and X-ray light curves rise with temporal in¬ 


dex, a R4 0.6 (Section 3.1) and p 2 (Section 5.1), this im¬ 
plies m « 1.6, which corresponds to i ~ 10. 

Lor a detailed analysis, we use the afterglow parameters for 
the maximum likelihood model listed in Section 5.1 as our 
starting model (which explains all available multi-wavelength 
data after the re-brightening). We then reduce the energy at 
earlier times as a power law with time, 

LK.iso.f, t > to 

m 

, ti <t <to (3) 

LK,iso,i = Lk iso.f 

where isK.iso.i and iiK.iso.f are the initial and final energy of 
the blastwave, respectively, with the latter fixed to the value 
determined by our MCMC analysis of the data at > 0.4 d. En¬ 
ergy injection into the blastwave at the rate E{t) cx f'" begins 
at t\ and ends at to, with to, t\, and m being free parameters 
in this model. We compute the spectrum at each time accord¬ 
ing to the instantaneous energy in the blastwave, adjusting the 
parameters to and t\ by hand until the resulting light curves 
match the observations. 

We find that the X-ray and optical light curves can be 
modeled well by a single period of energy injection, with 
LK.iso = 1.4 X 10^^erg(f/0.38d)* ‘^‘^ between fi = 4.3 x 10“^d 
and to = 0.38 d (Ligure 4). In this model, the blastwave energy 
at the start of energy injection is LK.iso.i sa 5.9 x 10^' erg. Thus 
« 96% of the final energy of the blastwave is injected during 
this episode. In comparison, E^^i^o = (3.15 ±0.12) X 10^2erg 
(Section 3). The blastwave Lorentz factor can be computed 

from the expression E = 3.65£'K.iso,52^^^«o*^^f~^^^(l + 2)^^* 
(Granot& Sari 2002) and is E Ri 19 at the start of energy 
injection, decreasing to E Ri 13 at 0.38 d. Interpreted in the 
context of the massive-ejecta model, the energy injection rate 
of ZM = 1.44 would imply s = 8.4, corresponding to an ejecta 
distribution with E{> E) cx E'“^ cx E"’ ^^ for 13 < E < 19. 


LK,iso(0 ~ ^ 


C,iso,f 


Table 5 

Energy injection in GRB afterglows 


GRB 

100418A 

100901A 

120326A 

120404A 

Redshift, ^ 

0.6235 

1.408 

1.798 

2.876 

Tgo (s) 

8.0 ±2.0 

439 ±33 

11.8±1.8 

38.7±4.1 

£^ 7 ,iso (erg; 1-10^ keV, rest frame) 

9.9!^ l X lO^o 

(8±l)x 10^2 

(3.2±0.1)x 10^2 

(9±4)x 10^2 



Best fit model* 


section in text 

6.1.2 

6.2.2 

5.1 

6.3.2 

circumburst environment 

ISM 

ISM 

ISM 

ISM 

p 

2.14 

2.03 

2.09 

2.06 

Ce 

0.12 

0.33 

0.33 

0.27 

eg 

1.1 X 10“^ 

0.32 

0.33 

0.16 

no (cm“3) 

1.4 

3.2 X 10-2 

0.27 

2.8 X 102 

£K,iso(10^^erg) 

3.36 

29.7 

13.7 

12.3 

Av (mag) 

<0.1 

0.1 

0.40 

0.13 

Zjet (days) 

16.9 

0.96 

1.54 

6.6 X 10-2 

(deg) 

20.4 

2.1 

4.7 

3.1 

£K(10®erg) 

21.1 

2.1 

4.5 

1.7 
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Other authors have defined this in terms of luminosity, with Loct or 


E ^. Thus m = 1-q. 
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Table 5 — Continued 


GRB 

100418A 

100901A 

120326A 

120404A 

erg; l-lO'^keV, rest frame) 

0.62 

0.5 

1.0 

1.3 

Peak time of re-brightening (d) 

0.75±0.13 

0.36±0.02 

0.40 ±0.01 

(2.8 ± 0.4) X 10-3 

^K.iso,initial 

1.0 X 10^“ 

5.4 X 10^2 

3.1 X 103‘ 

4.4 X 103‘ 

^K.initial (erg) 

6.3 X lO"** 

3.8 X lO"*^ 

8.7 X lO'^* 

6.4 X lO"** 



MCMC Results* 


p 

o 1 Q-t-0.02 
^••^■^-0.01 

2.05 ±0.01 

2.095 ±0.007 

2.07 ±0.02 


1.2)|j;3 X 10-‘ 

-0.05 

n 'T9Q+0.003 
u.jz^-0_oo7 

o 

to 

bb 

eg 

9.4+” X 10^3 

-0.08 

^•■^^-0.03 

0 

^•^■^-0.08 

no(cm“^^) 

-0.3 

7.2)|j X 10-3 

0 27^-^ 

^•"^'-0.03 

3.5!^| X 103 

log(no)t 

0.22 ±0.09 

—2 

^•^-0.2 

-0.57 ±0.06 

^•^-0.3 

£k,iso (10^2 erg) 

J.D_0 7 

31 

14.0 ±0.07 

1 j.j-2 0 

Av (mag) 

<0.1 

0.09±0.01 

0.40 ±0.01 

0.13±0.01 

tjet (days) 

17.3±1.0 

0.97 ±0.02 

1.55±0.06 

(6.9 ± 0.8) X 10-3 

^'jet (deg) 

20.9 ±0.5 

2.4 ±0.1 

4.6 ±0.2 

3.1±0.3 

£K(105'’erg) 


2 

^•°-0.4 


2 0"^-^ 

£-,(10^® erg; 1-10** keV, rest frame) 

0.66 ±0.33 

0.70±0.10 

1.0±0.10 

1.3±0.6 


Note. — * Note that the parameters for the best-fit model may differ 
slightly from the results from the MCMC analysis. This is because the for¬ 
mer is the peak of the likelihood distribution (and is appropriate for generat¬ 
ing model light curves), while the latter is summarized here as 68% credible 
intervals about the median of the marginalized posterior density functions 
for each parameter. The posterior density serves as our best estimate for 
the value of each parameter, and may be asymmetric about the value corre¬ 
sponding to the highest-likelihood model. 

t In instances where the measured value of mq spans more than about a factor 
of 2, log no is a more meaningful quantity. We therefore report both no and 
log no for all cases. 

6 . PANCHROMATIC RE-BRIGHTENING EPISODES IN 
OTHER GRBS 

Whereas individual flattening or re-brightening episodes 
have been seen both in optical and X-ray observations of GRB 
afterglows (Liang et al. 2007; Li et al. 2012), simultaneous re¬ 
brightening of the afterglow in multiple bands spanning both 
the optical and the X-rays as seen in GRB 120326A is quite 
rare. The X-ray and R-band light curves of GRB 970508 
exhibited a bump at around 60 ks, lasting for about lOOks 
(Piro et al. 1998). Multiple episodes of flux variations of 
~ 30% on time scales of At/t k. 1 between 1 and 10 d since 
the burst were detected superposed on a power law decline 
in multi-band X-ray and BVRI-bdsvA data for GRB 030329 
(Lipkin et al. 2004), while de Ugarte Postigo et al. (2007) de¬ 
tected an X-ray and optical re-brightening in GRB 050408 
around 260ks. The short burst GRB 050724 had a large 
re-brightening at around 42 ks, simultaneous with an optical 
flare (Berger et al. 2005; Malesani et al. 2007; Margutti et al. 
2011b). GRB 060206 exhibited a simultaneous X-ray 
and optical re-brightening at 3.8ks (Monfardini et al. 2006; 
Liu et al. 2008). GRB 070311 exhibited a re-brightening at 
Ri looks in the X-rays and in R-band, but the peak in the 
X-rays is not well-sampled and does not appear to have oc¬ 
curred at the same time in the two bands (Guidorzi et al. 
2007). Covino et al. (2008) report a simultaneous X-ray and 
RJHK-bdxvA re-brightening for GRB 071010A at « 52 ks. 
GRB 060614 exhibited a prominent brightening in the UBVR- 


bands simultaneous with an X-ray plateau (Mangano et al. 
2007). Similarly, the z Ri 6 GRB 120521C exhibited an X-ray 
plateau at the same time as a bump in the z-band light curve 
(Laskar et al. 2014). GRB 081028 was one of the first events 
with a well-sampled rising X-ray light curve following the 
steep decay phase with a concomitant optical re-brightening 
(Margutti et al. 2010b). An X-ray re-brightening simultane¬ 
ous with a re-brightening in Rvw/f/UVOT observations was 
seen for GRB 100418A and interpreted as energy injection 
by Marshall et al. (2011). Gorbovskoy et al. (2012) report 
a prominent multi-band re-brightening in GRB 100901 A, 
whose X-ray light curve is remarkably similar to that of 
GRB 120326A. GRB 120404 was found to exhibit a strong 
re-brightening in the optical and NIR, although X-ray obser¬ 
vations around the time of the re-brightening are sparse as this 
time fell within an orbital gap of Swift (Guidorzi et al. 2013). 

Of these instances of X-ray or optical re-brightening events 
seen in GRB afterglow light curves, we select those objects 
with multi-band (at least X-ray and optical) datasets that 
exhibit a simultaneous re-brightening in both X-rays (fol¬ 
lowing the steep decay phase) and optical/NIR. Since ra¬ 
dio observations are vital for constraining the physical pa¬ 
rameters, we restrict our sample to events with radio detec¬ 
tions; GRBs 100418A, 100901A, and 120404A. We perform 
a full multi-wavelength afterglow analysis for these events 
and compare our results with those for GRB 120326A in Sec¬ 
tion 7. 

6.1. GRB 100418A 

6.1.1. GRB properties and basic considerations 

GRB 100418A was detected and localized by the Swift 
BAT on 2010 April 18 at 21:10:08 UT (Marshall et al. 2010). 
The burst duration is Tgo = 8 .0 ± 2.0 s, with a fluence of 
Rly = (3.4 ±0.5) X 10“^ergcm“^ (15-150keV observer frame; 
Marshall et al. 2011). The afterglow was detected in the X- 
ray and optical bands by XRT, UVOT, and various ground- 
based observatories (Eilgas et al. 2010; Updike et al. 2010c; 








Energy injection in GRBs 


15 


Siegel & Marshall 2010), as well as in the radio by WSRT, 
VLA, the Australia Telescope Compact Array (ATCA), and 
the Very Long Baseline Array (VLBA) (van der Horst et al. 
2010a; Chandra & Frail 2010b; Moin et al. 2010). Spec¬ 
troscopic observations with the ESO Very Large Telescope 
(VET) yielded a redshift of z = 0.6235 (Antonelli et al. 
2010). By fitting the 7 -ray spectrum with a Band function, 
Marshall et al. (2011) determined the isotropic equivalent 7 - 
ray energy of this burst to be E-^.iso = ( 9 . 9 ;|; 34 ) x 10 ^*’ erg ( 1 - 
lO^keV; rest frame). 

The X-ray light curve for this burst gently rises to 0.7 d, 
while the optical light curves exhibit a slow brightening at 
the same time. Both the X-ray and optical light curves break 
into a power law decline following the peak. Using XRT and 
UVOT W/u'fe-band observations, Marshall et al. (2011) find 
that the X-ray and optical bands are located on the same part 
of the synchrotron spectrum after the peak at 0.7 d. By fit¬ 
ting the X-ray and optical decline rates and X-ray-to-optical 
SED, they find that both and are located below the op¬ 
tical, requiring p 2.3. In this regime, the light curves are 
insensitive to the density profile of the circumburst environ¬ 
ment. In the absence of radio data, Marshall et al. (2011) as¬ 
sume an ISM model and determine that the non-detection of 
a jet break out to 23 d requires a high beaming corrected en¬ 
ergy, Ek > 3 X 54 ^* 0 ^^- They investigate the off-axis 

model, and compute an intrinsic burst duration for an on-axis 
observer of Eg™ = 'D{9 = Oobs— ^jet)/I2(6> = OjTgo sa 10 ms, where 
V = [ 70(1 -/3ocos0)]“' is the Doppler factor, 6obs is the off- 
axis viewing angle, 0 jet is the jet opening half-angle, and Tgo 
is the observed burst duration. Based on this short on-axis 
Eg™ and large intrinsic energy output, they argue against the 
off-axis scenario. 

Marshall et al. (2011) also consider the energy injection 
model. For the magnetar model involving injection of energy 
into the blastwave via Poynting flux, they And q = -0.23^g [4, 
which is marginally consistent with the q = 0 case. For mod¬ 
els with a Lorentz factor distribution, they use their measured 
value of q to calculate s = 6.4j;[ g (ISM model), or i < 0 (wind 
model, which is unphysical). They argue that since the ISM 
value is too different from the value i « 2.5 indicated by X- 
ray plateaus in Swift observations (Nousek et al. 2006), this 
mode of energy injection is unlikely. 

Moin et al. (2013) reported multi-wavelength radio obser¬ 
vations of GRB 100418A from 5.5 to 9.0GHz from ATCA, 
VLA and VLBA. Using upper limits on the expansion rate of 
the ejecta derived from the VLBA observations, they report 
an upper limit on the average ejecta Lorentz factor of L < 7 
and an upper limit on the ejecta mass ofMej <0.5 x IO^Mq. 
They use this limit on L to suggest that the contribution of 
low-Lorentz factor ejecta to the blastwave energy must be 
negligible during the period of energy injection, and that a 
separate injection mechanism is required. We include their 
radio observations in our analysis (Section 6.1.2) and address 
the constraints derived from the VLBA data in Section 6.1.3. 

Marshall et al. (201 1) fit multi-band UVOT and XRT spec¬ 
tra at three different epochs (7 x 10<, 0.8, and 7.3 d) and 
found that a single power law with spectral index /3 Ri -1.0 
fit all three epochs well. We follow the procedures outlined 
by Evans et al. (2007), Evans et al. (2009), and Margutti et al. 
(2013) to obtain time-resolved XRT spectra for this burst. The 
spectra are well fit by an absorbed single power law model 
with Vn^Gai = 4.78 x 10^*^cm< (Kalberla et al. 2005), and in¬ 


trinsic hydrogen column Anint = x 10 ^'cm with 

r = 3.35 ±0.09 and L = 3.9 ±0.3 during the early (85-130s; 

= 1.74 for 123 degrees of freedom) and late (130s-180s; 

= 0.97 for 55 degrees of freedom) WT-mode steep decay, 
respectively, followed by Lx = 2.08 ± 0.18 during the plateau 
phase (2.3 x 10<d-1.8d; x^ = 0-75 for 79 degrees of free¬ 
dom), and Lx = 1.81 ± 0.20 for the remainder of the observa¬ 
tions (1.8-34.5 d; x^ = 0-70 for 89 degrees of freedom). We 
obtain a flux light curve in the 0.3-10keV XRT band follow¬ 
ing the methods reported by Margutti et al. (2013), which we 
convert to a flux density light curve at 1 keV using the above 
time-resolved spectra. We also analyze all UVOT photometry 
for this burst and report our results in Table 6 . 


Table 6 

Swift/UVOT Obsei-vations of GRB 100418A 


?-?0 

(days) 

Filter 

Frequency 

(Hz) 

Flux density 
(AiJy) 

Uncertainty 

(/iJy) 

Detection? 

(1 = Yes) 

0.00188 

White 

8.64e+14 

11.5 

2.47 

1 

0.00499 

White 

8.64e+14 

10.6 

2.81 

1 

0.00591 

ii 

8.56e+14 

8.41 

2.74 

0 

0.00753 

b 

6.92e+14 

74.7 

22.3 

0 

0.011 

White 

8.64e+14 

10.8 

3.25 

1 

0.0119 

V 

5.55e+14 

123 

35.3 

0 

0.0125 

uvwl 

1.16e+15 

19.9 

5.85 

0 

0.0126 

uvw2 

1.48e+15 

12.9 

3.62 

0 

0.0152 

u 

8.56e+14 

50.9 

14.8 

0 

0.0155 

b 

6.92e+14 

101 

31.5 

0 

0.0158 

White 

8.64e+14 

22.8 

7.49 

0 

0.0441 

uvm2 

1.34e+15 

10.1 

3.15 

0 

0.0465 

uvm2 

1.34e+15 

9.69 

2.95 

0 

0.0604 

V 

5.55e+14 

67.3 

19.7 

0 

0.0699 

b 

6.92e+14 

34.5 

10.2 

1 

0.0711 

uvm2 

1.34e+15 

7.29 

2.32 

1 

0.0723 

White 

8.64e+14 

15.8 

3.22 

1 

0.0735 

uvwl 

1.16e+15 

8.37 

2.71 

0 

0.0746 

uvw2 

1.48e+15 

7.68 

2.23 

0 

0.0757 

u 

8.56e+14 

18.3 

5.7 

1 

0.077 

V 

5.55e+14 

67 

21.9 

0 

0.602 

White 

8.64e+14 

33.5 

4.57 

1 

1.03 

White 

8.64e+14 

22.6 

3.15 

1 

1.77 

White 

8.64e+14 

11.7 

1.68 

1 

2.8 

White 

8.64e+14 

8.74 

1.3 

1 

3.45 

White 

8.64e+14 

5.67 

0.829 

1 

4.42 

White 

8.64e+14 

6.22 

0.951 

1 

5.22 

uvwl 

1.16e+15 

6.28 

1.76 

0 

5.22 

11 

8.56e+14 

13 

3.96 

0 

5.22 

b 

6.92e+14 

27.5 

8.56 

0 

5.23 

uvw2 

1.48e+15 

2.87 

0.819 

0 

5.23 

V 

5.55e+14 

60.2 

18.9 

0 

5.23 

uvm2 

1.34e+15 

2.92 

0.967 

0 

5.43 

White 

8.64e+14 

4.05 

0.735 

1 

7.64 

uvwl 

1.16e+15 

3.93 

1.28 

0 

7.64 

u 

8.56e+14 

5.55 

1.84 

0 

7.64 

b 

6.92e+14 

12.1 

4.02 

0 

7.64 

White 

8.64e+14 

3.74 

1.13 

1 

7.64 

uvw2 

1.48e+15 

2.98 

0.995 

0 

7.64 

V 

5.55e+14 

38.2 

12.5 

0 

7.65 

uvm2 

1.34e+15 

4.05 

1.27 

0 

8.64 

uvwl 

1.16e+15 

3.29 

1.06 

0 

8.64 

u 

8.56e+14 

5.05 

1.64 

0 

8.64 

b 

6.92e+14 

11.3 

3.68 

0 

8.65 

White 

8.64e+14 

2.89 

1.05 

0 

8.65 

uvw2 

1.48e+15 

2.61 

0.79 

0 

8.65 

V 

5.55e+14 

31.4 

10.3 

0 

8.65 

uvm2 

1.34e+15 

3.83 

1.14 

0 

9.58 

White 

8.64e+14 

2.89 

0.557 

1 

10.4 

White 

8.64e+14 

2.73 

0.482 

1 

11.4 

White 

8.64e+14 

2.61 

0.51 

1 

12.4 

White 

8.64e+14 

4.04 

0.635 

1 

13.4 

White 

8.64e+14 

2.74 

0.468 

1 

14.4 

White 

8.64e+14 

3.2 

0.545 

1 

15.4 

White 

8.64e+14 

3.11 

0.525 

1 

16.9 

White 

8.64e+14 

2.28 

0.425 

1 

17.6 

White 

8.64e+14 

3.63 

0.587 

1 







Figure 7. Broken power law fit to the X-ray (grey, solid) and Swiftl\]\OI 
White-hand light curves for GRB 100418A near the re-brightening ai'ound 
0.6 d. The X-ray fit includes points between 0.02 and 50 d, while the White- 
band fit includes points between 6 X 10“^ and 3.1 d. A joint fit yields the 
best fit parameters fb = O.SS^gj^ d, a[ = 0.42;^'jj, and a 2 = —l.lO^'ji" 
a fixed value of y = 5.0 (Section 6.1.1). 


Table 6 — Continued 


t-to 

(days) 

Filter 

Frequency 

(Hz) 

Flux density 

(fiJy) 

Uncertainty 

(/iJy) 

Detection? 

(1 = Yes) 

19.7 

White 

8.64e+14 

2.55 

0.464 

! 

22.4 

White 

8.64e+14 

2.07 

0.502 

1 

23.5 

White 

8.64e+14 

3.5 

0.62 

1 

24.8 

White 

8.64e+14 

2.79 

0.587 

1 

25.6 

White 

8.64e+14 

3.01 

0.591 

1 

26.7 

White 

8.64e+14 

1.82 

0.594 

1 

34.4 

White 

8.64e+14 

1.97 

0.401 

1 

The X-ray light curve up to 8 x 10 

d declines steeply as 


ax = -4.52 ± 0.14, transitioning to a plateau, at about 0.01 d, 
followed by a re-brightening. The X-ray observations be¬ 
tween 0.02 d and 5.5 d can be well fit with a broken power 
law model, with fb = 0.54 ± 0.18 d, /v x(fb) = 0.32 ±0.08 pJy, 
a\ = 0.46 ±0.20, and a^ = -1.15 ±0.10, with the smooth¬ 
ness parameter fixed at y = 5.0. A broken power law fit to the 
UVOT White-bdXiA data between 6 x 10“^ d and 3.1 d yields 
fb = 0.64 ± 0.17d, 7v.white(fb) = 32 ±4/xJy, a\ = 0.34 ± 0.09, 
and a 2 = -1.02 ± 0.18, with y = 5.0 (Figure 7). The decline 
rate in the White band following the re-brightening is con¬ 
sistent with the decline rate in the /?-band between 0.4 d 
and 10 d (ctr = -0.97 ± 0.08) and in the //-band between the 
PAIRITEL observations at 0.46 d and 1.49 d (cth = -1.0 ± 
0.1). Thus, the break time and rise and decay rates of the 
re-brightening in the X-ray and optical bands are consistent. 
Finally, we fit the X-ray and UVOT Mi/fe-band data jointly, 
where we constrain the model light curves in the two bands to 
have the same rise and decay rate and time of peak with inde¬ 
pendent normalizations. Using a Markov Chain Monte Carlo 
simulation using EMCEE (Foreman-Mackey et al. 2013), we 
find fb = 0.58;|{j-34d, a\ = 0.42j^ j 3 , and 012 = -l-lOj^^i, with 

The BYUAVMO R-band detection at 35.25 d reported by Moody et al. 
(2010) is about 0.5 mag brighter than observations from KPNO/SARA 
(Updike et al. 2010a) and VLT/X-shooter Malesani & Palazzi (2010) at a 
similar time, suggesting that the BYU data point may be plagued by a typo. 
We do not include the BYUAVMO /?-band data point at 35.25 d in our analy¬ 
sis. 


a fixed value of y = 5.0. 

We extract XRT spectra at 0.06-0.08 d (mean photon ar¬ 
rival 0.07 d), and at 0.58-1.81 d (mean photon arrival 0.88 d). 
We extrapolate the second SED to 1.5 d and show the com¬ 
plete optical to X-ray SED at 0.07 d and 1.5 d in Figure 8 . 
The spectral index between the UVOT B-band observation at 
0.07 d and the X-rays is /3opt-x = -0.93 ± 0.09, which is con¬ 
sistent with the X-ray spectral index of fix = -1.08 ±0.18 
over this period, suggesting that the optical and X-ray are 
located on the same part of the afterglow SED at 0.07 d. 
The spectral index between the NIR and X-rays at 1.5 d is 
/3nir-x = -0.94 ±0.04, which is consistent with the spectral 
index in the NIR and optical alone, /3NiR-opt = -0.8 ±0.1, 
once again suggesting that the NIR/optical and X-rays are 
located on the same part of the afterglow SED at 1.5 d and 
that intrinsic extinction is negligible. This implies that either 
< t'NiR with /7« 3 or t'ni, Vc < J^iR < vx with 

p Ri 2 at both 0.07 d and 1.5 d. If < i/nir < vx < Vc at 
1.5 d, the decay rate would be a r: -1.5 in the ISM model or 
a R -2.0 in the wind model, whereas in the latter we would 
expect a R -1 in both ISM and wind scenarios. The observed 
UV/NIR/X-ray common decay rate after the re-brightening of 
a R -1 indicates that < J^nir < t'x at 1.5 d and p !v2. 

The 8.46 GFlz light curve drops as f“^ ^ between the ATCA 
observation at 2 d and the VLA observation at 3 d, while the 
SMA 345 GHz light curve is flat over this period (Eigure 9). 
At R2d, the radio spectral index is r 0.7 from 4.9 GHz to 
8.46 GHz, and r 0.3 from 8.46 GHz to 345 GHz. The rapid 
decline in the light curve at a single waveband can not be ex¬ 
plained in the standard synchrotron model by forward shock 
emission. We consider an alternative scenario, in which the 
radio to millimeter emission at 2 d is dominated by a reverse 
shock. In that case, the peak frequency of the reverse shock 
must pass through 8.46 GHz before 2 d with a peak flux den¬ 
sity of at least r 1.3 mJy. Propagating the RS spectrum back 
in time assuming the peak frequency and peak flux evolve as 
f“'-^ and f“' respectively (Kobayashi & Sari 2000), we would 
expect a minimum /?-band flux density of r 2Jy at 120 s. 
This is more than five orders of magnitude brighter than Swift 
Mu'fe-band observations at this time. Thus the rapid decline 
and re-brightening of the 8.46 GHz light curve between 2 and 
20 d cannot be explained in the standard synchrotron frame¬ 
work. It is possible that these data suffer from strong ISS, and 
we therefore exclude these data from our analysis. 

6.1.2. Forward shock model at t >0.4 d 

The optical and X-ray light curves after 0.4 d are insensitive 
to the circumburst density profile since < J^ir < >x 

after the re-brightening. The best-sampled radio light curve 
is the composite 8.46 GHz formed from VLA, VLBA, and 
ATCA observations. However, this light curve exhibits sig¬ 
nificant scatter about a smooth power law evolution, possibly 
due to either ISS or inter-calibration issues across the three 
telescopes'^. Thus we are unable to constrain the circumburst 
density profile from the afterglow data. Eor the remainder of 
this section we focus on the ISM model, discussing the possi¬ 
bility of a wind environment in Appendix C. 

Since this event occurred at a relatively low redshift, 
z = 0.6235 (Antonelli et al. 2010), the host galaxy is 
detected in the Sloane Digital Sky Survey (SDSS) as 
SDSS J170527.10H-112742.5. We use the SDSS photome- 

To avoid the radio data from driving the fit, we apply an uncertainty 
floor of 15% to the cm- and mm-band observations. 
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Figure 8. Observed NIR to X-ray SED of the afterglow of GRB 100418A at 
0.07 d (orange circles) and after the re-brightening (1.5 d; black squares), to¬ 
gether with the best-fit ISM model (Section 6.1.2) including energy injection 
(Section 6.1.3) at 0.07 d (orange; dashed) and 1.5 d (black; solid). The data 
points and model at 0.07 d have been multiplied by a factor of 0.1 for clarity. 
The XRT SED at 1.5 d has been extrapolated from 0.88 d using the best-fit 
broken power law model to the XRT light curve (Figure 7 ; the con'ection fac¬ 
tor is Ri 0.5). The optical data have been extrapolated using the best fit to 
the SwiftfUNOI light curve (the corrections are small at < 5%). The cooling 
break (visible in the orange, dashed curve at Ri 3 X lO*'* Hz) is already below 
the optical bands at 0.07; the NIR to X-ray frequencies are therefore on the 
same part of the afterglow SED at 1.5 d. 

try (m = 24.55 ± 1.14, g = 22.92 ± 0.17, r = 22.45 ±0.17,/ = 
21.95±0.17,z = 22.54±0.83) as fixed a-priori measurements 
of the host galaxy flux in the griz bands. We compute the host 
flux density at B- andV-band usingB = g-l-0.33(g-r)-l-0.2 and 
V =g-0.58(g-r)-0.01 (Jester et al. 2005) and hold these val¬ 
ues fixed for our MCMC analysis. Since there is no a-priori 
measurement of the host flux density in the UVOT/Mi/fe- 
band and much of the UVOT data was taken in Mi/fe-band, 
we keep the host flux density in the White-hwA as a free pa¬ 
rameter, and integrate over the GRB spectral energy distribu¬ 
tion using the W/i/fe-band Alter function'®. 

Using our MCMC analysis as described in Section 4, we fit 
the data after 0.4 d, and our highest-likelihood model (Figure 
9) has the parameters p Ri 2.1, Ce ~ 0.12, eg w 1.1 x 10“^, yiq « 
1.4cm“^, Ek.Iso ~ 3.4 x 10®^erg, fjet ~ 17d, Ay ^ 0.1 mag, 
and Tv, host,White ~ 1.6/iJy, with a Compton y-parameter of 
« 2.9, indicating IC cooling is important. The blastwave 
Lorentz factor is T = 4.9(f/ld)“^/^ and the jet opening angle 
is ('jet ~ 20°. The beaming-corrected kinetic energy is Ek ~ 
2.1 X 10®* erg, while the beaming corrected 7 -ray energy is 
E^^6x 10"*® erg (1-10"* keV; rest frame), corresponding to an 
extremely low radiative efficiency of prad = E^/{Ej^ + Ej) « 
3%. These results are summarized in Table 5. We present his¬ 
tograms of the marginalized posterior density for all param¬ 
eters in Figure 10 and contour plots of correlations between 
the physical parameters in Figure 11 . 

In concordance with the basic analysis outlined above, we 
And that the break frequencies at 1 d for this model are located 
at iza « 9.0GHz, Ri 3.4 x lO" Hz, and ^ 2.3 x 10'^Hz, 
while the peak flux density is about 12mJy at zzm- The spec¬ 
trum transitions from fast to slow cooling at about 230 s af¬ 
ter the burst, and drops below zza at about 25 d. The 

http : //www.swift.ac.uk/analysis/uvot/filters . 

php 


high value of zza is required to suppress the flux density at 
8.46 GHz and 5 GHz (VLA), relative to the mm-band (PdBI 
and SMA). The ordering of the break frequencies relative to 
the observing bands, zz^, zz^, < zz^ir at 1 d, ensures that the op¬ 
tical and X-ray are on the same part of the synchrotron spec¬ 
trum at this time. The jet break at r:! 17 d is largely driven by 
the millimeter and radio data. The blastwave becomes non- 
relativistic at /nr ~ 40 d, resulting in a subsequent slow de¬ 
cline (Tv fO. 6 + 3 ^ ^ t with f3 = (l-p)/2 Ri-0.5; Frail et al. 
2000 ) in the radio light curves after the peak. 

Moin et al. (2013) use a VLB A limit on the size of 
the afterglow at 65 d to estimate the average ejecta 
Lorentz factor. However, we note that the appar¬ 
ent physical size of the afterglow is given by R± = 

[2'2-3*^(4-k)®-V(5-k)®-*] ~ 2.3rcfz, where 7 , is 

the ejecta Lorentz factor of the ejecta at time, E = fobs/(l +z), 
and k = 0 for the ISM model (Granot & Sari 2002); therefore 
the VLBI measurement actually represents an upper limit to 
the ejecta Lorentz factor at the time of the measurement, such 
that r(fobs = 65d) < 5.8. Since the blastwave becomes non- 
relativistic at around 40 days, our model does not violate the 
VLBA limit on the angular size of the afterglow. 


6.1.3. Energy injection model 

Taking the forward shock model described in Section 6 .1.2 
as a starting point, we And that the X-ray and UV/optical data 
before the re-brightening can be explained by two successive 
periods of energy injection. 


t > to = 0.45 d 
fi=0.05d<f <fo 

f 2 = L7 X 10“®d <t <ti 
t < t2 

(4) 

In this model, the energy increases by a factor of « 11 
from TiK.iso.i ~ 1.0 X 10®" erg at 1.7 x 10“® d to iiK.iso ~ 1 . 1 x 
10®* erg at 0.05 d followed by another increase by a factor of 
Ri 30 to TiK.iso.f 7a 3.4 X 10®^ erg at 0.45 d. The overall increase 
in energy is a factor of Ri 340, corresponding to an injected 
energy fraction of Ri 99.7% over this period. In comparison, 
± 7 ,iso = 9 . 9^3 4 X 10®** (Marshall et al. 201 1). The blastwave 
Lorentz factor decreases from L Ri 38 at the start of energy in¬ 
jection at 1.7 X 10“® to r Ri 11 at 0.05 d, and then to L Ri 7 at 
the end of energy injection at 0.45 d. The values of m derived 
above correspond to an ejecta distribution with £(> T) oc r*“* 
with i Ri 9.6 for 7 < r < 11 and i Ri 3 for 11 < r < 38. 



6.2. GRB 100901A 

6.2.1. GRB properties and basic considerations 

GRB 100901A was detected and localized by the Swift BAT 
on 2010 September 01 at 13:34; 10UT (Immleret al. 2010). 
The burst duration is 7go = (439 ± 33) s, with a fluence of 
± 3 , = (2.1 ±0.3) X 10“®ergcm“2 (15-150keV observer frame; 
Immleret al. 2010). The afterglow was detected in the X- 
rays and optical bands by XRT (Page & Immler 2010), UVOT 
(Pritchard & Immler 2010), and by multiple ground-based 
observatories (e.g. De Cia et al. 2010; Kurodaetal. 2010c; 
Hentunen et al. 2010; Updike et al. 2010b; Sahu et al. 2010b; 
Kopac et al. 2010; Gorbovskoy et al. 2010; Pandey & Zheng 
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Figure 9. X-ray, UV (top left), optical (top right), and radio (bottom left) light curves of GRB 100418A, with the full afterglow model (solid lines), including 
energy injection before 0.04 d. The X-ray data before 0.008 d is likely dominated by high-latitude prompt emission and we do not include these data in our 
analysis; the best fit power law to the X-ray data before 0.008 d added to the blastwave model is shown in the upper left panel (black, dashed). The dashed 
envelopes around the radio light curves indicate the expected effect of scintillation at the Icr level. The data indicated by open symbols are not used to determine 
the parameters of the forward shock (the MCMC analysis). Bottom right: blastwave Lorentz factor (green, dashed; upper sub-panel) and isotropic equivalent 
kinetic energy (red, solid; upper sub-panel) as a function of time, together with the energy distribution across ejecta Lorentz factors (black, solid; lower sub-panel) 
as determined from fitting the X-rayAJV/optical re-brightening at 0.7 d. 






Figure 11. Icr (red), 2cr (green), and 3cr (black) contours for correlations 
between the physical parameters, £’K,iso 5 ^e, and eg for GRB 100418A, in 
the ISM model from Monte Carlo simulations. We have restricted Ce < 
and CB < V3. See the on line version of this Figure for additional plots of 
correlations between these parameters and p, /jet, ^jet^ and 7^,host,White- 


Figure 10. Posterior probability density functions for the physical parame¬ 
ters for GRB 100418A in a constant density environment from MCMC sim¬ 
ulations. 

2010), as well as in the radio by the WSRT and the VLA 
(van der Horst et al. 2010b; Chandra & Frail 2010a). Spec¬ 


troscopic observations with Magellan yielded a redshift of 
z= 1.408 (Chornock et al. 2010). Using the BAT fluence, 
Gorbovskoy et al. (2012) determined the isotropic equivalent 
7 -ray energy of this burst to be F-^ iso = 6.3 x 10^^ erg (1- 
lO^keV; rest frame). However, the 7 -ray spectrum does 
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not exhibit a turnover within the BAT energy range, indi¬ 
cating that the peak of the 7-ray spectrum, Epeak, is outside 
the BAT energy range. Margutti et al. (2013) obtain a cor¬ 
relation between EnAT.obs (15-150keV) and E-^.iso: = 


10 - 3 ’ 


l.OSibO.Ol 


with a scatter of cr = 0.24 (68%), 


which we use to perform a K-correction and obtain E-yjso = 
(8 ± 1) X lO^’ erg. Note that the error bars do not include the 


uncertainty in the correlation itself. 


Table 7 

SwiftAJVOT Observations of GRB 100901A 


?-?0 

(days) 

Filter 

Frequency 

(Hz) 

Flux density^ 
(AiJy) 

Uncertainty 

(AiJy) 

Detection? 

(1 = Yes) 

0.00258 

White 

8.64e+14 

27.4 

5.65 

1 

0.00498 

u 

8.56e+14 

146 

11.4 

1 

0.00661 

b 

6.92e+14 

230 

75.6 

0 

0.00689 

White 

8.64e+14 

58.6 

18.9 

1 

0.00717 

uvw2 

1.48e+15 

67.9 

12.1 

0 

0.0454 

V 

5.55e+14 

155 

27 

1 

0.0478 

uvm2 

1.34e+15 

18.3 

4.42 

1 

0.0502 

uvwl 

1.16e+15 

37.1 

5.4 

1 

0.0597 

uvw2 

1.48e+15 

7.4 

2.36 

0 

0.0609 

u 

8.56e+14 

83.4 

6.21 

1 

0.062 

V 

5.55e+14 

161 

27.1 

1 

0.0632 

b 

6.92e+14 

134 

12.9 

1 

0.0644 

uvm2 

1.34e+15 

20.7 

4.63 

1 

0.0652 

White 

8.64e+14 

62.2 

9.12 

1 

0.0668 

uvwl 

1.16e+15 

17.2 

4.39 

1 

0.117 

uvw2 

1.48e+15 

3.45 

1.02 

1 

0.127 

V 

5.55e+14 

97.7 

12 

1 

0.137 

uvm2 

1.34e+15 

10.1 

2 

1 

0.183 

uvwl 

1.16e+15 

42.7 

3.03 

1 

0.194 

u 

8.56e+14 

129 

5.18 

1 

0.204 

b 

6.92e+14 

192 

10.2 

1 

0.25 

uvm2 

1.34e+15 

21.3 

2.22 

1 

0.261 

uvwl 

1.16e+15 

58.3 

3.6 

1 

0.271 

u 

8.56e+14 

160 

6.96 

1 

0.318 

uvw2 

1.48e+15 

11.4 

1.37 

1 

0.329 

V 

5.55e+14 

291 

15.5 

1 

0.338 

uvm2 

1.34e+15 

26.6 

2.92 

1 

0.384 

uvwl 

1.16e+15 

40.4 

8.67 

1 

0.534 

uvwl 

1.16e+15 

36.8 

2.88 

1 

0.541 

u 

8.56e+14 

107 

10.1 

1 

0.599 

uvw2 

1.48e+15 

6.15 

1.19 

1 

0.607 

V 

5.55e+14 

178 

28.3 

1 

0.647 

b 

6.92e+14 

130 

43.4 

0 

0.717 

uvm2 

1.34e+15 

8.27 

2.26 

1 

0.736 

uvwl 

1.16e+15 

24.4 

2.5 

1 

0.743 

u 

8.56e+14 

63.4 

13.4 

1 

0.786 

uvw2 

1.48e+15 

3.82 

1.01 

1 

0.796 

V 

5.55e+14 

117 

12.3 

1 

0.806 

uvm2 

1.34e+15 

12.7 

2.14 

1 

0.853 

uvwl 

1.16e+15 

16.8 

2.01 

1 

0.863 

u 

8.56e+14 

54.4 

3.36 

1 

0.873 

b 

6.92e+14 

64.8 

8.72 

1 

0.92 

uvm2 

1.34e+15 

6.86 

1.49 

1 

0.93 

uvwl 

1.16e+15 

19.1 

2.14 

1 

0.94 

u 

8.56e+14 

41.4 

4.63 

1 

0.983 

uvw2 

1.48e+15 

7.57 

2.11 

0 

1.06 

V 

5.55e+14 

74.5 

11.6 

1 

1.13 

uvm2 

1.34e+15 

6.64 

1.18 

1 

1.14 

b 

6.92e+14 

51.5 

8.67 

1 

1.15 

uvw2 

1.48e+15 

2.53 

0.828 

0 

1.16 

uvwl 

1.16e+15 

12.2 

1.35 

1 

1.17 

u 

8.56e+14 

31 

2.53 

1 

1.87 

uvw2 

1.48e+15 

1.57 

0.51 

0 

1.87 

uvm2 

1.34e+15 

2.87 

0.771 

1 

1.88 

uvwl 

1.16e+15 

5.66 

1.13 

1 

1.93 

u 

8.56e+14 

15.5 

1.74 

1 

1.94 

b 

6.92e+14 

26.9 

3.6 

1 

1.94 

White 

8.64e+14 

10.7 

1.69 

1 

1.94 

V 

5.55e+14 

54.1 

8.97 

1 

2.8 

b 

6.92e+14 

17.4 

5.28 

1 

2.8 

White 

8.64e+14 

5.82 

1.45 

1 


Table 7 — Continued 


t-to 

Filter Frequency 

Flux density^ 

Uncertainty Detection? 

(days) 

(Hz) 

(/rJy) 

(/rJy) (l=Yes) 


2.81 

V 

5.55e+14 

35.6 

11.7 

0 

2.86 

u 

8.56e+14 

6.97 

2.32 

0 

2.93 

uvw2 

1.48e+15 

2.38 

0.694 

0 

2.93 

uvm2 

1.34e+15 

3.5 

0.917 

0 

2.94 

uvwl 

1.16e+15 

4.22 

1.29 

0 

4.2 

ii 

8.56e+14 

8.36 

2.64 

0 

4.2 

uvw2 

1.48e+15 

3.99 

0.955 

0 

4.27 

b 

6.92e+14 

24.6 

7.94 

0 

4.27 

White 

8.64e+14 

5.56 

1.89 

0 

4.28 

V 

5.55e+14 

53.8 

17.3 

0 

4.34 

uvm2 

1.34e+15 

7.14 

1.64 

0 

4.35 

uvwl 

1.16e+15 

8.36 

2.24 

0 

4.8 

u 

8.56e+14 

8.03 

2.61 

0 

4.8 

uvw2 

1.48e+15 

4.13 

1.1 

0 

4.8 

b 

6.92e+14 

17.4 

5.63 

0 

4.8 

uvm2 

1.34e+15 

6.04 

1.48 

0 

4.81 

White 

8.64e+14 

3.98 

1.36 

0 

4.81 

uvwl 

1.16e+15 

7.16 

1.97 

0 

4.81 

V 

5.55e+14 

53.2 

16.9 

0 

5.67 

u 

8.56e+14 

33.1 

9.34 

0 

5.8 

uvw2 

1.48e+15 

10.1 

2.13 

0 

7.24 

uvw2 

1.48e+15 

1.39 

0.394 

0 

7.25 

uvm2 

1.34e+15 

1.83 

0.542 

0 

7.25 

uvwl 

1.16e+15 

2.41 

0.758 

0 

9.46 

V 

5.55e+14 

19.8 

6.4 

0 

10.2 

11 

8.56e+14 

3.96 

1.28 

0 

10.2 

b 

6.92e+14 

8.75 

2.87 

0 

10.2 

White 

8.64e+14 

2.05 

0.678 

0 

11 

11 

8.56e+14 

4.21 

1.36 

0 

11 

b 

6.92e+14 

9.31 

3.05 

0 

11 

White 

8.64e+14 

2.26 

0.743 

0 

11 

V 

5.55e+14 

23 

7.48 

0 

11.7 

u 

8.56e+14 

4.41 

1.4 

0 

11.7 

b 

6.92e+14 

9.97 

3.23 

0 

11.7 

White 

8.64e+14 

2.55 

0.859 

0 

11.8 

V 

5.55e+14 

30.4 

9.91 

0 

12.6 

11 

8.56e+14 

4.11 

1.33 

0 

12.6 

b 

6.92e+14 

8.96 

2.94 

0 

12.6 

White 

8.64e+14 

2.12 

0.703 

0 

12.6 

V 

5.55e+14 

23.1 

7.59 

0 

13.5 

b 

6.92e+14 

20.4 

6.56 

0 

13.5 

White 

8.64e+14 

5.14 

1.68 

0 

13.5 

V 

5.55e+14 

62 

19.5 

0 

13.6 

u 

8.56e+14 

8.12 

2.58 

0 

14.6 

11 

8.56e+14 

4.3 

1.42 

0 

14.6 

b 

6.92e+14 

9.41 

3.08 

0 

14.7 

White 

8.64e+14 

2.19 

0.721 

0 

14.7 

V 

5.55e+14 

23.5 

7.62 

0 

15.7 

11 

8.56e+14 

4.1 

1.31 

0 

15.7 

b 

6.92e+14 

8.93 

2.95 

0 

15.7 

White 

8.64e+14 

2.08 

0.693 

0 

15.7 

V 

5.55e+14 

22.7 

7.38 

0 

16.7 

u 

8.56e+14 

4.13 

1.32 

0 

16.8 

b 

6.92e+14 

9.21 

2.99 

0 

16.8 

White 

8.64e+14 

2.41 

0.794 

0 

16.8 

V 

5.55e+14 

24.8 

8.08 

0 

17.7 

u 

8.56e+14 

4.29 

1.38 

0 

17.7 

b 

6.92e+14 

9.39 

3.09 

0 

17.7 

White 

8.64e+14 

2.15 

0.707 

0 

17.7 

V 

5.55e+14 

21.8 

7.11 

0 

18.6 

u 

8.56e+14 

4.27 

1.37 

0 

18.6 

b 

6.92e+14 

9.32 

3.04 

0 

18.6 

White 

8.64e+14 

2.13 

0.701 

0 

18.6 

V 

5.55e+14 

21.3 

6.85 

0 


Note. — tjn cases of non-detections, we report the formal flux measure¬ 
ment at the position of the afterglow. 


XRT began observing the GRB during the 7-ray emission 
interval. The X-ray spectra from 1 ks to the end of XRT ob¬ 
servations are well fit by an absorbed single power law model 
with NH,Gai = 7.1 X 10’°cm"’ (Kalberlaet al. 2005), intrinsic 
hydrogen column An.mt = (3.1 ±0.7) x lO^'cm"’, and photon 
index, Tx = 2.15 ±0.06, with no evidence for spectral evo- 
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lution during this period. We convert the count rate light 
curve published by Margutti et al. (2013) to a flux density 
light curve at 1 keV using their time-resolved spectra. We also 


analyze all UVOT photometry for this burst and report our re¬ 
sults in Table 7, and a compilation of all photometry listed 
in GCN circulars or published elsewhere for this event in 
Table 8. 


Table 8 

Optical Observations of GRB 100901A 


t-to 

(days) 

Observatory 

Telescope / 
Instrument 

Filter 

Frequency 

(Hz) 

Flux density^ 

(z^Jy) 

Uncertainty 

(ply) 

Detection? 
(1 = Yes) 

Reference 

0.00131 

MASTER 


CR 

4.4e+14 

82.1 

436 

0 

Gorbovskoy et al. (2012) 

0.00178 

MASTER 


CR 

4.4e+14 

156 

170 

0 

Gorbovskoy et al. (2012) 

0.00246 

MASTER 


CR 

4.4e+14 

131 

143 

0 

Gorbovskoy et al. (2012) 

0.00396 

MASTER 


CR 

4.4e+14 

119 

107 

0 

Gorbovskoy et al. (2012) 

0.00494 

MASTER 


CR 

4.4e+14 

304 

78.6 

1 

Gorbovskoy et al. (2012) 

0.00615 

MASTER 


CR 

4.4e+!4 

125 

73.4 

0 

Gorbovskoy et al. (2012) 

0.00758 

MASTER 


CR 

4.4e+14 

88.4 

80 

0 

Gorbovskoy et al. (2012) 

0.0093 

MASTER 


CR 

4.4e+14 

109 

80.6 

0 

Gorbovskoy et al. (2012) 

0.0114 

MASTER 


CR 

4.4e+14 

204 

65 

1 

Gorbovskoy et al. (2012) 

0.012 

Okayama 

MITSuME 

Ic 

3.93e+14 

460 

93.1 

1 

Kuroda et al. (2010c) 

0.012 

Okayama 

MITSuME 

Rc 

4.56e+14 

336 

32.4 

1 

Kuroda et al. (2010c) 

0.012 

Okayama 

MITSuME 

g’ 

6.29e+14 

132 

77.1 

1 

Kuroda et al. (2010c) 

0.0136 

MASTER 


CR 

4.4e+14 

195 

62.1 

1 

Gorbovskoy et al. (2012) 

0.0159 

MASTER 


V 

5.55e+14 

292 

55.8 

1 

Gorbovskoy et al. (2012) 

0.0182 

MASTER 


CR 

4.4e+14 

232 

74 

1 

Gorbovskoy et al. (2012) 

0.0227 

MASTER 


V 

5.55e+14 

292 

55.8 

1 

Gorbovskoy et al. (2012) 

0.0251 

MASTER 


V 

5.55e+14 

351 

59.4 

1 

Gorbovskoy et al. (2012) 

0.0275 

MASTER 


CR 

4.4e+14 

228 

22 

1 

Gorbovskoy et al. (2012) 

0.0299 

MASTER 


CR 

4.4e+14 

224 

12.7 

1 

Gorbovskoy et al. (2012) 

0.0323 

MASTER 


CR 

4.4e+14 

202 

11.5 

1 

Gorbovskoy et al. (2012) 

0.0334 

Okayama 

MITSuME 

Ic 

3.93e+14 

265 

25.6 

1 

Kuroda et al. (2010c) 

0.0334 

Okayama 

MITSuME 

Rc 

4.56e+14 

212 

20.5 

1 

Kuroda et al. (2010c) 

0.0334 

Okayama 

MITSuME 

g’ 

6.29e+14 

209 

20.2 

1 

Kuroda et al. (2010c) 

0.0382 

MASTER 


CR 

4.4e+14 

186 

12.4 

1 

Gorbovskoy et al. (2012) 

0.0406 

MASTER 


CR 

4.4e+14 

195 

13 

1 

Gorbovskoy et al. (2012) 

0.043 

MASTER 


CR 

4.4e+14 

192 

12.8 

1 

Gorbovskoy et al. (2012) 

0.0454 

MASTER 


CR 

4.4e+14 

214 

12.2 

1 

Gorbovskoy et al. (2012) 

0.0478 

MASTER 


CR 

4.4e+14 

185 

12.3 

1 

Gorbovskoy et al. (2012) 

0.0502 

MASTER 


CR 

4.4e+14 

212 

14.1 

1 

Gorbovskoy et al. (2012) 

0.0526 

MASTER 


CR 

4.4e+14 

155 

13.4 

1 

Gorbovskoy et al. (2012) 

0.0548 

Okayama 

MITSuME 

Ic 

3.93e+14 

220 

21.3 

1 

Kuroda et al. (2010c) 

0.0548 

Okayama 

MITSuME 

Rc 

4.56e+14 

122 

11.8 

1 

Kuroda et al. (2010c) 

0.0548 

Okayama 

MITSuME 

g’ 

6.29e+14 

120 

11.6 

1 

Kuroda et al. (2010c) 

0.055 

MASTER 


CR 

4.4e+!4 

117 

13.6 

1 

Gorbovskoy et al. (2012) 

0.0574 

MASTER 


CR 

4.4e+14 

109 

15 

1 

Gorbovskoy et al. (2012) 

0.0599 

MASTER 


CR 

4.4e+14 

183 

14 

1 

Gorbovskoy et al. (2012) 

0.0623 

MASTER 


CR 

4.4e+14 

185 

14.1 

1 

Gorbovskoy et al. (2012) 

0.0647 

MASTER 


CR 

4.4e+14 

186 

14.3 

1 

Gorbovskoy et al. (2012) 

0.067 

MASTER 


CR 

4.4e+14 

193 

14.8 

1 

Gorbovskoy et al. (2012) 

0.0695 

MASTER 


CR 

4.4e+14 

164 

14.2 

1 

Gorbovskoy et al. (2012) 

0.0721 

MASTER 


CR 

4.4e+14 

206 

15.8 

1 

Gorbovskoy et al. (2012) 

0.0745 

MASTER 


CR 

4.4e+14 

167 

16.1 

1 

Gorbovskoy et al. (2012) 

0.0759 

Okayama 

MITSuME 

Ic 

3.93e+14 

152 

30.8 

1 

Kuroda et al. (2010c) 

0.0759 

Okayama 

MITSuME 

Rc 

4.56e+14 

134 

12.9 

1 

Kuroda et al. (2010c) 

0.0759 

Okayama 

MITSuME 

g’ 

6.29e+14 

120 

11.6 

1 

Kuroda et al. (2010c) 

0.0804 

MASTER 


CR 

4.4e+14 

131 

14 

1 

Gorbovskoy et al. (2012) 

0.0828 

MASTER 


CR 

4.4e+14 

no 

15.2 

1 

Gorbovskoy et al. (2012) 

0.0852 

MASTER 


CR 

4.4e+14 

101 

16 

1 

Gorbovskoy et al. (2012) 

0.0877 

MASTER 


CR 

4.4e+14 

61.2 

18.7 

1 

Gorbovskoy et al. (2012) 

0.0901 

MASTER 


CR 

4.4e+14 

106 

16.9 

1 

Gorbovskoy et al. (2012) 

0.0926 

MASTER 


CR 

4.4e+14 

141 

16.5 

1 

Gorbovskoy et al. (2012) 

0.095 

MASTER 


CR 

4.4e+14 

156 

18.3 

1 

Gorbovskoy et al. (2012) 

0.0969 

Okayama 

MITSuME 

Ic 

3.93e+14 

152 

30.8 

1 

Kuroda et al. (2010c) 

0.0969 

Okayama 

MITSuME 

Rc 

4.56e+14 

122 

11.8 

1 

Kuroda et al. (2010c) 

0.0969 

Okayama 

MITSuME 

g’ 

6.29e+14 

100 

9.65 

1 

Kuroda et al. (2010c) 

0.118 

Okayama 

MITSuME 

Ic 

3.93e+14 

167 

33.8 

1 

Kuroda et al. (2010c) 

0.118 

Okayama 

MITSuME 

Rc 

4.56e+14 

101 

20.5 

1 

Kuroda et al. (2010c) 

0.118 

Okayama 

MITSuME 

g’ 

6.29e+14 

83.2 

16.8 

1 

Kuroda et al. (2010c) 

0.131 

MASTER 


CR 

4.4e+14 

22 

25.7 

0 

Gorbovskoy et al. (2012) 

0.134 

MASTER 


CR 

4.4e+14 

74.9 

17.7 

1 

Gorbovskoy et al. (2012) 

0.136 

MASTER 


CR 

4.4e+14 

58.9 

18.8 

1 

Gorbovskoy et al. (2012) 

0.139 

MASTER 


CR 

4.4e+14 

54.3 

35 

0 

Gorbovskoy et al. (2012) 

0.139 

Okayama 

MITSuME 

Ic 

3.93e+14 

183 

17.7 

1 

Kuroda et al. (2010c) 

0.139 

Okayama 

MITSuME 

Rc 

4.56e+14 

134 

12.9 

1 

Kuroda et al. (2010c) 


The KPNO/SARA V-band detection (Updike et al. 2010b) is signif¬ 
icantly brighter than expected from interpolating the UVOT V-band light 


curve, implying a difference either in calibration or in the filter response be¬ 
tween the two instruments. We do not include the KPNO V-band data point 
at 0.78 d in our analysis. 
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Table 8 — Continued 


t-tQ 

(days) 

Observatory 

Telescope / 
Instnament 

Filter 

Frequency 

(Hz) 

Flux density^ 

(z^Jy) 

Uncertainty 

(/iJy) 

Detection? 
(1 = Yes) 

Reference 

0.139 

Okayama 

MITSuME 

g' 

6.29e+14 

no 

10.6 

1 

Kuroda et al. (2010c) 

0.141 

MASTER 


CR 

4.4e+14 

108 

31.8 

1 

Gorbovskoy et al. (2012) 

0.144 

MASTER 


CR 

4.4e+14 

152 

30.8 

1 

Gorbovskoy et al. (2012) 

0.146 

MASTER 


CR 

4.4e+14 

106 

32.5 

1 

Gorbovskoy et al. (2012) 

0.149 

MASTER 


CR 

4.4e+14 

185 

29.3 

1 

Gorbovskoy et al. (2012) 

0.151 

MASTER 


CR 

4.4e+14 

135 

22.9 

1 

Gorbovskoy et al. (2012) 

0.154 

MASTER 


CR 

4.4e+14 

193 

18.7 

1 

Gorbovskoy et al. (2012) 

0.16 

Okayama 

MITSuME 

Ic 

3.93e+14 

319 

30.7 

1 

Kuroda et al. (2010c) 

0.16 

Okayama 

MITSuME 

Rc 

4.56e+14 

193 

18.7 

1 

Kuroda et al. (2010c) 

0.16 

Okayama 

MITSuME 

g’ 

6.29e+14 

120 

11.6 

1 

Kuroda et al. (2010c) 

0.172 

Hanle 

HCT 

R 

4.56e+14 

267 

12.6 

1 

Sahu et al. (2010b) 

0.172 

MASTER 


R 

4.56e+14 

180 

92.3 

0 

Gorbovskoy et al. (2012) 

0.178 

Hanle 

HCT 

I 

3.93e+14 

356 

16.8 

1 

Sahu et al. (2010b) 

0.181 

Okayama 

MITSuME 

Ic 

3.93e+14 

349 

33.7 

1 

Kuroda et al. (2010c) 

0.181 

Okayama 

MITSuME 

Rc 

4.56e+14 

280 

27 

1 

Kuroda et al. (2010c) 

0.181 

Okayama 

MITSuME 

g’ 

6.29e+14 

191 

18.4 

1 

Kuroda et al. (2010c) 

0.182 

Terksol 

Z-600 

R 

4.56e+14 

282 

7.9 

1 

Andreev et al. (2010a) 

0.183 

MASTER 


CR 

4.4e+14 

202 

36.5 

1 

Gorbovskoy et al. (2012) 

0.183 

MASTER 


R 

4.56e+14 

230 

39 

1 

Gorbovskoy et al. (2012) 

0.186 

MASTER 


CR 

4.4e+14 

220 

14.7 

1 

Gorbovskoy et al. (2012) 

0.186 

MASTER 


R 

4.56e+14 

290 

30.9 

1 

Gorbovskoy et al. (2012) 

0.188 

MASTER 


CR 

4.4e+14 

239 

13.6 

1 

Gorbovskoy et al. (2012) 

0.188 

MASTER 


R 

4.56e+14 

277 

26.7 

1 

Gorbovskoy et al. (2012) 

0.19 

MASTER 


R 

4.56e+14 

321 

24.5 

1 

Gorbovskoy et al. (2012) 

0.193 

MASTER 


R 

4.56e+14 

239 

23.1 

1 

Gorbovskoy et al. (2012) 

0.195 

MASTER 


R 

4.56e+14 

358 

23.9 

1 

Gorbovskoy et al. (2012) 

0.198 

MASTER 


R 

4.56e+14 

290 

25.1 

1 

Gorbovskoy et al. (2012) 

0.2 

MASTER 


R 

4.56e+14 

321 

21.4 

1 

Gorbovskoy et al. (2012) 

0.202 

Okayama 

MITSuME 

Ic 

3.93e+14 

383 

37 

1 

Kuroda et al. (2010c) 

0.202 

Okayama 

MITSuME 

Rc 

4.56e+14 

255 

24.6 

1 

Kuroda et al. (2010c) 

0.202 

Okayama 

MITSuME 

g’ 

6.29e+14 

229 

22.1 

1 

Kuroda et al. (2010c) 

0.202 

MASTER 


R 

4.56e+14 

339 

22.6 

1 

Gorbovskoy et al. (2012) 

0.205 

MASTER 


R 

4.56e+14 

269 

20.6 

1 

Gorbovskoy et al. (2012) 

0.207 

MASTER 


R 

4.56e+14 

287 

19.1 

1 

Gorbovskoy et al. (2012) 

0.21 

MASTER 


R 

4.56e+14 

324 

21.6 

1 

Gorbovskoy et al. (2012) 

0.212 

MASTER 


R 

4.56e+14 

304 

20.2 

1 

Gorbovskoy et al. (2012) 

0.214 

MASTER 


R 

4.56e+14 

318 

18.1 

1 

Gorbovskoy et al. (2012) 

0.217 

MASTER 


R 

4.56e+14 

330 

18.7 

1 

Gorbovskoy et al. (2012) 

0.219 

MASTER 


R 

4.56e+14 

301 

17.1 

1 

Gorbovskoy et al. (2012) 

0.222 

MASTER 


R 

4.56e+14 

309 

17.6 

1 

Gorbovskoy et al. (2012) 

0.223 

MASTER 


R 

4.56e+14 

304 

38.6 

1 

Gorbovskoy et al. (2012) 

0.224 

Okayama 

MITSuME 

Ic 

3.93e+14 

383 

37 

1 

Kuroda et al. (2010c) 

0.224 

Okayama 

MITSuME 

Rc 

4.56e+14 

280 

27 

1 

Kuroda et al. (2010c) 

0.224 

Okayama 

MITSuME 

g’ 

6.29e+14 

191 

18.4 

1 

Kuroda et al. (2010c) 

0.224 

MASTER 


R 

4.56e+14 

321 

18.2 

1 

Gorbovskoy et al. (2012) 

0.226 

MASTER 


I 

3.93e+14 

383 

33.1 

1 

Gorbovskoy et al. (2012) 

0.226 

MASTER 


R 

4.56e+14 

355 

20.2 

1 

Gorbovskoy et al. (2012) 

0.227 

MASTER 


R 

4.56e+14 

324 

18.4 

1 

Gorbovskoy et al. (2012) 

0.228 

MASTER 


I 

3.93e+14 

397 

34.3 

1 

Gorbovskoy et al. (2012) 

0.228 

MASTER 


R 

4.56e+14 

293 

19.5 

1 

Gorbovskoy et al. (2012) 

0.229 

MASTER 


R 

4.56e+14 

400 

18.9 

1 

Gorbovskoy et al. (2012) 

0.23 

MASTER 


I 

3.93e+14 

359 

34.6 

1 

Gorbovskoy et al. (2012) 

0.23 

MASTER 


R 

4.56e+14 

333 

18.9 

1 

Gorbovskoy et al. (2012) 

0.231 

MASTER 


R 

4.56e+14 

327 

18.6 

1 

Gorbovskoy et al. (2012) 

0.233 

MASTER 


I 

3.93e+14 

432 

33 

1 

Gorbovskoy et al. (2012) 

0.233 

MASTER 


R 

4.56e+14 

352 

20 

1 

Gorbovskoy et al. (2012) 

0.233 

MASTER 


V 

5.55e+14 

410 

43.7 

1 

Gorbovskoy et al. (2012) 

0.234 

MASTER 


R 

4.56e+14 

321 

18.2 

1 

Gorbovskoy et al. (2012) 

0.235 

MASTER 


I 

3.93e+14 

369 

31.9 

1 

Gorbovskoy et al. (2012) 

0.235 

MASTER 


R 

4.56e+14 

324 

18.4 

1 

Gorbovskoy et al. (2012) 

0.235 

MASTER 


V 

5.55e+14 

371 

43.3 

1 

Gorbovskoy et al. (2012) 

0.236 

MASTER 


R 

4.56e+14 

365 

17.2 

1 

Gorbovskoy et al. (2012) 

0.237 

MASTER 


I 

3.93e+14 

369 

31.9 

1 

Gorbovskoy et al. (2012) 

0.237 

MASTER 


R 

4.56e+14 

372 

17.5 

1 

Gorbovskoy et al. (2012) 

0.239 

MASTER 


R 

4.56e+14 

382 

18 

1 

Gorbovskoy et al. (2012) 

0.24 

MASTER 


I 

3.93e+14 

366 

31.6 

1 

Gorbovskoy et al. (2012) 

0.24 

MASTER 


R 

4.56e+14 

330 

18.7 

1 

Gorbovskoy et al. (2012) 

0.241 

MASTER 


R 

4.56e+14 

330 

18.7 

1 

Gorbovskoy et al. (2012) 

0.242 

MASTER 


I 

3.93e+14 

448 

29.8 

1 

Gorbovskoy et al. (2012) 

0.242 

MASTER 


R 

4.56e+14 

330 

18.7 

1 

Gorbovskoy et al. (2012) 

0.243 

MASTER 


R 

4.56e+14 

312 

20.8 

1 

Gorbovskoy et al. (2012) 

0.245 

MASTER 


I 

3.93e+14 

456 

30.4 

1 

Gorbovskoy et al. (2012) 

0.245 

MASTER 


R 

4.56e+14 

352 

20 

1 

Gorbovskoy et al. (2012) 

0.246 

MASTER 


R 

4.56e+14 

342 

19.4 

1 

Gorbovskoy et al. (2012) 

0.247 

MASTER 


I 

3.93e+14 

412 

27.5 

1 

Gorbovskoy et al. (2012) 

0.247 

MASTER 


R 

4.56e+14 

355 

16.7 

1 

Gorbovskoy et al. (2012) 
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Table 8 — Continued 


t-tQ 

(days) 

Observatory 

Telescope / 
Instrument 

Filter 

Frequency 

(Hz) 

Flux density^ 

(z^Jy) 

Uncertainty 

(/iJy) 

Detection? 
(1 = Yes) 

Reference 

0.248 

MASTER 


R 

4.56e-i-14 

336 

19.1 

1 

Gorbovskoy et al. (2012) 

0.249 

MASTER 


I 

3.93e-i-14 

460 

30.7 

1 

Gorbovskoy et al. (2012) 

0.249 

MASTER 


R 

4.56e-i-14 

342 

16.1 

1 

Gorbovskoy et al. (2012) 

0.251 

MASTER 


R 

4.56e-i-14 

318 

21.2 

1 

Gorbovskoy et al. (2012) 

0.252 

MASTER 


I 

3.93e-i-14 

412 

27.5 

1 

Gorbovskoy et al. (2012) 

0.252 

MASTER 


R 

4.56e-i-14 

352 

16.6 

1 

Gorbovskoy et al. (2012) 

0.253 

MASTER 


R 

4.56e-i-14 

309 

20.6 

1 

Gorbovskoy et al. (2012) 

0.254 

MASTER 


I 

3.93e-i-14 

529 

30 

1 

Gorbovskoy et al. (2012) 

0.254 

MASTER 


R 

4.56e-i-14 

330 

15.5 

1 

Gorbovskoy et al. (2012) 

0.256 

MASTER 


R 

4.56e-i-14 

358 

20.4 

1 

Gorbovskoy et al. (2012) 

0.257 

MASTER 


I 

3.93e-i-14 

482 

27.4 

1 

Gorbovskoy et al. (2012) 

0.257 

MASTER 


R 

4.56e-i-14 

372 

17.5 

1 

Gorbovskoy et al. (2012) 

0.257 

Hanle 

HOT 

R 

4.56e-i-14 

355 

16.7 

1 

Sahu et al. (2010b) 

0.258 

MASTER 


R 

4.56e-i-14 

372 

21.1 

1 

Gorbovskoy et al. (2012) 

0.259 

MASTER 


I 

3.93e-i-14 

444 

25.2 

1 

Gorbovskoy et al. (2012) 

0.259 

MASTER 


R 

4.56e-i-14 

375 

14.1 

1 

Gorbovskoy et al. (2012) 

0.261 

MASTER 


R 

4.56e-i-14 

431 

20.3 

1 

Gorbovskoy et al. (2012) 

0.261 

MASTER 


I 

3.93e-i-14 

505 

23.8 

1 

Gorbovskoy et al. (2012) 

0.262 

Hanle 

HOT 

I 

3.93e-i-14 

465 

17.4 

1 

Sahu et al. (2010b) 

0.264 

MASTER 


I 

3.93e-i-14 

519 

24.5 

1 

Gorbovskoy et al. (2012) 

0.266 

MASTER 


1 

3.93e-i-14 

440 

25 

1 

Gorbovskoy et al. (2012) 

0.269 

MASTER 


I 

3.93e-i-14 

496 

28.2 

1 

Gorbovskoy et al. (2012) 

0.271 

MASTER 


I 

3.93e-i-14 

491 

27.9 

1 

Gorbovskoy et al. (2012) 

0.273 

MASTER 


I 

3.93e-i-14 

482 

27.4 

1 

Gorbovskoy et al. (2012) 

0.274 

Terksol 

Z-600 

R 

4.56e-i-14 

349 

6.48 

1 

Andreev et al. (2010b) 

0.276 

MASTER 


I 

3.93e-i-14 

452 

25.7 

1 

Gorbovskoy et al. (2012) 

0.278 

MASTER 


I 

3.93e-i-14 

469 

26.6 

1 

Gorbovskoy et al. (2012) 

0.281 

MASTER 


I 

3.93e-i-14 

405 

27 

1 

Gorbovskoy et al. (2012) 

0.283 

MASTER 


I 

3.93e-i-14 

448 

25.4 

1 

Gorbovskoy et al. (2012) 

0.286 

MASTER 


I 

3.93e-i-14 

534 

30.3 

1 

Gorbovskoy et al. (2012) 

0.288 

MASTER 


I 

3.93e-i-14 

448 

29.8 

1 

Gorbovskoy et al. (2012) 

0.29 

MASTER 


I 

3.93e-i-14 

482 

27.4 

1 

Gorbovskoy et al. (2012) 

0.293 

MASTER 


I 

3.93e-i-14 

514 

29.2 

1 

Gorbovskoy et al. (2012) 

0.295 

MASTER 


I 

3.93e-i-14 

376 

28.8 

1 

Gorbovskoy et al. (2012) 

0.296 

MASTER 


V 

5.55e-Hl4 

303 

54.6 

1 

Gorbovskoy et al. (2012) 

0.297 

MASTER 


I 

3.93e-i-14 

432 

28.8 

1 

Gorbovskoy et al. (2012) 

0.299 

MASTER 


V 

5.55e-Hl4 

341 

50.5 

1 

Gorbovskoy et al. (2012) 

0.3 

MASTER 


I 

3.93e-i-14 

465 

26.4 

1 

Gorbovskoy et al. (2012) 

0.302 

Golovec 

Vega 

R 

4.56e-i-14 

306 

29.6 

1 

Kopac et al. (2010) 

0.302 

MASTER 


I 

3.93e-i-14 

416 

31.8 

1 

Gorbovskoy et al. (2012) 

0.305 

MASTER 


I 

3.93e-i-14 

452 

34.6 

1 

Gorbovskoy et al. (2012) 

0.307 

MASTER 


I 

3.93e-i-14 

369 

35.6 

1 

Gorbovskoy et al. (2012) 

0.317 

THO 

C-14 

CR 

4.4e-Hl4 

333 

32.1 

1 

Hentunen et al. (2010) 

0.32 

THO 

C-14 

CR 

4.4e-Hl4 

304 

61.4 

1 

Hentunen et al. (2010) 

0.353 

Terksol 

Meade-35 

V 

5.55e-Hl4 

338 

15.9 

1 

Andreev et al. (2010e) 

0.355 

Hanle 

HOT 

R 

4.56e-i-14 

339 

12.7 

1 

Sahu et al. (2010b) 

0.37 

Terksol 

Zeiss-2000 

b 

6.92e-i-14 

232 

4.31 

1 

Andreev et al. (2010e) 

0.416 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

304 

96.6 

1 

Sposetti & Immler (2010) 

0.427 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

285 

90.6 

1 

Sposetti & Immler (2010) 

0.439 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

287 

91.4 

1 

Sposetti & Immler (2010) 

0.442 

THO 

C-14 

CR 

4.4e-Hl4 

277 

26.7 

1 

Hentunen et al. (2010) 

0.449 

Terksol 

Zeiss-2000 

b 

6.92e-i-14 

169 

3.15 

1 

Andreev et al. (2010e) 

0.45 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

277 

88.1 

1 

Sposetti & Immler (2010) 

0.458 

Terksol 

Zeiss-2000 

V 

5.55e-Hl4 

266 

4.94 

1 

Andreev et al. (2010e) 

0.461 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

246 

78.2 

1 

Sposetti & Immler (2010) 

0.464 

Terksol 

Zeiss-2000 

r’ 

4.56e-Hl4 

263 

4.89 

1 

Andreev et al. (2010e) 

0.464 

THO 

C-14 

Rc 

4.56e-i-14 

253 

24.4 

I 

Hentunen et al. (2010) 

0.468 

THO 

C-14 

b 

6.92e-i-14 

234 

47.3 

1 

Hentunen et al. (2010) 

0.471 

Terksol 

Zeiss-2000 

g’ 

6.29e-i-14 

163 

3.03 

1 

Andreev et al. (2010e) 

0.472 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

243 

77.5 

1 

Sposetti & Immler (2010) 

0.483 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

267 

85 

1 

Sposetti & Immler (2010) 

0.494 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

253 

80.4 

1 

Sposetti & Immler (2010) 

0.505 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

246 

78.2 

1 

Sposetti & Immler (2010) 

0.505 

NOT 

ALFOSC 

R 

4.56e-l-14 

277 

88.1 

1 

De Cia et al. (2010) 

0.516 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

232 

74 

1 

Sposetti & Immler (2010) 

0.527 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

210 

66.9 

1 

Sposetti & Immler (2010) 

0.539 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

218 

69.4 

1 

Sposetti & Immler (2010) 

0.551 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

216 

68.7 

1 

Sposetti & Immler (2010) 

0.562 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

218 

69.4 

1 

Sposetti & Immler (2010) 

0.574 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

206 

65.6 

1 

Sposetti & Immler (2010) 

0.585 

Gnosca 

0.4m 

CR 

4.4e-Hl4 

195 

62.1 

1 

Sposetti & Immler (2010) 

0.782 

KPNO 

SARA-N 

b 

6.92e-i-14 

135 

13 

1 

Updike et al. (2010b) 

0.785 

KPNO 

SARA-N 

V 

5.55e-Hl4 

234 

22.6 

1 

Updike et al. (2010b) 

0.787 

KPNO 

SARA-N 

R 

4.56e-i-14 

145 

14 

1 

Updike et al. (2010b) 

0.79 

KPNO 

SARA-N 

I 

3.93e-i-14 

148 

14.3 

1 

Updike et al. (2010b) 

0.856 

ISON-NM 

0.45m 

CR 

4.4e-Hl4 

147 

4.11 

1 

Elenin et al. (2010a) 
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Table 8 — Continued 


t-tQ 

(days) 

Observatory 

Telescope / 
Instnament 

Filter 

Frequency 

(Hz) 

Flux density^ 

(z^Jy) 

Uncertainty 

(fiJy) 

Detection? 

(1 = Yes) 

Reference 

1.05 

Ishigakijima 

Murikabushi 

Ic 

3.93e-Hl4 

121 

5.71 

1 

Kuroda et al. (2010b) 

1.05 

Ishigakijima 

Murikabushi 

Rc 

4.56e-i-14 

87.6 

3.29 

1 

Kuroda et al. (2010b) 

1.05 

Ishigakijima 

Murikabushi 

g’ 

6.29e-i-14 

58.1 

2.18 

1 

Kuroda et al. (2010b) 

1.07 

Hiroshima 

Kanata 

Rc 

4.56e-i-14 

82.1 

7.92 

1 

Yoshida et al. (2010) 

1.08 

Hiroshima 

Kanata 

Rc 

4.56e-i-14 

83.6 

5.57 

1 

Yoshida et al. (2010) 

1.09 

Ishigakijima 

Murikabushi 

Ic 

3.93e-i-14 

105 

5.94 

1 

Kuroda et al. (2010b) 

1.09 

Ishigakijima 

Murikabushi 

Rc 

4.56e-i-14 

82.1 

3.08 

1 

Kuroda et al. (2010b) 

1.09 

Ishigakijima 

Murikabushi 

g’ 

6.29e-i-14 

53 

1.99 

1 

Kuroda et al. (2010b) 

1.1 

Hiroshima 

Kanata 

Rc 

4.56e-i-14 

88.4 

3.32 

1 

Yoshida et al. (2010) 

1.14 

Ishigakijima 

Murikabushi 

Ic 

3.93e-i-14 

96.2 

6.41 

1 

Kuroda et al. (2010b) 

1.14 

Ishigakijima 

Murikabushi 

Rc 

4.56e-i-14 

69.6 

2.61 

1 

Kuroda et al. (201 Ob) 

1.14 

Ishigakijima 

Murikabushi 

g’ 

6.29e-i-14 

50.6 

2.38 

1 

Kuroda et al. (2010b) 

1.18 

Okayama 

MITSuME 

Ic 

3.93e-i-14 

96.2 

9.28 

1 

Kuroda et al. (2010d) 

1.18 

Okayama 

MITSuME 

Rc 

4.56e-i-14 

77 

7.43 

1 

Kuroda et al. (2010d) 

1.18 

Okayama 

MITSuME 

g’ 

6.29e-i-14 

63.1 

6.09 

1 

Kuroda et al. (2010d) 

1.22 

Hanle 

HCT 

R 

4.56e-i-14 

75.6 

7.29 

1 

Sahu et al. (2010a) 

1.24 

Hanle 

HCT 

I 

3.93e-i-14 

92.7 

13.7 

1 

Sahu et al. (2010a) 

1.3 

MASTER 


I 

3.93e-i-14 

90.2 

3.39 

1 

Gorbovskoy et al. (2012) 

1.38 

Tautenburg 

TLSI.34m 

Rc 

4.56e-i-14 

66.4 

1.86 

1 

Kann et al. (2010a) 

1.39 

Tautenburg 

TLSI.34m 

Rc 

4.56e-i-14 

61.7 

2.32 

1 

Kann et al. (2010a) 

1.42 

Tautenburg 

TLSI.34m 

Rc 

4.56e-i-14 

59.5 

1.67 

1 

Kann et al. (2010a) 

1.43 

Tautenburg 

TLSI.34m 

Rc 

4.56e-i-14 

55.3 

1.55 

1 

Kann et al. (2010a) 

1.43 

Tautenburg 

TLSI.34m 

Rc 

4.56e-i-14 

61.2 

2.29 

1 

Kann et al. (2010a) 

1.44 

Tautenburg 

TLSI.34m 

Rc 

4.56e-i-14 

65.2 

2.45 

1 

Kann et al. (2010a) 

1.46 

Terksol 

Z-600 

R 

4.56e-i-14 

62.9 

2.96 

1 

Andreev et al. (2010c) 

1.7 

Lightbuckets 

0.61m 

R 

4.56e-i-14 

46.8 

14.9 

0 

Ukwatta et al. (2010) 

2.13 

Ishigakijima 

Murikabushi 

Ic 

3.93e-i-14 

23.3 

3.95 

1 

Kuroda et al. (2010b) 

2.13 

Ishigakijima 

Murikabushi 

Rc 

4.56e-i-14 

31.8 

1.81 

1 

Kuroda et al. (2010b) 

2.13 

Ishigakijima 

Murikabushi 

g’ 

6.29e-i-14 

18.9 

1.44 

1 

Kuroda et al. (2010b) 

2.3 

CrAO 

AZT-II 

R 

4.56e-i-14 

28.7 

2.2 

1 

Rumyantsev et al. (2010) 

2.41 

Terksol 

Z-600 

R 

4.56e-i-14 

27.2 

3.18 

1 

Andreev et al. (2010d) 

2.43 

Terksol 

Z-600 

R 

4.56e-i-14 

23.7 

3.76 

1 

Andreev et al. (2010d) 

2.83 

ISON-NM 

0.45m 

CR 

4.4e-Hl4 

24.1 

6.25 

1 

Elenin et al. (2010b) 

2.88 

UKIRT 


K 

L37e-Hl4 

47.6 

9.64 

1 

Im et al. (2010) 

3.31 

Tautenburg 

TLSI.34m 

Rc 

4.56e-l-14 

18.8 

1.25 

1 

Kann et al. (2010a) 

3.55 

Tautenburg 

TLSI.34m 

Rc 

4.56e-l-14 

16.1 

0.604 

1 

Kann et al. (2010a) 

3.82 

ISON-NM 

0.45m 

CR 

4.4e-Hl4 

13.9 

2.81 

1 

Elenin et al. (201 Ob) 

4.26 

Ishigakijima 

Murikabushi 

Ic 

3.93e-l-14 

29.1 

9.25 

0 

Kuroda et al. (2010a) 

4.26 

Ishigakijima 

Murikabushi 

Rc 

4.56e-l-14 

9.26 

2.95 

0 

Kuroda et al. (2010a) 

4.26 

Ishigakijima 

Murikabushi 

g’ 

6.29e-i-14 

8.32 

2.65 

0 

Kuroda et al. (2010a) 

4.4 

Tautenburg 

TLS 1.34m 

Rc 

4.56e-i-14 

9.26 

0.436 

1 

Kann et al. (2010a) 

4.4 

CrAO 

AZT-11 

R 

4.56e-i-14 

11.1 

2.25 

1 

Rumyantsev et al. (2010) 

4.42 

Tautenburg 

TLS 1.34m 

Rc 

4.56e-i-14 

9.17 

2.48 

1 

Kann et al. (2010a) 

4.43 

Tautenburg 

TLS 1.34m 

Rc 

4.56e-i-14 

8.68 

1.1 

1 

Kann et al. (2010a) 

4.56 

Tautenburg 

TLSL34m 

Rc 

4.56e-i-14 

9.51 

0.918 

1 

Kann et al. (2010a) 

4.87 

ISON-NM 

0.45m 

CR 

4.4e-Hl4 

8.06 

1.81 

1 

Elenin et al. (2010b) 

5.55 

Tautenburg 

TLSL34m 

Rc 

4.56e-i-14 

5.28 

0.563 

1 

Kann et al. (2010b) 

6.3 

Maidanak 

AZT-22 

R 

4.56e-i-14 

4.64 

0.401 

1 

Voinova et al. (2010) 


Note. — ^In cases of non-detections, we report the formal flux measure¬ 
ment at the position of the afterglow. 

The X-ray data before the first orbital gap at around 10“^ d 
exhibits rapid flaring, ending in a steep decay. The earliest 
UVOT observations during this period detect a counterpart at 
multiple wavelengths. Owing to the contemporaneous 7-ray 
emission, indicating on-going central engine activity during 
this period, we do not consider data before the first orbital gap 
for our afterglow analysis. The X-ray re-brightening and sub¬ 
sequent fading between 0.08 d and 20 d can be well fit with a 
broken power law model (Equation 5), with t\, = 0.40 ± 0.02 d, 
fVx(fb) = 2.0±0.1/rJy, m = 0.77 ±0.15, 02 =-1-55 ±0.06, 
andy = 2.3±l.l (Figure 12). 

A broken power law fit to the UVOT t/-band data between 
0.02 d and 20 d yields fb = 0.40 ± 0.03 d, 7v,u(fb) = 144 ± 

7/iJy, ai = 0.45 ± 0.06, and a 2 = -1.63 ± 0.12, with y = 2.5, 
similar to the results from fitting the X-ray re-brightening. 
Over the same time range, a broken power law fit to the UVOT 
B-band data yields fb = 0.39 ±0.08d, /v.B(fb) = 201 ±24/iJy, 


ai = 0.40 ±0.06, and q ;2 = -1.46 ±0.21, similar to the X- 
ray and t/-band fits. Finally, we fit the X-ray, UVOT U-band 
and UVOT B-band data jointly, where we constrain the model 
light curves at the two frequencies to have the same rise and 
decay rate and time of peak, and use independent normaliza¬ 
tions in the three bands. Using a Markov Chain Monte Carlo 
simulation using EMCEE (Foreman-Mackey et al. 2013), we 
find fb = 0.36 ±0.02, ai = 0.57 ± 0.09, and 02 = -1.47 ±0.07, 
with a fixed value of y = 2.5. 

The spectral index between the UKIRT X-band observa¬ 
tion at 2.88 d (Im et al. 2010) and the X-rays is /3nir.x = 
-0.78 ±0.08. The spectral index between the NIR X-band 
and the Swift/U -band is slightly steeper, /?nir,uv = -1.0 ± 0.2, 
indicating that some extinction may be present. The spec¬ 
tral slope in the X-rays following the re-brightening is /3x = 
-1.15 ±0.06, leading to /3nir,x-/3x = 0.4 ±0.1 suggesting 
that f/opt <Vc<vx at 2.88 d. 

If the X-rays indeed lie above v^, we would have p = -2/3x = 
2.3 ± 0.12. The expected light curve decay rate is a = (2- 
3/7)/4 = -1.2±0.1 above v^: and a = 3(l-p)/4 = -1.0±0.1 
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Figure 12. Broken power law fit to the X-ray (grey, solid), Swift/WOT U- 
band (pink, dashed), and Swift/UYOT 5-band (blue, dotted) light curves for 
GRB 100901A near the re-brightening around 0.6 d. The X-ray fit includes 
points between 0.08 and 20 d, while the U- and 5-band fits include points 
between 0.02 d and 20 d. A joint fit yields the best fit parameters fb = 0-36 di 
0.02, ai =0.57 ±0.09, and 0^2 =-1.47 ±0.07 (Section 6.2.1). 

(ISM) or a = -1.5 ±0.1 (wind) below The t/-band light 
curve after the re-brightening can be fit either with a single 
power law with a = -1.5 or with a break from a « -0.9 to 
a ~ -1.8 at « Id and a = -1.0 ±0.1 between 0.2 d and 1.0 d. 
The former is consistent with a wind-like environment, while 
the latter suggests an ISM model and a jet break at Ri 1 d. We 
return to this point later. 

Extrapolating the Lightbuckets /?-band data at 1.7 d to the 
time of the SMA 345 GHz upper limit at 1.8 days using the 
fit to the optical t/-band light curve described above, we find 
a spectral index of /3mm-R ^ “0.6 from the millimeter to the 
optical at this time. Combined with the steeper spectrum of 
/3nir,x = -0.78 ±0.08 between the optical and X-rays, this 
suggests that the spectrum turns over above the millimeter 
band, indicating that > 345 GHz at 1.8 d. 

The spectral index between the VLA 4.5 GHz and 7.9 GHz 
bands at 4.92 d is /3 = 0.9 ± 0.3. However, the light curve at 
4.5 GHz declines as a = -0.17 ± 0.12 between 3 d and 12 d, 
and as a = -1.1 ±0.5 between 12 d and 20 d. The rising radio 
spectrum coupled with the declining light curve implies that 
the jet break has occurred before 5 d. The steepening of the 
4.5 GHz light curve at r; 12 d implies r 4.5 GHz at this 
time. Given that > 345 GHz at 1.8 d, this would imply an 
evolution of ~ r/m cx consistent with the expected evo¬ 
lution of I'm oc t~^ following the jet break. Together, this im¬ 
plies that the jet break must have occurred before 1.8 d, which 
agrees with the results from the X-ray and optical analysis 
above and indicates an ISM-like environment. We therefore 
focus on the ISM model for the remainder of this section, and 
discuss the wind model for completeness in Appendix D. 

To summarize, the NIR to X-ray SED exhibits mild ev¬ 
idence for extinction and suggests that i/nir < < I'x at 

R 3d, with p = 2.3 ±0.1. The radio observations indicate 
I'm ~ 5 GHz at R 12 d. Together, the X-ray, NIR and radio 
data suggest fjet Rid and an ISM-like environment. All ra¬ 
dio observations took place after this time, and are therefore 
insensitive to the density profile of the circumburst environ¬ 
ment. 

6.2.2. Forward shock model att> 0.25 d 



Figure 13. Observed NIR to X-ray SED of the afterglow of GRB 100901A at 
0.06 d (orange circles) and after the re-brightening (1.14 d; black squares), to¬ 
gether with the best-fit ISM model (Section 6.2.2) including energy injection 
(Section 6.2.3) at 0.06 d (orange; dashed) and 1.14d (black; solid). The data 
points and model at 0.06 d have been multiplied by a factor of 10. for clarity. 
The XRT SEDs at 1.14 d has been extrapolated from 3 d using the best-fit bro¬ 
ken power law model to the XRT light curve (Figure 12; the con'ection factor 
is R 4.5). The optical data have been extrapolated using the joint best fit to 
the SwiftKSH and UVOT light curves (the corrections are small at < 5%). 
The dotted curves indicate afterglow models with no extinction or IGM ab¬ 
sorption. Note that the Swi/f/UVOT tivwl-, uvw2-, and «vm2-band data lie 
blueward of Ly-o in the rest-frame of the host galaxy (1.0 X 10*^ Hz in the 
observer frame), and are likely subject to significant IGM absorption. For 
this reason, we do not include these bands (or the Swi/f/UVOT IWnYe-band, 
not shown here) to constrain our afterglow model. 


We employ our MCMC analysis to fit the multi-band data 
for GRB 100901A after t r 0.25 d. At the redshift of z = 
1.408, the UVOT White-, uvwl-, uvw2-, and uvni2-hanA 
data are affected by IGM absorption, and we do not in¬ 
clude these bands in our analysis. The highest-likelihood 
ISM model (Eigure 14) the parameters p r 2.03, R 0.33, 
cb « 0.32, iiQ R 3.2 X 10“^ cm“^, E^Aso ~ 3.0 x 10^^ erg, fjet R! 
0.96 d. Ay R 0.09 mag, and /v,host.r' R 4.1 /iJy, with a Comp¬ 
ton y-parameter of R 0.6. The blastwave Lorentz factor is 
r = 15.9(f/ld)“^/^ and the jet opening angle is 0jet R 2.1°. 
The beaming-corrected kinetic energy is Er ~ 2.1 x 10^°erg, 
while the beaming corrected 7-ray energy is r 5 x 10"^^ erg 
(l-10"^keV; rest frame). We plot histograms of the measured 
parameters in Eigure 15 and correlation contours between the 
physical parameters in Eigure 16, providing summary statis¬ 
tics from our MCMC analysis in Table 5. 

In this model, the break frequencies at 1 d are located at 
iZa R 0.7 GHz, R 2 X 10^2 Hz, and j/c R 9 x 10'^ Hz, while 
the peak flux density is about 4.2 mJy at I'm- The spectrum 
transitions from fast to slow cooling at about 350 s after the 
burst. The cooling frequency is located between the optical 
and X-rays and the jet break is before 1 d, both of which are 
expected from the basic analysis outlined above. The low 
value of p is driven by the shallow measured decline in the 
X-rays following the jet break, and the resulting spectrum re¬ 
mains consistent with the optical-to-X-ray SED (Eigure 13). 

6.2.3. Energy injection model 

Taking the forward shock model described in Section 6.2.2 
as a starting point, we hnd that the X-ray and UV/optical data 
before the re-brightening can be explained by a single period 












Energy injection in GRBs 


25 



Time (days) 



Time (days) 



10° io’ 

Time (days) 




Figure 14. X-ray, UV (top left), optical (top right), and radio (bottom left) light curves of GRB i00901A, with the full afterglow model (solid lines), including 
energy injection before 0.25 d. The X-ray data before 0.008 d exhibits strong flaring activity and we do not include these data in our analysis. The dashed 
envelopes around the radio light curves indicate the expected effect of scintillation at the Icr level. The data indicated by open symbols are not used to determine 
the pai'ameters of the forwai'd shock (the MCMC analysis). Bottom right: blastwave Lorentz factor (green, dashed; upper sub-panel) and isotropic equivalent 
kinetic energy (red, solid; upper sub-panel) as a function of time, together with the energy distribution across ejecta Lorentz factors (black, solid; lower sub-panel) 
as determined from fitting the X-ray/UV/optical re-brightening at 0.36 d. 







Figure 16. Icr (red), 2(j (green), and 3cr (black) contours for correlations 
between the physical parameters, and for GRB 100901A, in 

the ISM model from Monte Carlo simulations. We have restricted ce < 
and CB < V^- version of this Figure for additional plots of 

congelations between these parameters and p, (jet, ^jet^ ^K, -^v, and iv,host,r' • 
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Figure 15. Posterior probability density functions for the physical pai'ame¬ 
ters for GRB 100901A in a constant density environment from MCMC sim¬ 
ulations. We have restricted £'K,iso 52 < 500, Ce < i/3, and cb < 1/3. 
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In this model, the energy increases by a factor of Ri 5.5 from 
^Kasoj « 5.4 X 10^^ erg at 0.105 d tofiKjso.f« 3.0 x 10^^ erg at 
0.26 d, corresponding to an injected energy fraction of « 85%, 
while the Lorentz factor decreases from F « 30 to F « 26 
over this period. In comparison, E^^iso « 8 x 10 ^^ erg (Section 
6.2.1). The value of m « 1.88 corresponds to i r: 14.4 for 
26 < F < 30. 


6.3. GRB 120404A 
6.3.1. GRB properties 

GRB 120404A was detected and localized by the Swift 
BAT on 2012 April 04 at 12:51:02 UT (Stratta et al. 2012). 
The burst duration was Tgo = 38.7 ±4.1s, with a flu- 
ence of = (1.6 ±0.1) x 10“^ergcm“^ (15-150keV ob¬ 
server frame, Ukwatta et al. 2012). Swift and ground- 
based observatories detected an afterglow in the X-rays 
and UV/optical (Guidorzi et al. 2012; Xin et al. 2012b; 
Breeveld & Stratta 2012; Gorbovskoy et al. 2012), as well 
as in the radio (Zauderer et al. 2012). Spectroscopic obser¬ 
vations with Gemini-North yielded a redshift of z = 2.876 
(Cucchiara & Tanvir 2012). The isotropic-equivalent 7-ray 
energy for this event is E^ iso = (9 ±4) x 10^^ erg (1-10^keV, 
rest frame; Guidorzi et al. 2013). 

This burst has been previously studied in detail by 
Guidorzi et al. (2013), who interpret the optical re¬ 
brightening starting around 800 s in their well-sampled, 
multi-band optical light curves as due to the passage of the 
characteristic synchrotron frequency, z/m. They additionally 
invoke reverse shock emission to explain the flat (f0 0±0') 
portion of the optical light curve before the onset of the 
re-brightening. In the following, we propose an alternate 
model for the multi-band radio through X-ray light curves in 
the context of energy injection. 

The X-ray light curve before 700 s can be modeled as a 
power law decay with ax ~ “2. This light curve phase is 
likely part of the high latitude emission, and we ignore the 
data before 0.008 d in our analysis. The X-ray photon index 
at 0.12-0.24d, Fx = 2.3 ±0.3 (Guidorzi et al. 2013), imply¬ 
ing a spectrum, F^, oc jy-' ^±0.3^ consistent with the spectral 
slope between the optical /'-band and the X-rays at 0.07 d, 
/3opt-x = -0.91 ±0.04, suggesting that the optical and X-ray 
bands are on the same segment of the afterglow SED, al¬ 
though the large uncertainty in /3x leaves open the possibil¬ 
ity that z/c lies between the optical and X-rays. Addition¬ 
ally, the spectral slope within the optical (B- to z'-band) is 
/3opt = -1.3 ±0.2 at 0.07 d, indicating that extinction is present. 
The spectral index between the 19.2 GHz and 24.5 GHz obser¬ 
vations at 0.75 d is ^Sradio ~ 2, which indicates that z/a is located 
above 24.5 GHz at this time (Figure 17). 

The optical R-band light curve declines as ^±0.02 
0.13 d, consistent with the X-ray decline rate of ^±0.1 
ter 0.05 d. The steep decline of a R -2 is indicative of a jet 
break before r 0.1 d. A broken power law fit to the B-band 
light curve results in the parameters fb = (2.8 ±0.4) x 10“^d, 
7v,B(fb) = 252± lOpJy, ai = 1.74 ±0.58, a 2 = -1.71 ±0.17, 
and y = 0.78 ±0.43, making this the earliest re-brightening 
episode of the four events studied in this paper. 

6.3.2. Forward shock model at t >0.04 d 

We interpret the optical light curve peaks at around 2500 s 
as the end of a period of energy injection, after which the af¬ 
terglow evolves according to the standard framework with a 


flxed energy, TiK.iso.f- We use the data after 0.04 d for esti¬ 
mating the parameters of the blastwave shock and employ our 
MCMC tools described in Section 4 to model the afterglow 
after this time. 

The parameters of our highest likelihood ISM model are 
p R 2.06, Ce « 0.27, eg r 0.16, hq r 2.8 x 10^cm“^, isK.iso ~ 
1.2 X 10^^ erg, fjet R! 6.6 x 10“^ d, and Ay r 0.13 mag. The 
Compton y-parameter for this model is r 0.9. The blast- 
wave Lorentz factor is F = 4.1 (f /1 d)“^/* and the jet opening 
angle is 0jet R 3.1°. The beaming-corrected kinetic energy 
is Ex ~ 1-7 X 10^’’erg, while the beaming corrected 7-ray 
energy is E’.y r 1.2 X 10^*^ erg (1-10"^ keV; rest frame). The 
MCMC analysis yields an isotropic equivalent kinetic energy 
of £’K,iso= 1-3 !!!o 2 X 10^^erg, which is similar to the value of 
1.9;!;oj X 10^^ erg derived by Guidorzi et al. (2013). The high 
circumburst density of log(zzo) = 2.5 °q 3 is also consistent with 
the value of log(no) = 2.4^g 2^ determined by Guidorzi et al. 
2013), and is driven by the low flux density and steep spec¬ 
trum in the radio. However, unlike Guidorzi et al. (2013), our 
low value of p r 2.1 allows us to match the NIR to optical 
SED and the normalization of the X-ray light curve. We note 
that the high circumburst density also results in a low cool¬ 
ing'^ frequency: in our highest likelihood model 
and the spectrum remains in the fast cooling regime through 
the duration of the observations (spectrum 4 of Granot & Sari 
2002), with z^ac R 5.9 x 10'°Hz, z^sa R 2.0 x 10" Hz, z/m R 
9.1 X 10" Hz, andi/c r6.3 x 10'“Hz at 1 d. 

We plot the posterior density functions for the all param¬ 
eters in Figure 18 and correlation contours in Figure 19; we 
list our measured values for the physical parameters in Table 
5. In the wind and ISM models, light curves at all frequencies 
become indistinguishable following a jet break. Since the jet 
break occurs early, soon after the start of the data used for 
deriving the parameters of the blastwave, we expect a viable 
wind model to exist as well. We discuss this model in Ap¬ 
pendix E. 


6.3.3. Energy injection model 

We model the light curves before to = 0.04 d by injection of 
energy into the blastwave shock. We use the afterglow param¬ 
eters for the highest-likelihood model (Section 6.3.2) as the 
final parameters following the end of energy injection, with 
FK,iso,f = 1.2 X 10^^ erg. We And that the optical and X-ray 
light curves can be modeled well by two subsequent periods 
of energy injection, beginning at fa = 7 x 10“^ d (Figure 20): 



f > fo = 0.04d 
fi = 2.2 X 10“^d < f < fo 
f2 = 7.0 X 10“^d < f < fi 
f < fa 

( 6 ) 


In this model, isK.iso increases by a factor of R 18 from 
4.4 X 10^'erg to 7.8 x 10^“ erg between 7.0 x 10“^ d and 
2.2 X 10“^ d, and then by r 50% toFK.iso.f R 1.2 x 10^^ erg be¬ 
tween 2.2 X 10“^ d and 0.04 d. Thus r 95% of the Anal kinetic 
energy is injected into the blastwave between 7.0 x 10“^ d and 


Once again, Uc = yfFpFn, with z /3 R 8.0 x 10*** Hz, and an R 5.0 x 
10*° Hz at 1 d (see Footnote 12). 
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Figure 17. Left: Observer-frame SED of the afterglow of GRB 120404A at 0.0664 d, together with the best fit model (spectrum 4 from GS02: black, solid). 
Right: The same as the left panel, but at 0.78 d. The radio detections with a steeply rising spectrum suggest that the self-absorption frequency is near or above 
the radio band, requiring a high circumburst density. Both SEDs are corrected for the effect of inverse Compton cooling, which is significant, with the Compton 
y-parameter 0.9 during the fast cooling phase. We do not show the inverse Compton radiation, since even at its peak (near the SwiJt/XRT band), it is four 
orders of magnitude weaker than the synchrotron component. See Appendix A for a discussion of the inverse Compton effect in the context of modeling GRB 
afterglows. 









Figure 18. Posterior probability density functions for the physical parame¬ 
ters for GRB 120404A in a constant density environment from MCMC sim¬ 
ulations. We have restricted £K,iso ,52 < 500, ee < and ea < 
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Figure 19. Icr (red), 2a (green), and 3a (black) contours for correlations 
between the physical parameters, £K,iso. f>o, ee, and eg for GRB 120404A, in 
the ISM model from Monte Carlo simulations. We have restricted £K.iso,52 < 
500, ee < 1/3, and ea < ^/3. See the on line version of this Figure for addi¬ 
tional plots of correlations between these parameters and p, fjet, 0 jet, E^, and 
Av. 


r « 13.7 at the end of energy injection at 0.04 d. The value of 
m derived above corresponds to s ss 3.7 for 13.7<r<16.2 
and 5 Ri 40 for 16.2 < T < 17.4. 

7. DISCUSSION 
7.1. Parameter distributions 

We now turn to the question of how GRBs that exhibit si¬ 
multaneous optical and X-ray re-brightening episodes com¬ 
pare with each other, and with events that do not exhibit 
such a feature. For this purpose, we use the compila¬ 
tion of measurements of E^, 0jet, Ek, and no (or A*) from 
the literature (Panaitescu & Kumar 2002; Yost et al. 2003; 
Friedman & Bloom 2005; Ghirlanda et al. 2007; Cenko et al. 
2010, 2011) reported in Laskaretal. (2014). This sample 
includes GRBs from the pre-Swift era, as well as Swift and 
Fermi events. 

The radio to X-ray observations of all four GRBs presented 
here can be fit by constant density ISM models. The best fit 
models yield densities from a; 3 X 10“^ to Ri 3 X 10^ cm ^and 
final beaming-corrected kinetic energies from Ri 2 x 10 ^*^ erg 
to Ri 2 X 10^* erg, spanning the full range of values inferred 
from GRB afterglow modeling (Figure 21). We constrain the 
jet break time and hence the opening angle of the GRB jet in 
each case, and find 0 jet ~ 2 °-21 °, spanning the range inferred 
from the comparison sample. 

The beaming-corrected 7 -ray energies in the 1-10'* keV rest 
frame energy band of the events in our sample range from 
5 X lO'^^erg to 1.3 x 10^** erg, while the median'® beaming 
corrected 7 -ray energy of the comparison sample is = 
(8.1)^ }) X 10^" erg (95% confidence interval. Figure 22). 
Therefore, the observed values of F-y for the events in our 
sample are all smaller than the best estimate for the median of 
the comparison sample. 

To further quantify this effect, we compute the radiative ef¬ 
ficiency, 7] = Ej^/{Ej-\- Ejf) for each GRB and in the com¬ 
parison sample. Since both £k and have associated un- 


0.04d. In comparison, iso « lO^-* erg is similar to the final 
kinetic energy. The blastwave Lorentz factor decreases from 
r Ri 17.4 at 7.3 X 10“^ to T R3 16.2 at 2.2 x 10“^ d, and then to 


The uncertainty on the median is computed using Greenwood’s formula 
for the variance of the Kaplan-Meier estimate of the cumulative distribution 
function. This method accounts for both upper and lower limits, which exist 
in the data. 
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Figure 20. X-ray, UV (top left), optical (top right), and radio (bottom left) light curves of GRB i20404A, with the full afterglow model (solid lines), including 
energy injection before 0.04 d. The X-ray data before 0.008 d is likely dominated by high-latitude prompt emission and we do not include these data in our 
analysis; the best fit power law to the X-ray data before 0.008 d added to the blastwave model is shown in the upper left panel (black, dashed). The dashed 
envelopes around the radio light curves indicate the expected effect of scintillation at the Icr level. The data indicated by open symbols are not used to determine 
the parameters of the forward shock (the MCMC analysis). The t/-band data are strongly affected by IGM absorption and are not included in the fit. Bottom right: 
blastwave Lorentz factor (green, dashed; upper sub-panel) and isotropic equivalent kinetic energy (red, solid; upper sub-panel) as a function of time, together 
with the energy distribution across ejecta Lorentz factors (black, solid; lower sub-panel) as determined from fitting the X-rayAJV/optical re-brightening at 0.04 d. 


certainties, while the expression for rj is non-linear in these 
two quantities, a proper accounting of the final uncertainty re¬ 
quires a Monte Carlo analysis. We generate 10^ Monte Carlo 
realizations of Ek and for each burst, assuming a uniform 
distribution^® between the Icr error bars. We then compute 
and summarize the resulting distribution of p using the me¬ 
dian and 68 % credible intervals (Figure 23). We find that the 
GRBs in our sample have systematically lower radiative effi¬ 
ciencies. This is consistent with the energy injection scenario, 
if the prompt 7 -ray radiation is dominated by emission from 
the fast-moving ejecta, while a significant amount of kinetic 
energy is carried by slow-moving ejecta. 

From this comparison, we conclude that events exhibiting 
simultaneous, multi-wavelength re-brightening episodes also 
span the same wide range of circumburst densities, jet open¬ 
ing angles, and beaming-corrected kinetic energies inferred 

We choose a uniform distribution instead of, say, a Gaussian distribution 
because the uncertainties on these parameters are frequently large compared 
to the mean, and a Gaussian distribution in linear space results in a tail of un¬ 
physical, negative values. A more detailed analysis requires the full posterior 
density functions for both and Ek for every GRB, which are not available. 


from previous studies of events at z ^ 1 and z > 6 . However, 
the events reported here have smaller beaming-corrected 7 - 
ray energies, and as a result lower radiative efficiencies than 
the comparison sample, suggesting that the prompt radiation 
is dominated by ejecta at high Lorentz factors, while the bulk 
kinetic energy is dominated by slow-moving ejecta at least in 
these cases. 

7.2. Absence of reverse shock 

When energy injection into the blastwave is caused by a 
distribution of ejecta Lorentz factors, the reverse shock (RS) 
from the initial interaction of the leading edge of the ejecta 
with the circumburst medium is expected to continue to prop¬ 
agate through the ejecta until the end of the period of energy 
injection (Sari & Meszaros 2000). During the period of en¬ 
ergy injection, the afterglow SED is expected to be composed 
of contributions from both forward (FS) and reverse (RS) 
shocks, each with its three characteristic frequencies and flux 
normalizations. The spectral characteristics of the two shocks 
are expected to be related by /v^niax,Rs(0 = F^,max,Fs(t) x r(f), 
l^m.Rs(t)= l^m,Fs(0/r^(f). ^nd r'c,Rs(0 = t'c ,Fs(0. while the two 
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Figure 21. Beaming-corrected kinetic energy (left) and circumburst density (right) for both ISM (black circles) and wind-like environments (grey squares). The 
four GRBs in our analysis 100418A (blue), 100901A (red), 120326AA (green), and 120404A (purple), do not appear distinct from the comparison sample (grey 
and black; Panaitescu & Kumai' 2002; Yost et al. 2003; Chandra et al. 2008; Cenko et al. 2010, 2011 ; Laskar et al. 2014). 
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self-absorption frequencies should be related by ,Rs (0 = 
t'a.FsCO X r^(f), where A = 3/5 when both shocks are fast cool¬ 
ing and A = ^/s when both shocks are slow cooling. The re¬ 
verse shock is expected to last through the period of energy 
injection, whereupon all ejecta have been decelerated to a 
common Lorentz factor, and the two shocks decouple in their 
subsequent evolution. 

A calculation of the reverse shock SED based on the above 
considerations, together with the Lorentz factor of the blast- 
wave at the end of the period of energy injection (fdec) and 
the SED of the forward shock at that time indicates that 
the reverse shock generically peaks around the millimeter 
band at fdec (Sari & Meszaros 2000). Eor GRB 100901A and 
120404A, there are no data in the millimeter bands at the 
relevant time. Eor GRBs 120326A and 100418A where we 
have millimeter observations at Ri fdec, we find that the RS 
light curves over-predict the observations by factors of Ri 2 - 
10. Additionally, the reverse shock would also contribute a 
flux density comparable to the forward shock in the X-rays 
and optical for all four events during energy injection, which 
would require suppressing the forward shock before the re¬ 
brightening by invoking even lower blastwave energies (and 
consequently requiring an even greater rate of energy injec¬ 
tion) before the optical/X-ray peak. Thus the data are incon¬ 
sistent with the presence of a strong reverse shock for these 
events. 

We note that in our model there is a gap before the begin¬ 
ning of energy injection. If the energy injection is caused by 
a shell of material with a distribution of Lorentz factors (the 
‘injective shell’) catching up with the initial shell (‘the impul¬ 
sive shell’), then a long-lasting reverse shock in the injective 
shell is only expected when the two shells collide violently, 
with a relative Lorentz factor greater that the sound speed in 
the injective shell (Zhang & Meszaros 2002). If the injective 
shell is released from the central engine at roughly the same 
time as the impulsive shell, the collision between the two will 


be mild, and a reverse shock will not occur. This appears to be 
the case for the GRBs studied here, suggesting that the shell 
collisions in these events exhibiting X-ray and UV/optical re¬ 
brightening episodes are gentle, resulting simply in a transfer 
of energy to the blastwave. Even more fundamentally, this 
implies that the shells are emitted at the same time from near 
the central engine and thus the engine need not be on during 
the re-brightening episode. 

7.3. Energy injection: a ubiquitous phenomenon? 

The phenomenon of short-lived plateaus at an early time is 
ubiquitous in SwiftlXKY afterglow light curves (Liang et al. 
2007). Using a complete sample of SwiftlXKY light curves 
through 2010, Margutti et al. (2013) find about 37% of long- 
duration GRBs exhibit a shallow decay phase, with -1 a < 
1. If these plateaus are caused by injection of energy into the 
forward shock, the events with X-ray/optical re-brightenings 
discussed in this paper might be the extreme tail of a distribu¬ 
tion in energy injection factor, duration or rate. A rigorous 
exploration of these possibilities requires multi-wavelength 
fits to the data, but such data are generally not available. We 
therefore compute energy injection fractions based on the X- 
ray light curves alone, and compare the results for the objects 
where we do have multi-wavelength observations as reported 
in this paper. 

Margutti et al. (2013) measured the timing of X-ray 
plateaus in SwiftlXKY light curves and quantified them by 
their start time, fi, end time, fa, and rise rate^' during the 
plateau phase, a. Eor our X-ray plateau analysis, we se¬ 
lect the 96 events in their study where the reported uncer¬ 
tainty in a is < 0.3. We repeat their light curve decompo¬ 
sition analysis on the X-ray light curves of GRB 120326A 
and GRB 120404A (which are not included in their sample), 
and find h = 5.7 x 10-^ d, fa = 0.62d, and a = 0.27 ±0.03, 

Margutti et al. (2013) use the convention F{t) oc , which is opposite 
to the convention in this paper. 
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Figure 22. Beaming-con'ected 7 -ray energy (left) and jet opening angle (right) for the events in our analysis, GRB 100418A (blue), 100901A (red), 120326AA 
(green), and 120404A (purple), together with a comparison sample of long GRBs (grey; Friedman & Bloom 2005; Ghirlanda et al. 2007; Cenko et al. 2010, 2011 ; 
Laskar et al. 2014). The isotropic-equivalent 7 -ray energy for GRB 050904 is taken from Amati et al. (2008), and for GRB 090423 from Salvaterra et al. (2009). 
The four GRBs exhibiting X-ray and optical re-brightening episodes do not appeal* distinct from the comparison sample in but appear to all reside at lower 
values of £7 than the median. 


for GRB 120326A and h = 9.8 x 10“^ d, fz = 2.9 x IQ-^ d, and 
a = 0.76 ± 0.23 for GRB 120404 A 22 
We define the plateau duration, AT = t2-t\ and the frac¬ 
tional duration, ^ = (t 2 -t\)/(t 2 + t\) € [0,1). Note that ^ 1 

when t 2 ^0, which corresponds to the case that the plateau 
lasts much longer than its onset time. Margutti et al. (2013) 
computed the start time of the plateau using the intersection 
of the steep decay phase with the best fit plateau model; hence 
they do not report an error on this quantity. We can estimate 

The X-ray light curve of GRB 120404A is not well-sampled near the 
peak of the re-brightening and we therefore fix the peak time in the fit to the 
value inferred from the optical light curve (p^ 0.3 d; Section 6.3.3). 


the uncertainty on ^ using uj = at 2 d^/dt 2 = 2fitT,,/(0 +( 2 )^, 
which is dimensionless as desired. 

Assuming the X-ray band is located above the cooling fre¬ 
quency for all cases (which is indeed the case for the five 
events considered in detail here), and that p = 2, we can com¬ 
pute an effective energy injection rate, E oc f™, where m = a+l 
is the rate required to bring the measured light curve slope to 
a from the theoretically-expected value^^ of f(2-3p)/4 ^ 
no energy injection. For the sake of simplicity and unifor- 

In the general case, m = ^ ™ 

(ISM) and m = (wind) when Uc> ux- 
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Figure 23. Radiative efficiency, rj for the events in our analysis, GRB 
100418A (blue), 100901A (red), 120326AA (green), and 120404A (purple), 
together with a comparison sample of long GRBs (grey; Friedman & Bloom 
2005; Ghirlanda et al. 2007; Cenko et al. 2010, 2011; Laskaretal. 2014), 
with 68% credible regions (error bars) about the median (points) computed 
from £'k and using a Monte Carlo procedure. The four GRBs exhibiting 
X-ray and optical re-brightening episodes have lower radiative efficiencies 
than the median. 

mity, we use p=2 even where we have other measurements 
of p, such as for the events reported in this paper. The ra¬ 
tio of the energy at the end of the plateau phase to the en¬ 
ergy at the beginning of the plateau phase is then given by 

T=£2/£i = feAir. "" 

We plot the plateau fractional duration, ^ against the plateau 
slope, a in Figure 24, scaling the area of the symbols by T. 
The slopes during the plateau phase range from a Ri -1 to 
PS 1.8. As evident from the kernel density estimate of a, most 
light curves exhibit a gentle decay during the plateau phase 
(Figure 24). However, all four events in our sample are ex¬ 
ceptions to this rule"^. These four events also exhibit the 
greatest rise in their plateau phase light curves for a given nor¬ 
malized duration, Two out of these four (GRBs 120326A 
and 100418A) have the largest fractional changes in the ki¬ 
netic energy of the entire sample of 98 events. These are 
also the events with the largest value of ^ in our sample. For 
GRB 120404A, we note that the X-ray light curve around the 
peak of the re-brightening is missing due to a Swift orbital 
gap. Therefore any results that rely solely on the X-ray data 
of this burst underestimate the value of a and the fraction of 
energy injection relative to the multi-wavelength analysis we 
carried out in Section 6.3.3. 

The uncertainty in this quantity can be estimated using 
= (dT/dt 2 fal +(dT/dmfal 
= (mTcr,^/t2f + (Tln(t2/ti)!^af, 

where we have taken a,,] = era. This assumes that and ol are independent. 
However, these quantities are expected to be correlated, and a more complete 
analysis of the uncertainty in T requires the full covariance matrix between 
OL and t 2 . 

Since we selected GRBs that exhibit an X-ray re-brightening, this is 
partly by sample construction. 


120326A 100418A 



Figure 24. The fractional duration of X-ray plateaus from Margutti et al. 
(2013) as a function of the light curve rise rate, a (with the convention 
F cc r“). Errors in both directions are computed according to the formulas 
described in section 7.3. The vertical dashed black line indicates the canon¬ 
ical light curve decay of a 5^ —1 expected for p 2. 21% of events with 
plateaus exhibit re-brightenings (a > 0). Events analyzed in this paper are 
shown in red. The ai'ea of each symbol is linearly proportional to the ratio of 
the energy at the end to that at the beginning of the plateau phase, T. The 
blue point aX —1.5 is GRB 081028, which also exhibits an X-ray/optical 
re-brightening; however, this event does not have radio data and we therefore 
exclude it from our multi-wavelength analysis. A probability density func¬ 
tion of the distribution of a (computed from a kernel density estimate using 
(Ta as a varying kernel bandwidth) is shown in cyan. 


In Figure 25 we plot the ratio of final to initial energy, T 
against the normalized plateau duration, scaling the ra¬ 
dius^® of the symbols by a. We find that for a given (nor¬ 
malized) plateau duration, the events in our sample have the 
largest fractional change in blastwave kinetic energy during 
the plateau phase. At the same time, they also possess the 
steepest rise rates in the sample, which is simply indicative 
of our selection criteria for inclusion in the present analy¬ 
sis. Finally, these events have the shortest normalized plateau 
durations for a given fractional change in energy. Together, 
these observations suggest that a large amount of injected en¬ 
ergy is not sufficient to cause simultaneous -ray/optical re- 
brightenings, but that it jnust be done in a relatively small 
amount of time, and thus that the defining characteristic of 
these events is a high rate of energy injection. Physically, this 
translates to a steep distribution of ejecta Lorentz factors over 
a small range of T. Hascoet et al. (2015) recently suggested 
that the interaction of low-Lorentz factor ejecta with the re¬ 
verse shock can explain flares in the X-ray light curves. Our 
observations provide supporting evidence for this hypothesis 
in the form of significant ejecta energy down to the requisite 
low Lorentz factors, T ~ 10. 

7.4. Energetics 

Having compared the plateau duration, light curve rise rate, 
and injected energy fraction for our sample with a complete 
sample of Swift events using the X-ray data alone, we now 
turn to an analysis of the results from our multi-wavelength 
energy injection modeling in this context. In Figure 26 we 
plot the fractional change in energy determined from multi¬ 
wavelength modeling against the fractional change in energy. 

Since o can be negative while the area of the symbols in the plot is a 
positive definite quantity, we scale the radius of the symbols as q:+ 1.5, where 
the additive term accounts for the minimum value of a. 
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Figure 25. The ratio of energy at the end to that at the beginning of the 
plateau phase in X-ray light curves of 96 events from Margutti et al. (2013), 
as well as GRBs 120326A and 120404A from our analysis, as a function of 
the plateau fractional duration, with en'or bai's computed according to the 
prescription in section 7.3. The size of the symbols is scaled to a measure of 
the light curve rise rate during the plateau decay phase, a, with the radius of 
the symbols proportional to q; + 1.5. Larger symbols denote a greater depai- 
ture from a canonical light curve decay of a 5^ -1, indicating greater rates of 
energy injection during the plateau phase. 



Figure 26. Ratio of final energy to the energy prior to the first injection 
episode from the multi-wavelength analysis (y-axis) compared to the ratio of 
the energy at the end to that at the start of the ‘plateau phase’, T, computed 
in section 7.3 for the five events in our sample, together with the 1:1 line 
(solid, black). The area of the symbols is linearly proportional to the physical 
duration of the plateaus, — 


T, inferred from the X-ray-only analysis in section 7.3 for 
the four events in our sample. We scale the area of the sym¬ 
bols with the plateau physical duration, f 2 - fi. As expected, 
events with longer plateau durations have greater estimated 
energy injection fractions both from the X-ray analysis and 
from the full multi-wavelength study. The fractional change 
in energy from the X-ray analysis is higher than inferred from 
the full model for three out of four cases. This is due to the 
typically lower value of m from the X-ray-only analysis com¬ 
pared to the full model. Ultimately, this can be traced to the 
fitting procedure: Margutti et al. (2013) fit a sum of a steep 
decay and a rising light curve, the sum of which results in the 
X-ray plateau, whereas we do not subtract fits to the steep de¬ 
cay phase from the X-ray light curve prior to our multi-band 


modeling. In the case of GRB 120326A, where the X-ray 
analysis over-estimates T by an order of magnitude, the vari¬ 
ance between the two techniques is due to differences in the 
plateau start time (4.3 x 10“^ d for the multi-wavelength anal¬ 
ysis, compared to 5.7 x 10“^ d in the X-ray analysis). Based 
on this comparison, we conclude that a determination of T 
using the X-ray data alone yields a reasonable estimate of the 
energy injection factor on average, but that multi-wavelength 
modeling is essential to obtain the full picture. 

8 . SUMMARY AND CONCLUSIONS 

Plateaus and re-brightening episodes are frequent in GRB 
afterglow X-ray light curves. However, X-ray data by them¬ 
selves provide only a limited understanding of the physi¬ 
cal processes underlying these unexpected phenomena. We 
perform a thorough multi-wavelength study of all long- 
duration GRBs through 2012 featuring simultaneous X-ray, 
and UV/Optical re-brightenings and that have radio detec¬ 
tions, using a Markov Chain Monte Carlo analysis and a phys¬ 
ical afterglow model. Our analysis yields the first set of mod¬ 
els that explain the multi-wavelength afterglows for all of 
these events with the same unifying principle. In all cases, 
the afterglow light curves can be modeled with a standard 
forward shock model following the re-brightening episode, 
an ISM density profile, and a jet break. From our multi¬ 
wavelength analysis, we find that the circumburst densities, 
jet opening angles, and beaming-corrected kinetic energies for 
these events span the full range described by typical GRBs at 
z ~ 1 and z > 6 . 

We explore a range of possible models to understand the 
re-brightenings, including the onset of the afterglow, off-axis 
viewing geometry, and continuous energy injection. We are 
able to rule out the afterglow onset and off-axis jets, and find 
instead that injection of energy into the blastwave (the so- 
called ‘refreshed-shock’ scenario) provides a good explana¬ 
tion for all events. We interpret energy injection in the frame¬ 
work of the stratified Lorentz factor model, and find that our 
measured energy injection rates always obey the theoretical 
constraints relating the rate of injection and the distribution of 
ejecta Lorentz factors in the ISM model. 

We perform the first measurement of the ejecta Lorentz fac¬ 
tor distribution index, s, and find s ra 3^0, suggesting that a 
large amount of kinetic energy resides in the slowest-moving 
ejecta. This is supported by low radiative efficiencies for the 
events in our sample, indicating that the 7 -ray radiation is 
dominated by ejecta at high Lorentz factors, while the kinetic 
energy is dominated by slower-moving ejecta. We note that 
keeping the injection rates simple power laws allows us to di¬ 
rectly convert the injection rate to a Lorentz factor distribution 
index, but the true injection rate and also the true Lorentz fac¬ 
tor distribution is likely to be a smoother function of time and 
Lorentz factor, respectively. 

Finally, using a compilation of X-ray plateaus in GRB after¬ 
glows, we present a comparative discussion of this interesting 
sub-population of GRBs. We find that the phenomenon of en¬ 
ergy injection is ubiquitous is long-duration GRBs, with re¬ 
brightening episodes likely simply extreme injection events. 
In future work, we aim to fit the light curves before the re¬ 
brightening episodes in a statistical sense, allowing us to es¬ 
timate uncertainties on the rate and duration of the energy in¬ 
jection episodes. At the same time, radio monitoring of Swift 
events exhibiting re-brightening episodes will be crucial in 
multi-wavelength modeling of this interesting class of GRBs, 
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while ALMA observations will irrefutably establish the pres¬ 
ence or absence of reverse shocks, laying to rest the question 
of whether the energy injection process is violent or gentle. 
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APPENDIX 

A. INVERSE COMPTON CORRECTION 

The synchrotron photons produced in the GRB blast- 
wave can Compton-scatter off the shock-accelerated rel¬ 
ativistic electrons in the blastwave, producing a Comp- 
tonized spectrum at high energies (Panaitescu & Meszaros 
1998; Wei&Lu 1998; Totani 1998; Chiang & Dermer 
1999; Dermer etal. 2000a,b; Panaitescu & Kumar 2000; 
Blandford & McKee 1977; Sari & Esin 2001). We com¬ 
pute the IC spectrum by directly integrating the synchrotron 
spectrum over the electron Lorentz factor distribution (back- 
calculated from the synchrotron spectrum) using equation A2 


in Sari & Esin (2001). We find that IC emission contributes 
negligible flux compared to synchrotron radiation and we ig¬ 
nore this component in our analysis. However, the IC mecha¬ 
nism can provide a significant source of cooling for the shock- 
accelerated electrons and thereby dominate the total cooling 
rate, even when IC emission itself is not directly observable 
(Sari et al. 1996; Sari & Esin 2001; Zhang et al. 2007). This 
effects the synchrotron cooling frequency, as well as the self- 
absorption frequency (when r'c < i/m). Thus IC cooling should 
be taken into account when computing SEDs and light curves 
for the synchrotron component. 

The importance of IC cooling is determined by the Comp¬ 
ton y-parameter^^. 


y = (Al) 

where 77 is th? fraction of energy that has been radiated away 
due to synchrotron and IC radiation, such that 77 = 1 dur¬ 
ing fast cooling and 77 = (r'c/r'm) during slow cooling 

(Sari & Esin 2001). Writing where ly' is the 

cooling frequency of the synchrotron SED not corrected for 
IC cooling, we have 


ly'd+y) 


-21 


-{ p - 2)/2 


=// 


1 + 


p-2 


-{p-2)l2 


(A2) 


where H = is independent of 77 . We therefore 

obtain the following implicit equation for rj. 




l + y/l+4jjee/eB 
2 


= 0, 


(A3) 


which can be solved numerically using (for instance) the 
Newton-Raphson method. Finally, the effect of IC cooling 
can be accounted for by scaling the spectral break frequencies 
and flux densities of the synchrotron spectrum by the appro¬ 
priate powers of 1 +Y (Granot & Sari 2002). For convenience, 
we summarize these scaling relations in Table 9. 


Table 9 

Inverse Compton Corrections to Specti'al Break Frequencies 


Spectral Break 

Break Frequency 

Break type 

Frequency scaling 

Flux density scaling 

1 

l^sa 

Self-absorption 

(i+y)° 

(l+F)'’ 

2 

t^ni 

Characteristic 

(1+T)° 

(1+F)“ 

3 


Cooling 

(i+yr2 

(1+F)P-' 

4 


Characteristic 

(1 + F)° 

(1+F)° 

5 


Self-absorption 

(1 + F)° 

(1+F)“ 

6 

t^sa 

Self-absoiption 

(1+F)?S 

(1 + F)^ 

7 

^ac 

Self-absoiption 

(l+FpS/s 

(1+F)-6/5 

8 

t^sa 

Self-absorption 

(l+Fr'/s 

(1+F)-V6 

9 


Characteristic 

(1 + F)° 

(l+F)-* 

10 


Self-absorption 

(1 + F)l 

(1+F)' 

11 

I'c 

Cooling 

(i+Fr2 

(1+F)“ 


Note that this formula for Y does not take the Klein-Nishina correc¬ 
tion into account. This frequency-dependent cori'ection is expected to be im¬ 
portant only at high frequencies, u > 10^^ Hz at f > 1 d (Fan & Piran 2006; 
Zhang et al. 2007). Upon detailed investigation, we find that a consequence 


of this effect is to reduce the overall energy required to match the light curves 
by up to 25%. The uncertainties aiising from model selection as well as 
due to correlations between parameters are usually also of this order and we 
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Figure 27. Posterior probability density functions for the physical parame¬ 
ters for GRB 120326A in the wind model from MCMC simulations. We have 
restricted £K,iso ,52 < 500, tc < 1/3, and eb < 


B. A WIND MODEL FOR GRB 120326A 

In Sections 5.1 we discussed the ISM model for 
GRB 120326A. We now consider the possibility of a wind¬ 
like circumburst environment. 

The best fit model in a wind environment requires p Ki 2.52, 
Ce « 4.9 X 10“^, cb w 1.0 X 10“^, A* Ri 7.0 x 10“^, Ek,™ ~ 

3.6 X lO^'^erg, fjet ~ 19d, Ay « 0.46mag, and Ei,_host.r' ~ 
1.6/iJy. This model transitions from fast cooling to slow 
cooling at 3.5 x 10“^d. The spectral break frequencies at 
1 d are located at « 6.8 x 10* Hz, k, 2.2 x 10*"^ Hz, and 
t'c Ri 4.2 x lO'^ Hz. The peak of the spectrum {f^) is at 
with a flux density of « 23 mJy. 

From our MCMC simulations, we And p = 2.52 ± 0.02, 
EK.iso = (3. Its) X 10^4erg. A* = (4.8!”) x lO”, Ce = 
(4.1![;*) X 10”, eB= (2.8t;4) x lO”, andfjet= 18.7!2td. Us¬ 
ing the relation 6 »jet = 0.17 forthejet opening 

angle (Chevalier & Li 2000), and the distributions of Ek.Iso, 
tiQ, and fjet from our MCMC simulations (Figure 27), we 
find 0jet = 2.it 3 degrees. Applying the beaming correction, 
= £-y,iso(l - cos0jet), we find E-, = (2.1 ± 0.3) x 10^^ erg. 
The beaming-corrected kinetic energy is much larger, Ek = 
(2.0t g) X 10^' erg. In this model, z/a falls below 10**’Hz at 

1.7 X 10“^ d and is therefore not probed by any of the radio 
data. Consequently, the model exhibits a degeneracy in its 
parameters due to to the unknown value of i/a (Figure 28). 
This high value of Ek also implies a low radiative efficiency, 

?7rad = E^/{Ey^+E^) Ki 1%. 

We now investigate the effect of energy injection in causing 
an X-ray/UV/optical re-brightening. The X-ray light curve 
during the re-brightening is located above the cooling fre¬ 
quency. In the wind model, the flux density above is 
oc (Granot & Sari 2002). For p = 2.5, 

this reduces to E,y>y_, oc Ej( |*Qf“’ **. During energy injec¬ 
tion, E oc f"', such that oc fi i3;w-i.38 'j’jjg steep jjse 

therefore do not consider this effect further in this work. 




Figure 28. Icr (red), 2cr (green), and 3cr (black) contours for correlations 
between the physical pai'ameters, A*, Ce, and for GRB 120326A, in 
the wind model from Monte Carlo simulations. We have restricted F^^iso ,52 < 
500, Ce < 1 / 3 , and cb < V3. The dashed grey lines indicate the expected 
relations between these parameters when is not fully constrained: £'K,iso oc 

^K.iso OC A* oc e^, Fk.Iso oc A* oc , and oc , 

normalized to pass through the highest-likelihood point (blue dot), while the 
dotted magenta lines indicate expected relations for changes the value of Uc'- 

^K.iso OC ^K.iso OC Ee, A* OC el, Fk.Iso oc Eg’^^, A* oc ef^^, and Ee oc 

Eg’^^. i/a falls below the radio band before any radio observations took place, 
and is therefore unconstrained, lies between the optical and X-ray bands 
and is better constrained; the correlations between the parameters along the 
lines of varying values of Vc are likely indicative of the strong changes (over 
two orders of magnitude) in the Compton y-parameter along these curves. See 
the on-line version of this Figure for additional plots of correlations between 
these parameters and p, fjet, 0jct, Fk, Ay, and . 


(ai = 0.85 ±0.19; Section 3.1) requires m = 2.0 ±0.2. How¬ 
ever 0 < ;« < 1 is bounded (Section 5.2.4), which implies that 
in this model, energy injection due to a distribution of ejecta 
energy to lower Lorentz factors can not cause the X-ray flux 
to rise with time. 

Relaxing the requirement m < 1, our best solution for 
the multi-wavelength re-brightening for an energy injection 
model requires two periods of energy injection. In the first 
episode between 1.7 x 10” d and 2.5 X 10 ^d, EK,iso increases 
as growing by a factor of 3.9 from 1.6 x 10^’ erg to 
6.2 X 10^’ erg. In the second episode, EK,iso « from 
2.5 X 10 ” d to 0.4d, further increasing by a factor of over 
2000 to its final value of EK.iso ~ 3.6 x lO^'^erg in this pe¬ 
riod. The resulting light curves, which are optimized to match 
the UV and optical re-brightening, cannot reproduce the X- 
ray light curve prior to the re-brightening (a shallower rise in 
the X-rays in the ISM model was achieved by placing be¬ 
tween the optical and X-rays). Due to the extremely large 
injected energy coupled with the fact that the steep rise in 
the optical violates the bounds on m, the wind model is a 
less attractive solution for the multi-wavelength afterglow of 
GRB 120326A. 

C. A WIND MODEL FOR GRB 100418 A 

We apply our MCMC analysis described in Section 6.1.2 
to explore afterglow models with a wind-like circumburst en¬ 
vironment for GRB 100418A. The parameters of the highest- 
likelihood model are p R 3 2.1, Ce ~ 0.33, r; 0.33, A* r 0.15, 
EK,iso ~ 4.0 X 10^’ erg, and E,,,host,white ~ 2.3 pJy, with neg¬ 
ligible extinction. This model transitions from fast to slow 
cooling at 1.5 d. The spectral break frequencies at 1 d are lo¬ 
cated at z/ac ~ 5.8GHz, z^sa ~ 32GHz, r 6.4 x lO” Hz, and 
z.'tn !R 3.6 X 10 ” Hz at 1 d, with Emax ~ 50mJy at r'c at 1 d and 
a Compton y-parameter of 0.6. 

Like in the ISM model, the optical and X-ray bands are 
located above both and This model does not require 
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Figure 29. X-ray, UV (top left), optical (top right), and radio (bottom left) light curves of GRB 120326A in the wind scenai'io, with the full afterglow model 
(solid lines), including energy injection before 0.4 d. 



Figure 30. Posterior probability density functions for the physical param¬ 
eters for GRB 100418A in a wind environment from MCMC simulations. 
We have restricted £’K,iso ,52 < 500, Ce < Vs, and cb < V^- The last panel 
corresponds to the flux density of the host galaxy in the SwijilWhite band. 

a jet break, and we find fjet > 140 d. Thus we cannot con¬ 
strain 0jet in this model, nor correct iso or EK.iso for beam¬ 
ing. The summary statistics from our MCMC simulations are 
p = 2.12 ±0.01, £k,so = (4.2 ±0.4) x lO^' erg. A, = 0.16“, 
Ce = 0.328 );oqq 9, and es = 0.31^oo2' We plot histograms of 
the posterior density functions for these parameters in Figure 
30 and present contours of the joint posterior density for the 
physical parameters A*, Ek.™, Ce. and cb in Figure 31. 

It is challenging to fit both the X-ray and optical light 
curves before the bump together in the energy injection see- 


Figure 31. la (red), 2cr (green), and 3cr (black) contours for correlations 
between the physical parameters, Fk^so. ^e, and eg for GRB 100418A, in 
the wind model from Monte Carlo simulations. We have restricted £’K,iso,52 < 
500, Ce < V^’ < V^- The highest-likelihood model is marked with a 

blue dot. See the on line version of this Figure for additional plots of correla¬ 
tions between these parameters and p and F,y host,white- 

nario under the wind model. Fike for the wind model for 
GRB 120326A (Appendix B), we find that the optical light 
curves before the peak require a steeper injection rate than 
allowed by a distribution of Forentz factors in the ejecta. In 
particular, our best energy injection model that matches the 
optical well but slightly under-predicts the X-ray data before 
0.5 d (Figure 32), requires E oc between 1.8 x 10“^ d and 

1.5 X 10“^d, steepening to E oc f* between 1.5 x 10“^d and 
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0.5 d. In this model, the blastwave kinetic energy increases 
by a factor of 4.4 between 1.8 x 10“^ d and 1.5 x 10“^ d, and 
another factor of « 325 between 1.5 x 10“^ d and 0.5 d, for an 
overall increase by a factor of 1440. Due to the large injection 
of energy required to account for the optical re-brightening 
in this model, we do not consider the wind environment to 
be a likely explanation for the multi-wavelength afterglow of 
GRB 100418A. 

D. A WIND MODEL FOR GRB 100901A 

We apply ourMCMC analysis described in Section 6.2.2 to 
explore afterglow models with a wind-like circumburst envi¬ 
ronment for GRB 100901 A. The parameters of the highest- 
likelihood model are p k, 2.02, Ce ~ 0.33, cb « 0.33, A* « 
1.7 X 10“^, EK.iso ~ 3.8 X 10^^, and Ay <0.1 mag. This model 
transitions from fast to slow cooling at 5.4 x 10< d. The spec¬ 
tral break frequencies at 1 d are located at r'a ~ 3.1 x 10*^ Hz, 
r'tn « 1.8 X 10'^ Hz, and r'c ~ 1.7 x 10^^ Hz, with k , 
3.8mJy at 1 d, and a Compton y-parameter of 0.6. Like in 
the ISM model, I'cis located between the optical and X-ray 
bands. Since the light curves in a wind environment decline 
faster than in the ISM case, the jet break in this model is later 
(fjet Ri 3.5 d) compared to the ISM model (fjet Ri 1 d). The jet 
opening angle is 0jet ~ 1.9, resulting in a beaming-corrected 
kinetic energy of 2.0 x 10 ^*^erg. 

The summary statistics from our MCMC simulations are 
p = 2.027to®^ £k,.so = (3.9 ±0.5) x lO^^erg, A, = x 
10-2, Ce = 0.32!°;°^, ee = 0.27!«;04, Ay < 0.1 mag, fjet = 3.5 ± 
0.2 d, 6 ljet= 1.9 ±0.1°, £k = (2.1 ±0.2) x lO^erg, and =. 
We plot histograms of the posterior density functions for these 
parameters in Figure 33 and present contours of the joint pos¬ 
terior density for the physical parameters A^,, LK.isoj 
in Figure 34. 

Like for GRBs 100418A and 120326A, it is challenging to 
fit both the X-ray and optical light curves before the bump to¬ 
gether in the energy injection scenario under the wind model. 
Like for the wind model forGRB 120326A (Appendix B) and 
GRB 100418A (Appendix C), the optical light curves before 
the peak require a steeper injection rate than allowed by a dis¬ 
tribution of Lorentz factors in the ejecta. In particular, our 
best energy injection model that matches the optical well but 
slightly under-predicts the X-ray data before 0.15d (Figure 
35), requires E oc t^ '' between 5 x lO-^ d and 0.12d, steepen¬ 
ing to £ oc ^ between 0.12d and 0.26 d. In this model, the 
blastwave kinetic energy increases by a factor of 9 between 
5x1 0-2 d and 0.12 d, and another factor of 15 between 0.12 d 
and 0.26 d, for an overall increase by a factor of Ri 140. Due 
to the discrepancy in the X-rays, the wind model may be con¬ 
sidered a marginally viable model for GRB 100901 A. 

E. A WIND MODEL FOR GRB 120404A 

We apply the methods described in Section 6.3.2 to explore 
afterglow models with a wind-like circumburst environment 
for GRB 120404A. The parameters of our highest-likelihood 
model are p ~ 2.03, Ce ~ 0.33, Ri 0.30, A* r;;! 1.9, LK.iso ~ 
1.1 X 10^2 erg, fjet R^ 8.9 x lO-^d, and Ay sa 0.12. This model 
also remains in spectrum 4 (Figure 17) for the duration of the 
observations, with the ordering < v^a- The spectral 

break frequencies are located at Ra 2.4 x 10'°Hz, r'sa ~ 
7.9 x 10" Hz, i/n,Ri4.0x 1022 Hz, and^c R^ 1.5 x 10'° Hz, at 
0.1 d with fniax = 7v,sa ~ 27 mJy. The Compton y-parameter is 


0 . 6 . 

The summary statistics (median and 68% credible inter¬ 
vals) of the posterior density for these parameters are p = 
2.04 ±0.01, Ce = 0.28![Jg^, cb = 0.18 ±0.09, A* = 2.3!^^ 
FK,iso = (l-2!!]i) x 10^2 erg, fjet = (9.2 ±0.6) x lO-^d, and 
Ay = 0.13 ±0.01, corresponding to a jet opening angle of 
0jet = 3°.0 ± 0°.3 and a beaming corrected kinetic energy of 
Fk = (1-6°52^) x 102 °erg, the model match the data after the 
optical peak well (Figure 36). We present the correlation con¬ 
tours between the physical parameters in Figure 37 and the 
marginalized distributions for individual parameters in Figure 
38. 

The light curves before 0.04 d in Figure 36 are based on the 
same energy injection model as presented in Section 6.3.3. 
In this model, the energy increases by a factor of Ri 27 like 
in the the ISM case. However, the injection rate is not com¬ 
patible with the maximum rate expected for a wind-like en¬ 
vironment, similar to the other wind models (Appendices 
B, C, and D). We note that the optical and X-ray frequen¬ 
cies are located above both the cooling frequency and 
with < r'm < r'opt < vx in this case, and in this regime, 
the light curves are independent of the circumburst density 
profile. Thus the measurements do not allow us to distin¬ 
guish between a wind or ISM-like environment in the case of 
GRB 120404A. Clear detections of a wind-like circumburst 
environment in conjunction with a steep energy injection rate 
in the future will enable us to furthur probe the massive ejecta 
model of energy injection in GRB afterglows, and thereby to 
further constrain the mechanism responsible for plateaus and 
re-brightening events in GRB afterglows. 
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Figure 32. X-ray (top left), UV (top right), optical (bottom left), and radio (bottom right) light curves of GRB 100418A in the wind scenaiio, with the full 
afterglow model (solid lines), including energy injection before 0.5 d. 
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Figure 35. X-ray and UV (top left), optical (top right), and radio (bottom left) light curves of GRB 100901A in the wind scenario, with the full afterglow model 
(solid lines), including energy injection before 0.26 d. 
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Figure 36. Similar to Figure 20, but for a wind-like circumburst environment. The light curves before 0.04 d are based on the same energy injection model 
described in Section 6.3.3 as applied to the case of the wind medium. See Appendix E for a discussion. 
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